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Die’ asymptotischen Entwicklungen für die beiden 
parabolischen Funktionen M,, „2(z) und W,, „2(z) bei großen 
Werten von » und z für -»<z/4x<+ © 
Von Herbert Buchholz in Darmstadt 


“Inder vorliegenden Arbeit werden die asymptolischen Entwicklungen der parabolischen Funktionen M,, „2 (2) 
und W,, „2 (2) bei großen Werten von x und z bis einschließlich zu den Gliedern 2. Ordnung in 1/x her- 


geleitet, und zwar mittels der Sattelpunktsmethode. Die Darstellung beschränkt sich auf den zunächst wich- 
tigsten Fall, daß das Verhältnis z]4x einer beliebigen positiv oder negativ reellen Zahl gleich ist. z und x 
selbst können jedoch komplex sein. 

The paper gives the asymplotical ewpansions of the parabolic functions M,, „2 (2) and 1 2/u (2) for 
large values of x and z up to the terms of the second order in 1/x. The derivation is accomplished by aid of the 
„Saitelpunktsmethode‘‘. The formula is limited io the case that is till now the most important, namely that 
Ihe ratio 2/4 x» is any positive or negative real quantily while z and x themselves may be complex. 

Dans ce travail les developpements asymptotiques sont donnes des fonetions paraboliques M, 2 (2) 
et W,, „a (z) pour grandes valeurs de x et de z jusgqu’ü les termes.du deuxieme ordre en I|x, et cela a l’aide 
de la „‚Sattelpunktsmethode‘‘. L’auteur se borne au cas, qui est d’abord le plus important, c’est que le quotient 
2/4 x soit reel (positif ou negatif). 2 et x pourtant peuvent ätre choisis complewes. : 

B 3Tof pa6oTe naeTca ACHMITOTHYECKOE Pas.lIoikeHue mapaboNnyeckux hyaeunü M 2 (z) 
an W,ur @) upu ÖOJIBIIHX 3HAYCHHAX % U z BEIIYaA NO WIEeHOB BTOPOTO TMOPAANKA OTHOCH- 
TeI5OoH 1/x. 3TO PA31oskeHHe IIPOHSBONHTCH IPH TOMOIM MeToNa „ruepbonmbeckux Toyer‘“ 
MIA HANDO.IeE BAKHTO CUYYAA, KOTAA OTHOMCHHE 2/4x PaBuaeTca 1000K NOANOKHTEJIBHOH HIIH 
OTPUHATEIBHOH BeINeCTBEHHOH BEIIHYHHE, TOTAA KAK CAMH 2 U x MOTYT ÖBITb U KOMINIEK-CHEIMH 
BEJIHYHHAMH. 


1.1. Allgemeines über die parabolischen Funktionen undüber 

-die hier benutzte Integraldarstellung. Die vorliegende Arbeit beschäftigt 
sich mit der Herleitung der asymptotischen Entwicklungen (abgek.: as. Entw.) der beiden para- 
bolischen Funktionen M,, .jg (2) und W, us (2) mit den Definitionsgl. (1) und (2). 
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konfluenten hypergeometrischen Funktionen bezeichnet. Da sie jedoch in den Anwendungen 
auf die klassische Physik der Wellenfortpflanzung i in einem auf die Koordinaten des Drehparabols 
bezogenen Raum die gleiche Rolle spielen wie die Zylinderfunktionen in einem auf Zylinder- 
koordinaten bezogenen Raum, so sollen sie in, dieser "Arbeit parabolische Funktionen genannt 
werden. 

Für das Folgende ist daran zu erinnern, daß zwar nach (2) W, „.»(@) eine in u gerade Funktion 
darstellt. Die beiden Funktionen M, „,(2) und M, _, (2) sind aber nach (1) im allgemeinen 
verschieden und bilden für ein au=m=0, +1, +2... einlinear unabhängiges Lösungssystem 
der D. Gl. (3) der parabolischen Funktionen. 
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Umgekehrt stehen zwar nach Gl. (1) die beiden Funktionen M,„ „p2(2) und M_, „(2 e* =i) in 
einem einfachen Zusammenhang miteinander. Die beiden entsprechenden W-Funktionen sind 
aber im Gegenteil die beiden einzigen Lösungen dieser D. Gl., die für alle Werte von 2, x und u 
voneinander linear unabhängig sind. Sie entsprechen in dieser Hinsicht den beiden Hankel- 
schen Funktionen, die sie gemäß Gl. (4) für x = 0. als Sonderfälle enthalten. 

Die gesuchte as. Entw. wird nach der Sattelpunktsmethode (abgek. Sp.-Methode) her- 
geleitet werden. Für eine möglichst ergiebige Anwendung dieses Verfahrens ist es wesentlich, 
eine geeignete Integraldarstellung zu benutzen. Ein maßgebender Gesichtspunkt für ihre Aus- 
wahl ist der, daß zum mindesten der schnell veränderliche Bestandteil des Integranden im Sp.- 
Integral, der den Verlauf der Paßwege bestimmt, eine eindeutige Funktion von möglichst ein- 
fachem Aufbau ist und daß sich die Integraldarstellungen für die verschiedenen Lösungen der 
D. Gl. (3) höchstens in ihrer Wegführung, nicht aber in ihrem Integranden voneinander unter- 
scheiden. Dabei bildet dann das Ideal, daß sich die maßgebenden Integrationswege bis ins Unend- 
liche erstrecken, weil in diesem Falle die freie Verschiebbarkeit komplexer Integrationswege am 
besten ausgenutzt werden kann. 

Nun läßt sich zwar im vorliegenden Falle diese letzte Forderung nicht im wörtlichen Sinne 
erfüllen, wohl aber dem Wesen der Sache nach. Benutzt man nämlich die durch Gl. (5) beschrie- 


Ps) =e?* 2 @2 Ras. (nie) tgy=Tgs,'tgs, mit $ = arc (Cpjs)a 22H) 
. > 3239 +43 . 
bene Funktion F(s) als Integranden, in dem allein die nme den schnell ver- 


änderlichen Bestandteil bildet, so ist es in der Tat möglich, die auf der linken Seite der 
Gl.(6a...k) 
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3 aufgeführten 1 ae alle durch die 
- Integrale I... X über ein und dieselbe . 
Funktion F(s) von Gl.(5) darzustellen, 
_ wenn man dementsprechend die Integration 
längs den in Abb. 1 mit I... X. bezeich- 
neten Wegen erstreckt. - Dabei. scheinen 

“ allerdings einige dieser Wege, weil sie 
allein zwischen den im Endlichen gelegenen 

En Punkten + ni/2 der s-Ebene verlaufen, 

2 mit der zuletzt aufgestellten Forderung 

nicht i im klang = zu stehen. Diese beiden 
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. are d daß, wenn nur die Tafeerakone ur 
wege i in ‚geeigneter Richtung in sie hinein- 
laufen, sie sich in nichts von dem Ver- 
halten der Wege in der Nachbarschaft des 
© unterscheiden. Das Integral 
ırde bereits in dem Bericht des 
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und (8) Gebrauch macht, die ihrerseits aus (1) und (2) folgen. Es versteht sich, daß nicht alle 
Integrale in (6) für die gleichen Werte von u und x konvergieren. Auch müssen bei einigen von 
ihnen zı und x bestimmte Forderungen erfüllen. So muß z.B. bei Gl. (6e) R (x — (1 +n)/2) < 0 
sein. , 

Im Hinblick auf die späteren Angaben ist noch ein Hinweis auf den Phasenwinkel von 
Coj s in der Potenz (Eof s)-!+#) erforderlich. Setzen wir darüber ein für allemalfest, daß für 
einen Punkt ; = $-+185 £ dem Horizontalstreifen zwischen den Punkten +r:/2 im Hinblick 
auf die Gl. (5) |are (Coj s) |< /2 ist, so liegt der Phasenwinkel von @ojs für einen Punkt 
s=s, 4%°5 Ar zip um zp höher, wenn s, > 0 ist” und um xp niedriger, wenn s; <0 ist. 
Liegt mithin der Integrationsweg wie im Falle der Gl. (6b) ganz auf einer Seite der imaginären 
s-Achse, so kann man auch für er 2 (2) eine Integraldarstellung angeben, in der das Integral 
zwischen den Punkten +r:/2-(2p—1) und +n:/2-(2p +1) verläuft. Sie hat die Gestalt 
der Gl. (9), wenn darin der a von Coj s dem im Halbstreifen zwischen diesen beiden 
Punkten für s, > 0 entspricht. 


+ @p+1) | ; 
2 ee { z 
Ws Narren F(s)-ds REF RREERNN N Fer (9). 
+7 07-0 P=0,+41,+2...) ; 


Verläuft ein Integrationsweg, wie es gemäß Abb. 1 vorkommt, teils in der linken, teils in 
der rechten s-Halbebene, so hat das solange nichts auf sich, als das Hinüberwechseln von der einen 
auf die andere Seite im Fundamentalstreifen —ri/2... +ri/2 stattfindet, weil sich in diesem _ 
Streifen der arc (Coj) stetig ändert. Innerhalb aller anderen Streifen ist jedoch mit diesem Über- 
tritt ein Phasensprung verknüpft, es sei denn, man trifft ausdrücklich die Festsetzung, daß bei 
einem solchen Übertritt, wie er z. B. beim Integral (IV) oberhalb des Punktes 4/2 von I;inks 
nach rechts erfolgt, die Änderung des Phasenwinkels von Coj sim Sinne der stetigen Änderung 
vor sich geht. In Abb. 1 wird an diese Vereinbarung durch Strichelung des bis dahin voll aus- 
gezogenen Linienzuges erinnert. 
Wegen der hier benutzten Integraldarstellung vergleiche man auch die Arbeit von Th. 
Sex1]2] 

12. Die, Lage’der-Sp., die as, Entiw. der sp. Integrale und die 
Sp.-Wege. In dem Integral, das hier zur Definition der parabolischen Funktionen benutzt wird, 
hat.der Exponent im schnell veränderlichen Bestandteil seines Integranden das Aussehen 
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Wir führen darin für den Quotienten der beiden großzahligen Parameter z und x die ständige _ 
Abkürzung 
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ein. Dann läßt sich damit die Gl. (la) in der Form der Gl. (1b) schreiben. | 
lx2|-f()=2]2-[s-A—- Te) — U]... :--- ee (1b), 
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Aus ihr folgt im Hinblick auf (1c), daß die Integrale der Gl. 1.1 (6 ) die re vielen Gitter- 
punkte der Gl. (3) zu Sp. haben. A 
s=ı(+y))=+y+tzik (en a (3. 


Wir werden die Punkte + y als die Hauptsp. bezeichnen a die übrigen ebenfalls paarweise. 3 
auftretenden Punkte als die Nebensp. j 


Wir bedürfen dann vor allen Dingen einer Betechinngsvornt für die Größe I = 4 2 = > 
bei gegebenen Werten von 2/4#x = & +4. Sie wird durch das aus af 
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(4a, b) dargestellt, das sich explizite durch die Gl. (4a, ß) auflösen läßt. Wir heben als be- 
sondere Folgerung aus dieser letzten Gleichung hervor: Es ist für n=0 und 


.—- os <?s0, 2.0=€<1 | 3. 1I2E< 0 


Die beiden Hauptsp. + y, die für y 0 getrennt liegen, fallen demnach zusammen, wenn 2/dx =1 
ist. Der Nullpunkt der s-Ebene ist dann ein Sp. 2. Ordnung. Ferner merken wir an, daß in den 
Sp. (3) der Realteil von f (s), der die Betragshöhe des schnell veränderlichen Bestandteils der Sp.- 
integrale 1.1 (6) bestimmt, sich aus der Gl. (5) berechnet. 


R{f(Ey taik)} = —2nk-sinv + sinv-(Cof2a-sin2P —2P) ] 
 .cosv- (Sin 2a - cos 28 — 2a) “=0; +1, +2...) J 
Die Nebensp. liegen also nur für ein reellwertiges x mit » =0 auf der gleichen Höhe wie die 


Hauptsp. Hingegen liegen die beiden Hauptsp. -+ y selbst auf gleicher Höhe, wenn die Gl. (6) 
erfüllt-ist. Das ist z.B. für =» =0 der Fall. 


cosv - [2% — Sin 2a :cos2ß]=sinv- [28 —Eoj2a-sin2ß] :... - (6). 


Wir bringen die Integrale 1.1 (6) auf die für die Anwendung der Sp.-Methode geeignete 
Form. Wir sehen dabei zunächst von jeglichem Faktor vor dem Integralzeichen ab und nehmen 
die Integrale zwischen den Punkten 8, und S,, über deren Lage noch Genaueres zu sagen sein 
wird. Dann läßt sich das auf das Hauptsp.-Paar -+-y bezogene Integral in der Form der Gl. 
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schreiben. Entwickelt man hierin f (s) —f (y) in der Umgebung von s = gemäß Gl. (8) nach 
dem Taylorschen Lehrsatz, 


Fo) —fWy) = Fa: an 2 -Eoj?y.Tg4r2y. (sy)? 
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so berechnen sich für die ersten fünf Koeffizienten b; dieser Entwicklung die Ausdrücke (8a). 
Auch hieraus folgt erneut, daß im Falle y=0in s=y ein Sp. 2. Ordnung vorliegt. 
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Die Kenntnis von b, nach Gl. (8a) erlaubt, eine genaue Aussage über den Neigungswinkel 
der Tangente an den Weg steilsten Gefälles, die sog. Fallinie, im Sp. + zu machen. Wird näm- 
lich gemäß GI. (9a) der Phasenwinkel von Tg y mit 7, bezeichnet, so stehen mithin zur Berech- 
nung von 7, die Formeln 
Sin 2&-+i-sin2ß 
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zur Verfügung. 
Man rechnet mit ihrer Hilfe leicht aus, daß dann z.B. für n=0 und 
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n=0 n=+7 für Bo; n=0 


ist. Andererseits ist in erster Näherung mit sy = o.+r):e'? nach den Gl. (8) und (8a) am 
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: {fa fly} = +21e| Zar a. erw turen, 

emnach gibt es am Sp. + y zwei Richtungen steilsten Gefälles, die mit der positiv reellen Achse) 
Prieas Au Winkel (+) und (+) von Gl. (10a, b) bilden und ganz entsprechend am Sp. — y die 


er. | 


ET ng PER 
Ze 


138 


gleichartigen Richtungen mit den Richtungssiakein o&und 0 von G .(10%,ß). 

dann, daß der Anfangs- und Endpunkt S,und 8, des Integrals (7 ) für den Sp. ae sowie de ent- 
sprechenden Integrals für den Sp. — Y jeweils i in denjenjgen Unterlandsektoren liegen, deren 
W inkelhalbierende im Punkt s = + y Dar den Richtungen @, und = übereinstimmen. 2 
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Für den langsam randerlichen Bestandteil im Integral (7) gelten für die ee a 
a; in der Entw. A ae 
h(s—y) = {&oj ( —y) +Tgy: Sin (—y)}-Utn ae 
= + BY) ta ey? +r... net 
nach Potenzen von s—y die Formeln ES BT Sk Ft 
o-1l m=-—MtmTy m=zl +): et y I) 8 
a Hy BSH a 


1 ER 
= EBEN LEHNEN Er HR Re) 


Die entsprechende Entw. am Sp. —y unterscheidet sich hiervon- nur in dem Vorzeichen von 5 = 

Damit sind wir nun auch bereits in der Lage, unter Zuhilfenahme des bekannten Formel" 
materials aus der Methode der Sp., wie es z.B. die Arbeit von O. Perron [6] bereitstellt, 
die as. Entw. für die beiden Integrale Y„(-+y) von Gl. (7) anzuschreiben. Zu diesem Zueil = - 
empfiehlt es sich aber, vordem noch die beiden en (12a, b) einzuführen. ee 
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Mit ihrer Hilfe lassen sich die as. Entw. der beiden Wen Ya 3 == N von a (7) durch die beic n 
-Gl. (13a, b) wiedergegeben. BR: 
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“ Um nun späterhin den. Zusammenhang zwischen 
Pshölschen Funktionen definierenden Integralen von 
noch kurz auf den Verlauf der Sp.-Wege auch in größ 
Das geschieht am übersichtlichsten mittels der beide 
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in denen s—=s; +6: s, und wie bisher » = arc (x), 22 
ist. Bekanntlich bestimmt sich dann der Verlauf der üh nd 
steilsten Gefälles oder steilsten N in Bet ' 
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Die Fallinien unter ihnen findet man dadurch heraus, daß für den Punkt s = o einer solchen 
Linie ® (6, +7) < 0 werden muß. Für den Sp. s = + y selbst ist die Gl. (15a) identisch erfüllt. 
Auch die Funktion ® (+y, +y) hat daselbst den Wert Null, was in diesem Falle bedeutet, daß 
in jedem der beiden Hauptsp. infolge unserer Normierung der Integrale Y, (+ y) der schnell 
veränderliche Bestandteil des Integranden auf der Betragsfläche 1 liegt. 


Wir entnehmen den Gl. (14a, b) die folgenden allgemeinen Gesetzmäßigkeiten: 1. Spiegelt 
man eine zum Sp. + gehörende Fallinie in ihrem ganzen Verlauf am Ursprung der s-Ebene, 
was auf eine Vertauschung der Vorzeichen von s, und s, hinausläuft, so geht die Kurve in die 
zum Sp.—y gehörende Steiglinie über, denn mit 7 (s, + y) ist auch Y (—s,— y) gleich Null, wäh- 
rend ® (s, +7) und © (— s, — y) entgegengesetzte Vorzeichen haben. 2. Spiegelt man dagegen 
eine zum Sp. -+y gehörende Fallinie in ihrem ganzen Verlauf an der positiv reellen Achse der 
s-Ebene, so geht sie in eine zum Sp. y gehörende Fallinie über mit 2 und x an Stelle der bisherigen 
Werte zund x. Damit steht die Aussage der beiden Gleichungspaare (10a, b) und (10x, 8) durch- 
aus im Einklang. Mittels der Gl. (14b) kann man auch leicht die Bedingungen formulieren, die 
erfüllt sein müssen, damit etwa auf einer Fallinie, die zum Sp. + gehört, auch zugleich der 
andere Hauptsp. —y oder einer der beiden Nebensp. +y E ri liegt. 


Mit Hilfe der Gl. (14b) läßt sich auch die Frage nach dem Verhalten einer Fallinie beim 
Einlaufen in eine der beiden wesentlich singulären Stellen +7.2/2 beantworten. Setzen wir näm- 
lich in (14a, b) s, =o°cosp und 3 = +r/2 +o:-sinp mit 0 (0, so ergibt sich in erster 
Näherung: 


DE 


5 (Coj2&+cos2ß) . (14a), 


olefaeen a 


ee +9: eir; +7) =, YP(—y, Ei tr.cosv— ze -(60j2%+cos2ß) . (14ß) 
Das Einlaufen des Integrationsweges in die beiden wesentlich singulären Punkte + x i/2 muß also 
nach (14a) in Rücksicht auf die Konvergenz des Integrals (7) unter einem Winkel p gegen die 
Horizoutale erfolgen, der dem-Winkelbereich —-a/? 2? + <p< +n/2 +{£ angehört. Tut er 
dies, so konvergiert das Integral (7) auch dann, wenn die Beifunktion (Coj s)> +” mit R(u)> —1 
dabei über alle Grenzen wächst. Nun garantiert aber die Vorschrift (15a) ganz von selbst dafür, 
daß das Einlaufen des Integrationsweges in die singuläre Stelle in dieser Weise erfolgt, denn nach 
Gl. (14ß) erfordert sie, daß dabei ? —=£& -+0:const. ist, wobei sich dann die Konstante aus der 
Bedingung 
2 Y(—y7, +nta: cos v + Const. (Coj 2 DIESE) Y=09 
bestimmt. Der Einlaufwinkel @ in eine der beiden wesentlich stnenlären Stellen ist also auf einer 
Fallinie genau gleich dem Phasenwinkel von 2. Um diese besondere ausgezeichnete Einlaufs- 
richtung in den Bildern hervorzuheben, sind darin senkrecht zu ihr durch die beiden singulären 
Punkte gehende, kurze Geradenstücke EIESEEI BEL worden, die auf der Seite des Hineinlaufens 
‚eine Schraffur tragen. 
| ‚Die ins Unendliche ET anTERden Fallinien verhalten sich wie Gerade, die unter dem Winkel 
ek} ngegen die positive Richtung der reellen s-Achse geneigt sind. 


BE Die asymptötischen Entwicklungen für dieparabolischen 
Funktionen für Werte von 2/4x im Intervall— o<z/4x<-+ mw. Die 

bisherigen Ergehnisse setzen uns in den Stand, in dem praktisch besonders wichtigen Falle 
BD - re en Verhältnisses 2 2/4x den. noch fehlenden Zusammenhang ‚zwischen den 


he Dabei wird wohl gemerkt a die Reellität des Verhältnisses von 2/4 x voraus- 
gesetzt. Die Phasenwinkel von z2 und x brauchen nur der Bedingung zu genügen, daß ihre 
Mierkuz AR Null oder gleich +rist. Wir unterscheiden die drei Fälle — © 2 ar “<h, 


2 Be diesem Bereich der Sp. a2 Y auf de positiv reellen Achse der s-Ebene an der Stelle 
-& mi mit 0 < & S: ©. Der Verlauf der Fall- und Steiglinien, der nach Gl. E 2 (15a) nicht bloß 


4 Die Fallinie über den Sp. +r ist voll ausgezogen, die über a Sp.. _ a 
ktierte Linie BpeeERIgE 


u ye 
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Aus Abb. Ic ist zu ersehen, daß die Fallinie über den Sp. +& vollkommen: ne Wegführung 
des Integrals in Gl. 1.1 (6b) für die Funktion W,, „., (2) entspricht. Das trifft auch noch bei der 


zu Fall 1a) gehörenden Fallinie zu. Im Fall1b) acht sich jedoch die Besonderheit bemerkbar, 
daß die Fallinie über den Sp. + in ihrem oberen Verlauf auf den Sp. +& + zei aufstößt, für . ° 
den sie wegen der nied- | 
riger liegenden Betrags- 
fläche einen Kammweg. 
oder eine Steiglinie dar- 
stellt. Würde man sie 
über diesen Punkt hinaus 
ihren Wegaufderanderen 
Seite fortsetzen lassen, so 
würde dabei 
Rt) fa t+ml} 
wieder ansteigen. Deshalb 
muß von nun an die Fall- RR 
linie in bezug auf den Sp. 
+x-+zi weiter verfolgt 
- werden, und wenndies auf 
demlinkenTeil geschieht, 
so gelangt auch in diesem 


sm [am [ocH Tacn Fall der Integrationsweg 
I zum wesentlich singu- 
+ n/2 


#420 |. rich, + a2 ) + 2, . lären Punkt +:/2. Da 
—3n/d +n — n/4 37 . „ 
Tara: | 2808: [dw Hemer | me Krale aber die Betragsfläche 


und außerdem in Bild 2b: R{f(—a + ri)} = —2n. ‚dieses zweiten Passes auf 

geringerer Höhe liegt, so 

kann sein Beitrag zum Integralwert außer Betracht bleiben. Für die Winkel = v = —.n/4 

oder —r/2 liegen die Verhältnisse mutatis mutandis auf Grund der u: 2) aus. 
Gl. (15a) ganz analog. 

Es braucht mithin der Ausdruck 1.2 (13a) mit y = «& gemäß Gl. 1. 3 (6b) ı nur noch mit dem 


Bil a ® > 2 le Zu a me in 2 


0 


Re Te 3 Ser nl ul get ren Luz ans ma ie ae eo. 


Faktor 
2 we .. E 
PSP 2) = r 
mb : r (4) 
Br, rs Ze € re re) No Par a 
24 » Aei ER ERBET: ET er 


multipliziert zu werden, um die as. Entw. für W, „„(2) zu hen Dabei empfiehlt es sich, den 


dritten Summanden im Argument der Funktion exp als besonderen Faktor abzuspalten, ihn in 

-eine Reihe zu entwickeln und diese mit der as. Reihe im Ausdruck für Y,(+y) zu vereinigen. 

” Dann tritt an die Stelle des Koeffizienten P, (y;a) von Gl.1.2 (12b) der ar ie % aa 
der Gl. 


Gta’ 
ty) = Ds 
+6: (121 — 282) Tgty — 924: Tgry + 385} = - P,(W; en 


und die endgültige as. Entw. für W, „,j (2) nimmt damit die Gestalt der Gl. &) an. a ist glich. 
mäßig in dem ganzen Winkelbereich |&»]<r/2 eng: ER BE 
W > (2), = 5 i2 zu a ., x. 


{9 (4 — 9) Au —1)- 27 +36 rs) Er y Fe 
y 


27 rl en . 

=W „(4x Coj?a) = (Tg) ? .exp Ix- I (&) +2 ena]+- 
2 

1 Perg les, = SLaRR. en 2) U <ldn Copa < c 


An ihr ist besonders ee daß auch die rechtsatehendstas 
die angeschriebenen Entwicklungsglieder zu beurteilen ‚gestatte 


selbst gegen einen Vorzeichenwechsel von u unempfind 


# Zr %, 


a de 15 Te r ArBREE DE ZEN 
ie 
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Die über den Sp. s = — y hinwegziehenden Fallinien verlaufen in den Fällen a und e von 


Bild 2 wie der zum Integral 1.2 (6h) gehörige Integrationsweg. Multiplizieren wir also 
die Gl. 1.2 (13b) mit dem Faktor 


14 
% 2 e+ri* 
24 . 
T (* er! (3) 
ar), 
I+u a 1 ) 
2 +rix— (#44) 10% Br: ESS | 
2 3 Al ) N 
aan? a Te a 0) 


so stellt die entstehende Entw. mit vorläufiger Ausnahme des Wertes v = + /2 as. die Funktion 

N ua(z e-”‘) dar. Dabei muß in dem obigen Faktor das Zeichen — aufgenommen werden, 
weilsich nach Abb. 2im Integral 1.1 (6h) der Integrationsweg von — oo nach dem Punkt —rz u/2 
erstreckt, der Aufstellung der Gl. 1.2 (13b) aber die Annahme zugrunde liegt, daß der Sp. in 
der Richtung > — ©(7 überschritten wird. Was nun den Fall 1b) mit v=-+-x/2 anbetrifft, 
so führt er, wie ein Vergleich des Fallinienverlaufs mit demWeg des Integrals in Gl. 1.1 (6d) lehrt, 
zwar ebenfalls auf die Funktion W (2° e”"t), jedoch ist diesmal mit dem Faktor 


— 2%, u]2 


} a 1 1 . . . 
Eher ye2 H Rn et , -rix+Z (4m 
(2.6192) RE ETF ie 


RE 2 Rü : ‘ et. 
R am P7 1 / 
z Ks r(r 22) 2.cosaz (# +6) 


zu multiplizieren. Da aber wegen &(x) > 0 und |#|> ©2-cos (nr (x +u/2)) » exp 
(—rix —ni u/2) ist, so sind in den drei Fällen la, b,c die obigen Faktoren gleich, und es 
gilt daher für W_, „„ (2° e””') die as. Entw. (4a) mit dem charakteristischen Winkelbereich 


0=zL,v<-+n/2. 


s a — ni 
W_ um ei W_, nt Coftaom)o(üga) "exp = 1) 
1 
ee : (4a) 
En 3 Q, e: ER  Qeta; u) EN 
Be Sin 2u| — —— I. u + 2 +0 3) 
"ld <azlax= &a<0;Ct=1;0ZLv<H+B) 
Um den Verlauf der Integrationswege in Bild 2 für die Winkel £, v = —/4 oder —ın/2 


überblicken zu können, braucht man sich lediglich die Wege in Bild 2 b, ce für &,v» = +n/4 oder 
+/2 an der reellen s,-Achse gespiegelt zu denken. Das Integral 1.2 (13b) stellt dann jetzt die 
Funktion W__, „(2° e*”') dar, und durch ganz gleichartige Überlegungen wie oben gelangt man 


dadurch für diese Funktion zu der as. Entw. (4b) mit dem charakteristischen Winkelbereich 


—n2=. 569 <0. 
» 7 - a | 
2 2 
1 
0; 12) ‚9216; (4b) 
u?) 1\] 
4 a eye +0(3)1 


1<24xr=bofa< w;i=v; —FZ Lv <0) 


Im Falle £,v = 0 hat man sich in der zu diesem Fall gehörenden Abb. 2a die Fallinie über den 
Sp.—a durch die obere Hälfte der über den Sp. +& zum N 8 gehenden Fallinie er- 
 gänzt zu denken. 


 —  Bildet man mit den Gl. (3) und (4a, b) die Gl. 1.1 (8) nach, so ergibt sich für die Funktion 
iM zZ eu die as. une - 
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ö Fr 
M „@ M „(4x:boj?o) Fi ur 
% 2 43 2 en - Sa 48 x v x 
Je 2%-% & .eo ee . y N 
Tut) u) u FE RE 
5 < sn Tr Iru_ ): + x (&i g Q% Q, | Er Bi 
2 sin 2 * = = * 'r x Se “e u 
REN) ee er 
Fair 5 =) - = 2° N 
Qı2 = Qı,. We Re P E% ; RK ee iu “ iS 
(1<2/Ax =bofo< 0; a ob. Vorz.: 0 (£, a <s+t3_ 5 
unt. Vorz.: —Z=(&,»)s0) ) © = 


Gemäß den dazu gemachten Angaben tritt das Doppelvorzeichen in (5) nur im Falle =v = 0 SE 

auf. Jedoch ist die darin liegende Zweideutigkeit nur scheinbar und ein Beispiel für das bekannte 
Phänomen von G. Stokes. Da nämlich für jedes «> 0 auch &in 2& > 2x ist, so ist für = 
x > 0 der zweite Summand in der geschweiften Klammer vernachlässigbar klein gegenüber dm 
ersten. Solange also*bei reellenWerten von «und z das Verhältnis 2/4 x > list, hat die Funktion 
M,„uj2 (z) auf Grund der Gl. (5) gewiß keine Nullstelle. Die beiden Summanden in dieser Klam- 3 
mer sind überhaupt nur für ein % (x) =0 von gleicher Größenordnung. Im Grenzfall z=i: 15 38 
mit 7 > 0 ergibt sich z. B. im besonderen aus (5) die as. Entw. 7378 


M „(-4rbo%o) mi M'7,e- Ar Vopa)in 102 a 
Der ee Bi > ri | R 
P(il+R) Ti+Rn) E 
9 18 oh sn | . ? 
(u) Zr, op[4ar+ = H ne Az N: ; 


. lcos [r(&in2«— 24) Si = | -% +0 9) 


It 


+ sin E (Sin2&—2u) ik (24 + 0 «)) 


ar 7>0; =.) Br 
Aus ihr kann geschlossen werden, daß die links stehende Funktion bei festen Has von & un- ; 
endlich viele Nullstellen in bezug auf 7 hat. Sie sind identisch mit den ‚Eigenwerten von Mi: iz le 
(* const) hinsichtlich r. 


Ein anderer interessanter Grenzfall von (5) liegt vor, wenn em: ES a er WR ist . mit "ri 
n=0,1,2... Dann entartet nämlich die Funktion 4 


„Hl z Sl Gl. ACER in das ag er 
Polynom io ET a: 3 
# (x an) Erle M, n m a ER 2 
I (2) a er eh zZ 2 Uyez Re Er r zu. e ( BR, a Diem 
IR £ Ti+n 1 SUR DE 1 
| r(r ER SE) | d+m) BI 
Es hat in dieser Definition im Hinblick auf die Gl. 1. 1 a) als Glied mi öch 


z den Summanden (—2)"/n! Bei großen Werten v von n N 2: = Od e 


tive ganze Zahlist.. Auf Grund von ae ist nach Gl. 1 


Polynom as. gleich dem Ausdruck. er 
AN er re a 
N ma Tga BB Bir ur Be 
P.(6; ER : 
+06 (on 


Geht hierin 2/4 x > ©, so wird angenähert 24 In 
aus (8) den besonderen as. REDEN 
E£ I+ ER 


ren ee. 
Für sehr großes n entspricht er dem Gli 


=v=+n/2. Der, en Sg, 
für &,v=—n/4 und 
x ebenfalls klar. un 
_ Unterschied in dem = 
Verlauf der Fallinien 


von Abb. 3 gegenüber 
- Abb. 2 besteht in dem 
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u R< 1. Für diesen wertiereich liegen die beiden Sp. + y auf der 
imaginären Achse der s-Ebene, ‚ und zwar für 2/4% = 1 in der Nähe des Nullpunktes und für ein 
z/4x =0 in der Nähe der beiden wesentlich singulären Stellen + 2/2. Der Verlauf der Fallinien 
istin Bild 3 wiedergegeben, und zwar wieder unterschieden nach den drei Unterfälln =v=0, 
C=v=H+nl und 7. 


Verlauf-der Fallinien 
= — r/2 ist am 


Der auffälligste 


Fehlen einer unmittel- 
bar vom Punkt —:/2 


2 i -3 
zum Punkt + ne/2 
übergehendenFallinie. 
Vielmehr zieht sich j a N b o 
hierindenFällen3a,c Bild3. Der Verlauf der Fallinien durch die Sattelpunkte s = + iß für ein 2/4# = +0,5. 
die vom Punkt —i/2 — - 
- : : Es ist in Bild =» Rlft+iB)| (+ Y) | wi (+ >> —y) 
ausgehende Fall- - s - are 
linie sofort über den >88 BR a 1 U Hal — a +3 a 
. . . B 7 70, —n — n/2 + .nj2 
Sp. — iß in die linke 3c il | #0404 | a8 | +m8 | ans | +5 2/8 
3 s-Halbebene hinüber - und Außerdem in Bild 3b: Ri (+iß+ zi)} = + 0,571 —2r 
und läuft dort in zu- und RICH iB—aU} = — 0,571 +4 2 
nehmendgeraderRich- - 5 
tung unter dem Winkelgo = —v» +rz mit v=0 oder r/& gegen den Punkt ©. In der 


gleichen Weise vollzieht sich der Verlauf der vom Punkt cı/2 ausgehenden Fallinie über den 


Sp. +iß. Der Zusammenhang mit dem früheren Verhalten der Fallinien wird aber sofort ver- 


‚ständlich, wenn man die Sache so ansieht, als handle es sich auch hier um einen einzigen vom 
Punkt — ri/2 zum Punkt + r«i/2 verlaufenden Weg, der jedoch unter den besonderen hier vor- 
- liegenden Bedingungen den Übergang von dem ersten zum zweiten singulären Punkt nicht direkt 


vollziehen kann, sondern wegen einer dort vorgelagerten hohen zaliche nach links aus- 


weichen muß. 


Die Fallinien zu Bild 3a, c en nun offenbar, solange Mansie äls neoht zusammenhängend 


s ansieht, als die Integrationswege der Integrale 1.1(6g, h) aufgefaßt werden. Das führt dann 
mittels der Formeln 2.1 (3°) für die Funktion W_ 


ER e e*"ö) sofort zu den beiden as. Entwick- 


lungen. en 


1 


Ma er) W__u(4r cost Bet) (tgß) ” 


| 1 


52: ir 
F RZ . = x 1 
a. ze E. Er Fr SEHE u ir 
AESEREN resp v2: Ep-ein2n + Tr - In Be . . (la), 
> ps ECHESIZ 1) ) 4 Se 10( «3 
BR (Sp: on | 
t e a 
RE 1. AFS URN 
| e ) 9: 8.%. il 


Bes +8 üd=itgh ’ al 


E; 
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% en IR} en 
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die damit bereits für den Winkelbereich | £, » | < x/2 als gültig nachgewiesen sind. Was nun den 
Fall von Abb. 3b anbetrifft, so stellt hier offenbar die Fallinie über den unteren Sp. —:ß, die 
in den oberen Sp. +ß von unten her kommend als Steiglinie einläuft und dann auf der Fall- 
linie dieses Passes den singulären Punkt +7 i/2 in einem nach rechts: offenen Halbkreis umläuft, 
gemäß Gl. 1.1 (6d) ebenfalls die Funktion W__ BES CH e”‘) dar. Von den drei Pässen, über die 
dieser Integrationsweg nach Abb. 3b verläuft, kommt jedoch für die Auswertung nur der Beitrag des ; 
untersten in Betracht. Die zugehörige as. Entw. entspricht dann aber wieder genau der Gl. (1b). 7 
Andererseits stimmt der Weg über den oberen Sp. in Bild 3b im Hinblick auf die Gl. 1.1 (6k) 
mit dem der Funktion M,, BE) überein. Er kann aber auch, da der Beitrag des Sp. — ip + ni 
gegenüber dem des Sp. +Bi vernachlässigt werden kann, mit dem Weg des Integrals 1. 1 (68) 
in Beziehung gesetzt werden. Der Geltungsbereich von (la) und (1b) hat sich damit im Hinblick 
auf die frühere Bemerkung über die Spiegelung der Wege an der reellen Achse auch auf die beiden 
Grenzwinkel £, 9» = +r/2 ausdehnen lassen. 
Zieht man nun noch die beiden Gl. 1.1 (7a, b) heran und berücksichtigt die as. Formel, 


1+ ee | j | 
te).r(a+ Ze) 2n:exp 12%- © 3 al 
r(* _ ): n(r -- 5 ) 2n-exp 3 2x-In 3 + 1% +0 (x) - Är, 
((x| > ©; Jarc (=)| < =) | 
so ergeben sich aus (1a, b) für W, „„(2) und M,, .jz (2) sofort die beiden weiteren as. Entw. 


a 


1 
5 % W—z 
en (2) = RS, cos? ß) ud: (tg P) - eXp “in (*) + 3% +0 (273) 
RR. (8% 
sin D (28 —sin2ß) En (i we #9) 
+ cos *26—sin2) +2] (Gi ee) R 
M» 2.121. M> 7, (4% 00528) 1 u 
ar »7 RÄDER LEER | AN wege A i 
= PR , - ——— 0 a 
Ti+n Titan e u: | Be 


‚sin [x (a —2ß-+5sin2ß) — 2 (n- 3): (1+ 0 ( — a RT LA). 


—cos |x( (a —2ß + sin 2 B)— "(n— ): (: aut in +0 (7 »)) 


0<2dx=cstß<1; C=v; O<PB<Y; —F<Z 5: <+?2) 

Bei einer imaginären Komponente von x darf natürlich in diesen Beziehungen der sin und cos 
gemäß der Eulerschen Formel durch eine der beiden Exponentialfunktionen ersetzt werden. 
Ohne die Glieder in x! und x? wurde die Formel (4) auf anderem Wege zum ersten Male von 
A.Erdelyi [4] angegeben. 

Mittels der Gl. 2.1 (7) kann die obige Gl. (4) wieder dazu benutzt werden, ı um auch für das 
Laguerre-Polynom im Falle großer Werte vonnundz2=&> O0 einen as. Näherungsausdruck 
zu gewinnen. Es ergibt sich diesmal dafür die Gl. ; 


- + u(®—1) 


(It en a2. Li (a) m (2: cosß)=r (mir sin2p) i u) 


sin |» (28 —sin2ß) + +2» (1 yelniron ;) 
+ cos [+(28—sin 29 “ {i UBER Loge) 


(<= osB<1; En ureell; E=v=0; 0<ß<z ar 


Es versteht sich, daß es im Hinblick auf diesen Fall allein zweckmäßiger wäre, nicht mit dem hie 

benutzten Hilfswinkel ß zu arbeiten, sondern etwa mit dem durch die Gleichung cos® P' = x/4 
definierten Winkel 5’. Man vergleiche im übrigen mit der Bor (öyei eine Be al de E 
Arbeit [5] von F. Tricomi., } ER | 


wie a A et a 47 nn a A EB s 
ne ea - Re Da er e BAR . 
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Aus der Gl. (4) lassen sich weitere bemerkenswerte Schlußfolgerungen ziehen. Sind x und 
2 = «rein reell und positiv und genügt das Verhältnis @/4 x der zu Gl. (4) angegebenen Forde- 
rung, so hat die Funktion M, „, (x) nach (4) offensichtlich unendlich viele Nullstellen in bezug 
auf x. Bezeichnen wir die als groß vorausgesetzte n-te Nullstelle dieser Funktion mit Ku,n, SO 
läßtsich aus der Gl. (4) ohne viel Mühe die asymptotische Formel (6) herleiten. 


1 4 
= "(n +54) ip SA —1)-tß +6 +5, 1 
FETTE ET ERTEENNE PER : 73 
x —2ß+sin2ß: 9%  (n +44) en +82) 6). 
4 Pa! 
(n groß ganzz., co®?ß=z2/4x=const, O<P<E) 
Ist andererseits in (4) im besonderen «= — i: tmitr> Oundz=— ii: nmitn> 0, so 
ist in allererster Näherung 
Re) ; = Be Sue en 1 
MEIN nn AR ER 29 RES le 
Tut) & (6) DT e ein |r( (7 ß-+sin2P) An 5/\ 7). 
(u reell) 


Daraus kann geschlossen werden, daß die Funktion M_, zug (in) bei einem.,n > 0 für alle 


t > n/4 gewiß nullstellenfrei ist. Dahingegen besitzt bekanntlich die Funktion De un.) 
mit 7 < 0 unendlich viele reelle Nullstellen in bezug auf r. 


2.3. Der Fall— © <z/4x<0. Ist das Verhältnis 2/4x gleich einer beliebigen nega- 
tiven Zahl, so liegen die beiden Hauptsp. +y in der s-Ebene an den Stellen + (& + 1/2) mit 
0<a< o. Diein den Gl. 1.2 (13a, b) auftretenden Hyperbelfunktionen von y stellen sich dann 
als Funktionen von & gemäß Gl. (1) dar. Die Differenz der beiden Phasenwinkel & und » von zist 


Siny=s3-&oja...(da) Lofy=i-Eina...(b) Toy=Üsa..... (1e) 
Gy=Toa. . . (ld) -S&in2y—2y=—- 2a +6Smn2a+tnmi .... (le) 
bei dieser Lage des Sp. +y im Hin- 
blick auf die Gl.1.2(2) durch dieGl.(2a) 
gegeben. Hätte man y=oa—ni/2 
gesetzt, so wäre C — = —rn zu 
setzen. -Wir bleiben jedoch im folgen- 
den bei der ersten Festsetzung. 


N EEE WERE ET (2a) 
a rare Haar ea ee (2b) 
ED RE EEE (2c) 
WEN gr en he te Zen el (2d) 


Die Frage nach dem Zusammen- 
hang der beiden Sp.-Integrale 1.2 
(13a, b) mit den Integraldarstellungen 
1.1 (6) für die parabolischen Funk- 
tionen entscheiden wir auch im vor- 
liegenden Falle auf anschaulichem 
Wege. Zu diesem Zweck sind wieder 
in Bild 4 für einen bestimmten, aber 
- willkürlich gewählten Zahlenwert des 
Verhältnisses 2/4 x getrennt nach den 


ild4. Der Verlauf der Fallinien durch die Sattelpunkte s = 
drei Unterfällen eines & = 0, 1 [4 +(&+ ni/2) für ein 2/4%x = — 2,019. Die Zahlenangaben für 


- ° En, Ä einzelne Punkte der imaginären Achse sind mit © multipliziert 
und + z/2 die Fallinien durch die zu denken. 


beiden Hauptsp.aufgezeichnetworden. zsıstm Bild 
Aus dem Verlauf der zum Sp. +y ge- „771 


ör ini i : 4a an +n/2 + 32/2 
hörenden Fallinie auf den Bildern 4a, c 4 met A alikaak 
geht dann sofort hervor, daß das Sp.- _ 4e —z/8 |+ 7 nj8 | +30/8 |+ 11218 


' Integral 1.2 (13a) gemäß Gl. 1.1 (6e) 
die Funktion W, „2 (@) mit entgegengesetztem Vorzeichen darstellt. Im Bild 4b entspricht der 
Verlauf der Fallinie der Wegführung des Integrals 1.1 (6b), das mit dem gleichen Vorzeichen 


ebenfalls ‚der Funktion W, ‚up (2) entspricht. ‚Jedoch ist es hierin mit einem anderen Faktor 
geechen: Nun ist aber das Verhältnis dieser Faktoren 


oe 
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und da im Falle von Bild 4b $ (x) <O und |x | sehr groß ist, so sind die beiden Faktoren. asymp- 
totisch einander gleich. Wir könnten demnach sofort die as. Entw. von Wu (@) anschreiben, 


und nach den eben gemachten Feststellungen hat sie im Hinblick auf die Spiegelungsregel i in 32 
bezug auf » beieinem & —» = +.den charakteristischen Winkelbereich —-3r2Z» s—n]. 
Für die praktischen Belange ist es aber wesentlich bequemer, den vorderen Parameter x allemal A 
auf den Winkelbereich | are (x) | <#/2 zu beschränken und in einem Falle wie dem vorliegenden Se 
dann lieber dem vorderen Parameter ein negatives Vorzeichen zu geben, haben wir ja doch auch 


er 
schon früher neben der Funktion W, „‚ (2) regelmäßig die Funktion W_ 2,2 (2) betrachtet. Wir = = gi 


tun daher besser daran, in der hier anzuschreibenden as. Entw., den Parameter x gemäß Gl. (2b) 


durch x’ - e”* zu setzen, und erhalten dann: : EA 
. # J A - 


2 1 R Be 
Won) W nr Sina) (Tg) ” = 
2 7 2 = 

ae - . Vest se 

-exp?— x In (2) x +20 + ein er (3) 3 
en 

EZ + 0 Sl 


0(<a<o;, =; —ZZzV SH) 


Etwas schwieriger gestaltet sich die Entscheidung der Frage, wie der Verlauf der über den 
SA VENEN 1/2 hinwegziehenden Fallinie auf den Bildern 4a, b, c mit den Integrations- 
wegen in den Gl. 1.1 (6) in Übereinstimmung gebracht werden kann. Am schnellsten dürfte hier 
die "folgende Auslegung zum Ziele führen. Wählen wir den Gang der Fallinie in Bild4c durch 
den Sp. — y zur Grundlage eines Vergleichs, so stellt offen- 
bar der Verlauf der zu dem gleichen Sp. gehörenden Fallinien 
auf den Bildern 4a,b je einen Grenzfall zu dem Fallinien- 2 
verlauf in Bild 4e dar. In Bild4a kann man nämlich vn 
dem Verlauf der Fallinie auch sagen, sie: 'entferne sich, nach- ar 
dem sie als Kammlinie von oben und unten her in den 
Sp. + — i/2 eingelaufen ist, auf einem doppelt belegten & 
Weg in horizontaler Richtung nach rechts ins Unendliche. In = 
Bild 4b kann man den Verlauf so deuten, als zöge sich hier 
der Integrationsweg nach dem im Punkte s= — 3z/2 er- 
folgten Übertritt von der linken zur rechten s-Halbebene i in. 
Richtung —t. © ins Unendliche. ‚Das ‚könnte lich .B. in 


Nebensp. + —p-° zil2 ‚mit p= ll 
alle diese Sp. haben in der rechten ee niedrig 
Bild5. Bei stetiger Abnahme des Phasen- Betragsfläche als ‚der Sp — & —nij2, | so daß ee 


winkels von z drehtsich die schraffierte Ge- gehörigen as. Entw. infinitär kleiner 
rade im Punkt s = — ri/2 von Bild1 im a $ 
Bamelgersinn. Een en ahme von arc (2) erwähnten Sp. Wir, 'beh 
um den Winkel 2” verläuft dann der We 

(VI) in Bild 1 wie in diesem Bilde Pe Bild 4a, b, e gleichzeitig, i 
geben. parabolische Funktion aueh das ı 


Nun befindet sich ein Integral mit dieser Wegtührangt nicht ur nitt 
wiedergegebenen Integrationswegen. Es ist aber darin:enthalten 
nutzen, in dem Integral der Gl. 1.1(6f) den Phasenwinkel 
. |&4 < n/2 heraus stetig zu verringern, bis er um 2% BES yunE 
das schraffierte Geradenstück ebenfalls im Uhrzeigersinn u 
Integral 1.1 (6f) gehörende Integrationsweg hat dann gen 
lichen Verlauf wie in Bild 4c. Dieser neue Integratic 
— x — 1/2 geführt werden. Im tlinblick auf die Gl. 
. die Darstellung 


n 


-", © + arc(-v—n) 


ur za m) 


E ers 22 28, ds 


a w@- ee) — 


= 2 


"in der der Integrationsweg dem Bild 5 entspricht. Nun darf aber nicht der Einfluß des Phasen- 
winkels von &oj s übersehen werden. In (4a) hat &oj sin der Halbebene s, > O in dem horizon- 
 talen Streifen zwischen den Punkten — 1/2 und — 3 mi/2 einen Phasenwinkel y aus dem Bereich 
an <yp< — a2. Auf dem Teil des Integrationsweges, der der Halbebene s, < 0 angehört, 
- hat dann Cof s einen Phasenwinkel, der sich nach dem Überschreiten der imaginären s-Achse 
unterhalb des Punktes —rij2 aus dem eben angegebenen durch stetige Fortsetzung ergibt. Auf 
dem Geradenstück « sı<0, gs, = —r/2 ist gemäß dieser Vereinbarung in (4a) arc (Coj s) = 
— 37/2. In dem Sp.-Integral von Gl. 1.2 (13b) hat aber Coj s im Punkte — «— i/2 den Phasen- 
‚winkel +7/2. Wünschen wir also für die Funktion auf der linken Seite von (4a) eine Integral- 
darstellung zu haben, in der der Phasenwinkel von &of sim Sp. — x — 1/2 diesen Wert +/2 
hat, so Ne in (4a) Eoj s-exp (— 2ri) an Stelle von Eoj s treten, so daß jetzt die Gl. (4b) gilt. 
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Ih (4b) läßt sich nun das Integral auf der rechten Seite ohne weiteres durch den asympto- 
- tischen Näherungsausdruck der G1. 1.2 ( 13b) darstellen, und zwar für alle Werte von are (2) aus 
_ dem Intervall |© | <x/2. Vollzieht man gleichzeitig den Übergang von x zu x’ gemäß Gl. (2b), 
3 80 ea sich schließlich die Formel .(5). 
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ı® nl identifiziert, so entsteht die weitere as. Entw. (6) und in der er 
1.(7b) die as. Entw. (7) 
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In der zuletzt genannten Mberwiegt im allgemeinen der eine oder der Aare Se in der ge- 
schweiften Klammer. Sollen beide von gleicher Größenordnung sein, so muß schon gerade die 


Beziehung 
4 ‘ r 
2% + Sin 2a ; 
ER tg v’ — en ("= +ix,) 
0.4 - 2 


gelten. In der Gl. (6) tritt der Fall gleicher Betragshöhen für beide Summanden dann ein, wenn 
im besonderen #- = i- mit r > Oist. Die Gl. (6) nimmt in diesem Fall die Gestalt von Gl. (6a) 
an. 


n ni u(u2 —1) 
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® 
Aus ihr kann geschlossen werden, daß die Funktion M_ BE (W-4 r- Sin?o) bei festen Werten 
von & unendlich viele Nullstellen in bezug auf 7 hat. Das erste Glied in der Entwicklung von 
(6a) wurde zuerst von A. Erdelyi[4]angegeben. Bezeichnet 7,,„ eine der großen Nullstellen £ 
von (6a) in bezug auf r, so berechnet sich dafür aus (6a) der as. Näherungsausdruck (7). Ten 


ni en | 

an +3 1) To nn +5 1 5 
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3. Schlußbemerkungen. Auf den noch ausstehenden Fall eines 24» =1L,in = 
dem der Sp. von der zweiten Ordnung ist, soll an dieser Stelle nicht eingegangen werden, obwohl 
er sich auf dieselbe Weise erledigen läßt und in der formalen Durchrechnung sogar der bei weitem 
einfachste ist. Die praktische Verwertbarkeit der as. Reihe, die sich auf diese Weise dafür uf- 
stellen ließe, zeigt aber die gleichen Mängel ‚wie die entsprechende Reihe bei den Zylinderfunk- 
tionen H(; 2) (2) oder J, (2) im Falle eines v #z. Der Anwendungsbereich ist nämlich äußerst be- 
schränkt, denn die Abweichung des Verhältnisses 4 x/z von 1 muß stets die sehr scharfe For- 
derung erfüllen, daß für ein 2> oo der Quotient (1 — 4 x/z) : (2/6)? —0O.strebt. Zu einer we- 
sentlich besseren Näherung führt hier jedoch der nämliche Kunstgrift, den G.N. Watson 
bei den Zylinderfunktionen angewendet hat. Diese Behandlungsweise verdient außerdem noch 
aus dem anderen Grunde den Vorzug, weil sie sich auch in dem Falle anwenden läßt, daß außer 

x und zauch noch der hintere Parameter u sehr groß ist und für alle drei Veränderliche der Aus- 
druck ((1 + u)/2)? +1 —4x/z in der Nähe von O liegt. Wegen der dabei benutzten anderen 
een: soll dieser Fall den Gegenstand einer DE Arbeit bilden. EZ EI PER 
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flüssige Halbraum bei einer mechanischen Beeinflussung 
seiner Oberfläche. (Zweidimensionales Problem.) 
ee Von Fritz Sauter in Göttingen 


Als Vorstudie für die Behandlung des wesentlich verwickelteren Problems der mechanischen Beein- 
B h flussung eines elastischen Halbraums !) wird hier (in akustischer Näherung) die Wellenausbreitung in 
 -' einem flüssigen Halbraum untersucht, auf dessen Oberfläche ein Druck wirkt, der außer von der Zeit nur 
E a von der einen (cartesischen) Koordinate in seiner Oberfläche abhängt. Dann läßt sich das Geschwindigkeits: 
En potential in seinem räumlichen und zeitlichen Verlauf angeben als Doppelintegral über den mit einer be- 
B- 0 stimmten Einflußfunktion multiplizierten Oberflächendruck. 
Bi: S As preliminary work for the en of a problem essentially more complicated concerning the mechanical 
en. Ir influence on the elastic semi-infinite solid!), Ihe author analyses the propagation of waves (in acoustie approzi- 
ee mation) in a liquid semi-infinite body. T'he pressure that is acting on the Surface of the liquid depends (upon 
BF the time and) only upon one of the (Cartesian) coordinates in the surface. Then a veloeity potential can be 
$. r given as double integral of the pressure on the surface multiplied with a certain Green’s function. 
= ER > L’auteur traite la propagation d’ondes (approximation acoustique) dans le demi-espace fluide, dont la 
surface est exposee ü une pression, qui depend du temps seulement et d’une coordonnede swperficielle (carte- 
FE sienne). Ce travail est une &iude preliminaire pour le traitement du m&me problöme dans le cas consid£rable- 
Ri -- ment plus compliqu& du demi-espace solide elastique!). Alors le potentiel de vitesse dans son dependance du 
temps et de l’espace se pr&sente comme une iniegrale superficielle sur la pression mültiplide par une certaine 
fonction de Green. _ 
- B kayectse mpenBapureisHoß paboTk 1A MaByyenuAa SHAYNTEIBHO DoNee CIIOKHOH IIPO- 
u- : ÖJeMbI MEXAHHYECKOTO BO3BNSÜCTBUA HA ‚yupyroe MOAIYHPOCTpaHcTBo!) 3BNecb HCCHeNyerch 
e (B ‚‚aKYCTHYecKoM‘‘ IIPMONMSKEHUN) PACHPOCTpaHeHue BONH B IKUHAKOM HONYIPOCTPHETBE, 
= HA TOBEePXHOCTb KOTOPOTO AEeÄCTBYeT MaBIeHme, 3aBACAIee, KPOMe KAK OT BPeMeHHu, eIme 
TOIbBKO OT OAHOÄÜ u3 HPAMOYTOMBHEIX KOOPAHHAT HA HOBEPXHOCTH HONYIPOCTPAHCTBa. 
DIpocrpauctBensas uU BPeMeHHAA BABHUCHMOCTB TIOTEHLMAAIA CKOPocTeii MOskeT ÖBIT B TAHHOM 
cayyae IpencTaBlleHa B BUNe ABOÄHOTO HHTETPAIa OT HOBEPXHOCTHOTO AABIeHHA, YMHO3KEH- 
HOTO HA HEKOTOPyIO OUpenesIeHHyI EYHRINIO BIIHAHHA. 2 


1. Problemstellung. 


Der Gegenstand der folgenden Untersuchungen ist eine unendlich ausgedehnte ursprünglich 
ruhende Flüssigkeit mit anfänglich ebener Oberfläche, welche von einem bestimmten Zeitpunkt 
ab unter der Wirkung eines zusätzlichen Oberflächendruckes 


BEN ee Na. all) 


R- steht; dabei soll die @-y-Ebene mit der Oberfläche zusammenfallen und die z-Achse in die 
Flüssigkeit hineinweisen. Gefragt wird nach der Wirkung dieses Zusatzdruckes auf die Druck- 
und Bewegungsverhältnisse im Innern der Flüssigkeit. : - 
% Zur rechnerischen Behandlung dieses Problems sind zwei vereinfachende Annahmen 
. erforderlich. Und zwar wird erstens von der Schwerkraft, der Oberflächenspannung der Flüssig- 
keit und ihrer inneren Reibung abgesehen. Zweitens soll es sich um so kleine Störungen handeln, 
daß man mit den hydrodynamischen Grundgleichungen in der sogenannten akustischen Näherung 
_ rechnen darf. Dann läßt sich der jeweilige Zustand der Flüssigkeit durch ein Geschwindigkeits- 
potential @(x, 2,t) beschreiben, das der Wellengleichung 
ur Ders 129 
ee, Re 


« 


j ie n E = > E F 3 .. 
_ genügt, und aus dem sich die Strömungsgeschwindigkeit v und der Überdruck p nach den Be- 
_ ziehungen > 


ae) 


Es \ 0 
ws E 3 k v — grad od, pen a = ne Te er ne Sat Wert WR (1,5) 


> 
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berechnen läßt. Dabei bedeutet g die mittlere Dichte der Flüssigkeit. 

Die mathematische Aufgabe besteht nun darin, diejenige Lösung ® der Wellengleichung 
u "welche den durch das Problem vorgeschriebenen Randbedingungen genügt. Und zwar 
nuß erste: 
orgegebenen Funktion q übereinstimmen. . Und zweitens muß die Flüssigkeit in hinreichend 

‚Ber Entfernung vom Störbereich noch in Ruhe, d.h. im Normalzustand ® = 0 sein. 

i a hysikalisch wohl kaum realisierbare, mathematisch sicherlich bereits irgendwie 
em hier nochmals ausführlich behandelt werden soll, hat seinen Grund darin, daß 
‚ösungsmethode entwickelt wird, welche auch auf andere ähnliche Probleme an- 
"allem auf die physikalisch recht wichtige, rechnerisch wesentlich kompli- 
ngew. Math. Mech. Bd. 30 (1950), Heft 7. $ 
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zierter zu behandelnde Frage nach der Ausbreitung von Störungen in einem elastischen Halb- 
raum infolge irgendeiner Beanspruchung seiner Oberfläche. Über dieses letztere Problem wird 
in einer anschließenden Untersuchung berichtet; und die folgenden Ausführungen sind nur als 
eine Art Vorstudie hierzu gedacht. 


2. Die formale Lösung des Problems. 
Die allgemeinste, nur von x, z und t abhängige Lösung ® der Wellengleichung kann in 


der Form 
er + 
D(a,2z,i)= j do | dk f(k, w) e 
& Fan SE 


angeschrieben werden. Dabei ist wegen der zweiten Randbedingung das Vorzeichen der Wurzel 
im Exponenten so zu wählen, daß der Integrand eine in die Flüssigkeit hineinlaufende bzw. 
eine in dieser Richtung exponentiell abklingende Störung darstellt. Die erste Randbedingung 
verlangt wegen (1,1) und (1,5) 


oo 


rer sa) 000000 


+ + oo aan 
aa!) =vo jf do [ dkwf(k, w) Et 1 ee 
woraus man durch Fouriermäßige Auflösung nach f(k, ) und Einsetzen in (2,1) die Lösung 
des Problems in Form eines vierfachen Integrals erhält: 


-H oo + + oo \ 
1 do i (kaan ar V&-r 2—a(t- 0) - 
D(z, 2,1) = —— dk dx’ dt’ g(&,t')e (259) 
4n?io {00} x 


Doch lassen sich, wie im folgenden gezeigt wird, die Integrale über ® und % allgemein auswerten, 
so daß man schließlich die Lösung unter Verwendung einer „Einflußfunktion‘‘ quellenmäßig in 
der Form 


. 


> + 
da,2.,1)= [.da [dt gast) Fle—-0,31—t) 2» 22.22.24 


—no 


darstellen kann. Die Aufgabe reduziert sich somit auf die Berechnung der Funktion 


4 + = 
' 1 la ilk&+ VS, -R z—ot 
Fo) = a, N" ) De Ber 


Zunächst eine Bemerkung wegen der in diesem Integral auftretenden Singularität bei 
©—=0: Diese ist dadurch bedingt, daß das Geschwindigkeitspotential aus dem Druck durch 
eine Integration nach der Zeit hervorgeht und daher aus dem Stördruck g grundsätzlich nur bis 
auf eine zeitlich konstante Lösung der Wellengleichung bestimmbar ist. Nun unterscheiden sich 
der Cauchysche Hauptwert des o-Integrals von den Integralwerten, die man durch Umgehen 
des Punktes = O0 in der komplexen &-Ebene erhält, nur um Ausdrücke, welche nicht von der 
Zeit abhängen. Da die Flüssigkeit ursprünglich ruhen soll (d.h. ®=0), hat man unter den 
verschiedenen Möglichkeiten denjenigen Weg für die o-Integration zu wählen, der kein zeitlich 
konstantes Glied liefert. Doch ist es nicht notwendig, sich schon jetzt auf einen solchen Inte- 
gralionsweg festzulegen, vielmehr genügt es, späterhin solche zeitlich konstanten Zusatzglieder 
von ® nachträglich zu streichen. Überhaupt brauchen bei der folgenden Auswertung von (2,5) 
nicht alle Schritte im einzelnen mit der ganzen mathematischen Strenge durchgeführt zu werden, 
da schließlich die so gewonnene Schlußformel doch noch besonders auf ihre Richtigkeit unter- 
sucht werden soll. 

Bei der Auswertung von (2,5) ist es zweckmäßig, die durch die Abstrahlbedingung wegen 
des Vorzeichens von YIw/c)? — /* erforderliche Fallunterscheidung &$ 0 durch eine Umformung 
“ nach dem Schema ‘ . 


+» e: oo , 
[ do G(o) Pu: = 2#e| [doGlo) ie ern ie (2,0 
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zu umgehen. Diese Umformung ist möglich, da wegen der Realität von ® und damit auch von F 
naturgemäß @(— 0) = @*(o) sein muß. Dabei bedeutet das Sternchen den Übergang zum 
Komplexkonjugierten. Man kann sich daher im folgenden durchweg auf positive w-Werte 
beschränken. N 
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x N e 
| Ferner ist es angezeigt, den Wurzelausdruck im Exponenten von (2,5) durch die Sub- 
stitution Er 
we w* © 
k=-sina, EEE 00 ne een) 
c a 6 


zu beseitigen. Dabei bedeutet die neue Veränderliche a, sofern sie reell ist, offenbar den Neigungs- 
winkel der ebenen Wellen, aus denen sich das Integral (2,5) zusammensetzt, gegen die Oberfläche. 
Setzt man außerdem noch zur Vereinfachung 


EN RE A 0 15.0.2 0 ee RE 5:3} 
so erhält man für F statt (2,5) den Ausdruck 
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Wegen der Definition von & und der Abstrahlbedingung hat man, um zu einer für alle 9-Werte 

zwischen —r/2 und + x/2 gültigen Formel zu kommen, ursprünglich geradlinig von —rz/2 + 

' nach —r/2, von da nach +/2 und schließlich nach +r/2 —i ® zu führen (vgl. den gestrichelten 

Integrationsweg in Bild 1). Er kann aber, da der Integrand nur im Unendlichen singulär ist, 

auch irgendwie anders, also z. B. in.der Art der ausgezogenen Kurve von Bild 1 geführt werden, 

ER sofern er nur bei —r/2 +t oo aus dem Unendlichen kommt und bei -+7/2 — ı oo wieder ins 
Unendliche führt. 
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Dieser Verlauf des Integrationsweges könnte nun in Hinblick auf die Analogie mit der 
Sommerfeldschen Integraldarstellung der Hankel-Funktionen nahe legen, erst die «-Inte- 
gration auszuführen. Man kommt aber hier und im besonderen bei komplizierteren Problemen 

_ viel schneller zum Ziel, wenn man in (2,9) erst die &-Integration ausführt. Man gelangt so un- 

mittelbar zu der Formel | 

- 1 & CHR 
ä 1 dacos«a 

j R cos («—d) — ct E10) 
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sofern man den Intefrationsweg in der &-Ebene so führt, daß dauernd Im{R cos («— #9) — ct} > 0 

st. Da letzteres nur im zweiten und vierten Quadranten eines durch den Punkt & = ö gelegten, 

zu den bisherigen Achsen parallelen Koordinatensystems gilt, hat man den Integrationsweg 
“ jetzt in der in Bild 2 durch die ausgezogene Kurve dargestellten Weise zu führen'). 

ch verdanke diese einfache Begründung für den Übergang von (2,9) nach (2,10) Herrn Dr. G. Leibfried. 


atte ursprünglich das w-Integral nur bis zu einem großen Wert Q erstreckt und dann das die Größe %2 
iltende Glied hinsichtlich der a-Integration mit der Sattelpunktmethode ausgewertet (Sattelpunkt bei 
Ö , wobei i dieses Glied im. Limes Q > & wie 1/2 verschwindet. 
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b: % 
Natürlich kann man nun den Integrationsweg wieder weiter deforiiieren, hat ab Setzt. 
auf die neu aulgetretenen Pole bei «= ® +arecos (ct/R) zu achten. Diese liegen für ct <R $ 
symmetrisch zum Punkt & — d auf der reellen Achse (?, und P,in Bild 2), für et > R auf der Ex 
Vertikalen durch diesen Punkt (9, und Q, in Bild 2). esichtlich muß der Integrationsweg ZzUu- Kar 
nächst so geführt werden, daß P, bzw. Q, rechts, P, bzw. Q, links von ihm liegen. Re ER 
Das Auftreten von Polen im Integranden läßt "vermuten, daß sich das Integral (2, 10) a, 
irgendwie auf den Umlauf um einen dieser Pole zusammenziehen und so nach der Residuum- 
methode auswerten läßt. Doch stört hierbei zunächst die Tatsache. daß der Integrationsweg 
aus dem positiv-imaginär-Unendlichen kommt und ins negativ-imaginär-Unendliche läuft. Da 
aber ein Integral mit reellem Integranden nur ins Komplex-konjugierte übergeht, wenn man 
den Integrationsweg an der reellen Achse spiegelt, und da in (2,10) der Realteil zu nehmen ist, 
so daß eine solche Spiegelung des Integrationsweges oder eines Teils desselben den Wert von F 
nicht ändert, kann man den in Bild2 ausgezogenen, bisherigen Integrationsweg durch den 
sestrichelten ersetzen, welcher sich von dem bisherigen nur durch das Hinaufklappen seines 
unterhalb der reellen gelegenen Stückes unterscheidet. 
Jetzt kann man den neuen Integrationsweg nach oben ins Unendliche ziehen und erhält von 
dort wegen lim cosa/(R cos («—9)—ct) = e-'?/R einen zeitlich konstanten Beitrag zum In- 


a>P-+tie 
tegral, der einem stationären Strömungsfeld entspricht und daher gemäß den früheren Bemer- 


kungen weggelassen werden muß. Einen endlichen Beitrag erhält man daher nur, sofern einer der 
neuen Pole oberhalb des gestrichelten Integrationsweges der Abb. 2 liegt. Dies ist nur der Fall, 
wenn ct> Rist, in welehem Fall man aus (2,10) durch Residuenbildung an der Stelle @, zum Wert 
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kommt. Für et < R verschwindet F (x, 2,t)identisch. Setzt man dies Resultat in die Formel (2,4) 
ein, so erhält man wegen (2,8) als Lösung des Problems ze: 
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3. ÜUberorüfung der Lösung und Dirkossiar- 


Die eben gefundene Lösung des Problems läßt klar erkennen, daß sich die Wirkung einer . 
Störung von der Oberfläche in das Wasser hinein mit Schallgeschwindigkeit ausbreitet. Im 
besonderen ist eine zur Zeit 0 punktartig im Nullpunkt einsetzende Störung zur Zeiti>0. 
gerade bis zur Zylinderfläche «® + 2? = et? gekommen ( Stoßwellenkopf). Mit dieser Feststellung 
ist die Randbedingung im Unendlichen als erfüllt erwiesen. Daß die Lösung auch die. Rand- 
bedingung an der Oberfläche befriedigt und außerdem der Wellengleichung ‚genügt, also tat- 
sächlich die richtige Lösung darstellt, kann man etwa folgendermaßen einsehen: 

Setzt man die obige Lösung in die Wellengleichung ein, so wirken sich die in ihr enthaltenen : 
Singularitäten beim zweimaligen Differenzieren recht störend aus. Setzt. nsaber PARERUSE. = a 


x = —zsinhu, !=t—- Corte ie Kae 
6 


so wird der Integrand im Endlichen durchwegs regulär: a Se SR 
t 4 ++ >) r N = 
DR ol ER 
Er | a) dv coshv AE: Pau Ba 


—o 


Setzt man diesen Ausdruck in die Wellengleichung ein, so kommt maı 
zu einem Doppelintegral, dessen Integrand sich als vollständiges -Difl 
das Integral verschwindet daher, sofern q(x,t) für u + , > 
t— — oo hinreichend stark gegen Null geht, was hiermit nochmals. a 
werden möge, Bil: 
Auch zur Überprüfung der Randbedingung an der Öberfiäc 
in'der Lösung transformieren, Setzt man hier 


iur — URN, ee A 


RR EN 
üs sige Halbraum bei einer mechanischen ei Hlussimg 153 
so geht die Lösung über in 55 y 

en: | Dia 2,0 = at iR - q(® i 3,4 
er De el 3 2% er oc - E mu it? En +7) q RA) — s/c) ur ( s P% 
nz “ = zVırr 

t Hier läßt sich der Grenzübergang z—>0 ohne weiteres uhren, und es gilt die einfache Be- \ 
Br 2 ziehung 2 
rn, am um... 00, 
te Ti 
Be - Daraus findet man für den Zusatzdruck an der Oberfläche 

Fr, Br 02.1 2, — es 

a > a. a 24, t— s/c) —— [a IN _ 9 (3,6), 
# also talsächich den geforderten Wert. 

T Aus der Formel (3,5) kann man übrigens sofort die &-Komponente der Geschwindigkeit 

‚eines Punktes an der Oberfläche ermitteln: 

K 2 .96(x,0,t) 1 5 0g9(2,t— s/c) : 
E: Sy | ; 77 (2,0, t) = se Free = en en a ae series (3,7). 
er = Viel diffiziler ist die Berechnung der 2-Komponente, RE die man nach einer längeren, hier nicht 
 wiederzugebenden Rechnung den Ausdruck 
BEN ee 

2; SEER q(z, t) 1 | j cos?p 9 

er t — — > ud ’ t — . * . ’ 

| v,(%, 0, t) SE "oe, ds Ä Br: Se 5,1% ssinp s/e) (3,8) 

D rg Pu : 
findet. 
Um schließlich noch einen Überblick über die Verhältnisse im Innern der Flüssigkeit zu z 


erhalten, dürfte es genügen, die Lösung für den Fall zu untersuchen, daß die Druckeinwirkung 
auf die Oberfläche nur in einem ganz schmalen Bereich erfolgt und in zeitlicher Hinsicht als 
° Rechteckstoß verläuft: 
DIL. ES 0 für t<0O 
- g(2, 1) = 


33 - 195(2) für 1>0 


_ 


93,07 


bei zeitlich konstantem Q. Dann erhält man für das Geschwindigkeitspotentia] aus (2,12) im 
let, > Tr +22? den Ausdruck 
Q 2 


Ö(z, Er eure ae dee, Ze Did {3; 10), 


& = während es für ne < J&? +2? verschwindet. Hieraus findet man für die Be 
ee im Bereich hinter. der Stoßwellenfront 


es ae 1. | 
Si . ee ee me 7, 22 292 — 22 
DAR Vx ts +2 Jer 22 E an | % 
N Z Ba 2 ; : 2111), 
x RS . i ® BE ES BET EEE TER EUR ) 
i V, = meol®+7 en Fa rl° % zZ 
sich für. den Druck die Beziehung | | 
RA u 2 | Ge er Br. Er 3,12 
X a ER a +2 ara Beuel e j5 6... ° (3,12) 


ittelbar. hinter der Stoßwellenfront werden Druck und Geschwindigkeit unendlich 
dann gegen das Störzentrum hin ab. Dabei überlagert sich über einen vorwiegend 
Ste ‚Bwellenfront wirksamen, radial gerichteten Strömungsanteil ein in der N ah- 
ngsstelle wichtig werdender zweiter Strömungsanteil, der von einer dipol- 
natio! on einer Quelle und einer Senke herzurühren scheint und die durch den Zu- 
Ausweichströmung re ‘An der Grenzfläche selbst verschwindet in 
‚ die Flüssigk Bet sich hier Be Aentzecht zur Oberfläche. 
r ei as“ } ‚ > 8 a Bags 
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Morphologisch homogene Funktionen und ihre Erzeugung 
durch statistische Superposition von Elementfunktionen. I _ 
I. Teil: Allgemeine Theorie 


Von Hanswalter Gresekus in Göttingen 


Die Menge der morphologisch homogenen Funktionen (mhF) und die sie charakterisierenden Trans- 
formationen ‚‚Spektrum‘“ und „Produktdurehschnitt‘‘ werden allgemein definiert. Dann wird aus dieser 
Menge diejenige Teilmenge behandelt, die man durch statistische Superposition quadratisch integrierbarer 
Elementfunktionen erzeugen kann, und zwar in der Weise, daß zuerst eine ganz spezielle Klasse sehr ein- 
facher Funktionen betrachtet und diese dann sukzessive verallgemeinert wird: Spektrum und Produktdurch- 
schnitt werden als Funktionen der entsprechenden Transformationen der Elementfunktionen dargestellt. 


The class of morphologically homogenous functions (mhF) and the transformations „spectrum‘‘ and 
„produet-average‘‘ characterising the mhF are introduced. The special class of mhF which can be generated 
by statistical superpösition of element-functions, quadratically integrable, is discussed, first by considering 
very simple functions, and Ihen by successive generalisation: Spectrum and product-average are represented 
as functions of Ihe corresponding transformations of the element-functions. 


On donne la definition de la classe des fonctions de morphologie homogene (mhF) et des transformations 
„‚specire‘‘ et „produit moyen‘‘, caracterisant les mhF. Alors la classe speciale des mhF', qui sont construites 
par superposition statistique de functions d’element, est traitee de maniere qu’on commence ü consid.erer 
des fonctions, Ires simples, qui ensuite sont g£neralisdes successivement. Specire ei produit moyen se 
presentent comme fonctions des iransformations correspondantes des fonctions d’&l&ment. 

Jaerca oömee ompenenenne MHOKEeCTBA. MOPPoNOoruyeckH OAHOPOAHEX byHkumü (mhF) 
u npeoÖpasoBaHunä ‚„CIEeKTPp‘“ m ‚„‚cpenHee ITPousBeeHue‘, XApakTepusylollux TU PyHEOun. 
BarTeMm U3 DTOTO MHOMKECTBA ÖepeTcHd TAKOE CHEMHAJIBHOE MHOIKECTBO, KOTOPO® MOKET ÖBITb 
IO1yyeHO CTATUCTUGECKUM HAJOKCHUCM KBAAPATHYECKH UHTErPHPYEeMBIX YHRIIMÄ-HIIEMEHTOB, 
IpuyeM CHAyala PACCMATPUBAWTCH COBCeM HPOCTHEe PyHKINM, KOTOPEI® IHOTOM IOCIHIENOBA- 
TEIBHO 0000MAMTCH! „‚CHERTP‘‘ MH „cpenHee IPonsBeNeHne‘‘ IPeNCTABIAMTCH Kak PYHKIUU 
COOTBETCTBEHHBIX IPeOoOpasoBaHnft PYyHEIMH-9IeMEHTOB. ' 


I. Einleitung > 


In der Physik, insbesondere in der Astrophysik, Geophysik und Meteorologie begegnen 

uns häufig Funktionen — ihre Abszisse ist fast immer die Zeit — die durch Überlagerung einer 
Menge von Elementfunktionen mit mehr oder weniger variabler Gestalt gebildet werden, derart, 
daß diese Elementfunktionen auf der Abszisse irgendwie statistisch verteilt sind. Liegen solche 
Funktionen als Zeitreihen vor, so nennt man diese wohl ‚‚statistische Reihen mit Erhaltungs- 
tendenz‘‘ (unter Umständen noch dazu ‚‚mit Wiederholungstendenz‘‘). Als einige der bekannte- 
sten Beispiele seien genannt — aus der Elektrodynamik: die Stromstärke in einer Elektronen- 
röhre!); aus der Optik: jeder Lichtstrahl (atomistisch betrachtet)2); aus der Astrophysik: die 
integralen Maße der Flecken-, Fackel- und Koronaintensität, sowie die Intensitäten der Ultra- 
violettstrahlung, der Korpuskularstrahlung und der erst vor kurzem entdeckten elektromagne- 
tischen Ultrakurzwellenstrahlung der Sonne; aus der Geophysik: die Schwankung des erd- 
magnetischen Feldes und die Nordlichthäufigkeit; aus der Meteorologie: die zeitlichen Varia- 
tionen des Luftdrucks, der Temperatur, usw. Solche Funktionen — wir nennen sie mathematisch 
„morphologisch homogene Funktionen‘ (mhF) — sollen hier ganz allgemein behandelt werden, 

"und ünsere Aufgabe wird darin bestehen, eine Methode zu finden, welche die Einflüsse der 
Statistik eliminiert und uns dadurch Aussagen über die mittlere Gestalt der sie erzeugenden | 
Elementfunktionen zu machen gestattet. | 

Die folgende Theorie ist zwar ganz allgemein mathematisch durchgeführt, jedoch steht 
hinter den einzelnen Schritten immer die Anschauung der wirklich vorkommenden Funktionen 
und des Mechanismus ihrer Entstehung. Insbesondere ist dabei an die solaren Funktionen 
gedacht, wie die Theorie in erster Linie auf eine Anwendung in der Physik der solar-terrestrischen 


Beziehungen zugeschnitten ist, ohne jedoch hierauf irgendwie beschränkt zu sein. e E 
II. Die Definition der morphologisch homogenen Funktion. Spektrum und Produkt- 
durehschnitt ee 


$1. Die mhF und ihr Spektrum = 


is werde eine in jedem endlichen Intervall definierte und darin quadratisch integrier- 
bare Funktion F(t) betrachtet. Wir ersetzen sie durch MM i 


3 _[F(t) im Intervall (7,, 7,) = 74% 
F 12(d u 0 außerhalb 1 2. . . $ On) . . ” . * . »2 (1), 
1) Das bedingt den sog. Schroteffekt. ’ ne 


°) Das ist ja die Ursache der fundamentalen Eigenschaft des Lichts, daß nur solche Strahlen, die von 
ein und derselben Lichtquelle ausgehen, interferieren, und auch dann nur, wenn die Differenz der Lichtwege 
nicht wesentlich größer als die Kohärenzlänge ist. > Ä . 


‘ 


Z. angew. Math. Mech. . . s £ 
Ei sy Nr. 5/6 Mät/Juni 1950 Giesekus, Morphologisch homogene Funktionen und ihre Erzeugung 155 
a a E05. ee... 


——— 


und können diese Funktion dann durch ein Fourier- Integral darstellen: 


r_ iot 
Fu) [: Bd 
12 ) N 12 V>r7 (2) 
mit 
co 7, 
7 ( [r n e-iet $ [ru Besen 5 
0, = N die 
12 1) 5 Vz (3) 
Dabei gilt die Vollständigkeitsrelation: 
. T, Re nd 
Kal er JS Fo)" do Er rer (AR 


Ist F(t) {) nicht nur in jedem endlichen Intervall, sondern auch in der Grenze 7,—— x, T,— © 
quadratisch integrierbar, so ist der Grenzübergang ohne weiteres auszuführen und # or selbst 
strebt gegen 


m —iot J 
on dog) 2 en 


wie die Theorie des F ourier-Integrals zeigt. 
Wir fragen nun, ob ein Grenzübergang auch dann möglich sein kann, wenn fr t) dt 


nicht existiert, beschränken uns dabei aber auf den Fall, Hal PD 
F? = lim [ro SR EEE 
‚im TH. dt el) 
T,— © 


einen endlichen Wert hat. Dann muß nach Gl. (4) auch 


A RE im | Ro) do - ae NER N 


= Si F 
e lim j 
4 VT,— 


existieren, woraus aber nicht folgt, daß dies auch für 


K(o) = lim K,,(® 0) = tim I 7 ER RAR 2 
selbst der Fall ist, und erst recht nicht für 
i F (©) 
lim en 
Y7T,—T, 


Wie wir später noch zeigen werden, existiert der letzte Ausdruck in der Grenze nur in ganz 
besonderen Fällen, nämlich dann, wenn F(t) periodisch oder fastperiodisch ist,.d.h. aus einer 
Summe periodischer Funktionen mit nichtrationalem Periodenverhältnis besteht®). Solche 


Funktionen, für welche K (w) existiert und eine integrierbare Funktion darstellt, müssen offenbar 


die notwendige und hinreichende Bedingung erfüllen, daß sie im Sinne der Statistik homogen 
sind, d.h. daß .bei Unterteilung der. Abszisse in genügend große Teilintervalle (7,1, 7,) die 


TR »(®©) um eine mittlere Funktion herum schwanken, wobei die Größe dieser Schwankung 
beliebig klein gemacht werden kann, wenn man nur die Intervalle genügend groß wählt®). Wir 


nennen,diese Funktionen deshalb allgemein „morphologisch homogene Funktionen‘ und K (®) 
ihr er a 


3) Wir Herden im folgenden die fastperiodischen Funktionen oft zu den periodischen Funktionen hin- 
zurechnen, 


4) Als einfaches Beispiel, für welches lim wi K(0) do existiert, nicht aber K(w) — lim K,.(®), sei die 


ER Funktion F(t) = sin (1?) genannt. Dagegen ezistiert Kt @) wohl für solche Funktionen, die wie sin t:cos 1/t oder 
cos t+ sin 1/t in ihrer Homogenität nur derart gestört sind, daß der Einfluß dieser singulären Störung in der 


Grenze verschwindet. Wir werden von solcher Art von Störungen i in der Folge absehen, da diese auf unsere 
„berlegunde en keinen Einfluß haben. 
ir sagen ‚Spektrum‘ (in der Sprache ‚der Optik: Intensitätssprektrum), weil wir den Ausdruck 


N Feriodogamm' (= Amplitudensprektrum) für F(o) und entsprechende Größ&n verwenden wollen. 


_ 
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Unsere Definition muß aber noch in einer Hinsicht präzisiert werden. Sie ist nämlich nur 

dann vernünftig, wenn wir den Begriff der integrierbaren Funktionen so weit fassen, daß er die 1 
uneigentliche Diracsche ö-Funktion Br 


5) = 55 | er du u N 


und Summen dieser Funktionen mit einschließt. K (©) wird nämlich stets solche ö-Funktionen ; 
enthalten, wenn in F(t) periodische Bestandteile vorkommen, im einfachsten Fall schon, wenn 
F-=0ist. Es wird im folgenden jedoch « = zweckmäßiger sein, diese 6-Funktionsanteile aus 


K(o ) herauszunehmen und in einer Folge RS (or) zusammenzufassen, Leise daß wir an Stelle 
von Gl.(4) ausführlicher schreiben können: 


EN Veen RE Sn 


a a) 


Wir haben also K (») in ein kontinuierliches Spektrum K® (») und ein diskretes Spektrum 


K®) (@„) zerlegt. K® (w) ist dann nicht nur integrierbar, sondern sogar eine quadratisch 
integrierbare Funktion. Für eine gewöhnliche periodische Funktion der Periodenlänge © 


© 2zinZ A 
= 2 Bez). . N en Eee ll) 
mit ; 
; 0/2 : 
3 97 —2r2inz S 3 k 
B(&n)-5 | Pe ne een 
— 9/2 n 
sind die 


2. 


h 


712), 


Keen 


also das Absolutquadrat der normalen Fourier-Koeffizienten. Im Falle der a | 


Funktion 7 2 | 
Pi „2. Ploy)eiont, wo DA er ee: (102) a 
und 5 | | = 
2 1 T/2 EN Re Re 
Y (o,) = lim F | PER Se ar 

x FE Ara ERBE" 

gilt entsprechend: 
N -- < =‘ £ 

Ko) = | Yo, Re. ee lea) 

Ri: . Das Spektrum der gewöhnlichen periodischen Funktionen ist in diesem letzten Fall natür- 


Jich als Sonderfall eingeschlossen, nämlich als Spektrum mit äquidistanten Linien. ah 
Derartige periodische Bestandteile werden wir also in den mhF im allgemeinen immer vor- 
finden. Natürlich gehört auch die Klasse der rein a Funktionen, d.h. der Funktionen, 


die nur ein Linienspektrum besitzen, für die also R° '(o) = Oist, hinzu. Doch wollen wir diese = 
besonders einfache Unterklasse im folgenden nicht nz betrachten. Peer = a ER Rare 3 


$2>Der Produktdurchsenaeie IR ee 


Pi. Wenn K” (©) als kontinuierliches Spektrum quadratisch Integrie 
ein Linienspektrum ist, so wird man ja durch Umkehr der Fourier-Tra 
Funktion K(r) finden können, wel ns ebenso einen quadratisch integrierk D 
und einen periodischen Anteil K(P (r) besitzt. Wir können sieleicht. 


an F (wo) 12 es 
Bu | re If 
12 ) Te T; FIREETR De Tr ei Srworue N 
Substitution von !’ durch = t— f' Be B-: EB 


oo 


., 2 EP. 
K.(®) u er Fe —T, Iruor /12 vo: 
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Wir können nun diese Transformation kehren: 5 
a Be = [Rn — 
Bi: ge En Vor 
Br da der Grenzwert unter dem Integral (im oben erläuterten Sinne) existiert, hierin zur Grenze 
Br - übergehen: € : 


ee Ka) = - [ko on DE [er (0 B IK ;..: (5). 

Nach Gl. (14) ist Bank. | 

Be SEEN Ki) “2 1 lim { [r (it dt 16 

= RR? | Pr nt) -12 ) ] 12 ( —7) . . . . . - . . ( ) 
oder auch 

Er - 7,+1/2 

A 1 BUT 1 


5 $ K(r) = — Ft) Fit — — lg Fi)Fie—v)di®) .. . (16a). 
N. y2r urees Im Fr/2 
Wir nennen diese eiforrähion den ‚„‚Produktdurchschnitt‘ von F(t). Sein quadratisch 


3 integrierbarer Anteil KU (7) ist durch den ersten Summanden, sein periodischer Anteil XP (r) 
E: durch den zweiten Summanden von Gl. (15) Begeben. Einige wichtige Eigenschaften des Pro- 
8 duktionsdurchschnitts seien kurz zusammengestellt: 

E:» 1. K(r) =K(—r), d.h. K(r) ist (wie in der Anmerkung zu -Gl. (16a) schon bemerkt 


wurde) eine symmetrische Funktion. 

2. K(r) hat für x = oein absolutes Maximum, welches an keiner anderen Stelle wieder an- 
genommen wird. (Eine Ausnahme bilden natürlich die rein periodischen, nicht aber die fast- 
periodischen Funktionen.) Dies bestätigt man unmittelbar, wenn man z = O0 in Gl. (15) einsetzt, 
da ja K (w) seinem SED, nach stets positiv ist. Anschaulich folgt das auch sofort aus der Tat- 

sache, daß 


DE K (0) = ee re N 4 (17) 
gilt. | 


Yan BE, T2 # 
aus 1 se 
FRE 3. K(d)=lm — [ko a - BE 22 {1B), 
= e | T>x Pin Y2 


wie man durch rzerätlon von Gl. (14) und Grenzübergang sofort erhält. 


Wie die vorangehenden Betrachtungen zeigen, sind Produktdurchschnitt und Spektrum 
durch Fourier-Transformation umkehrbar eindeutig zugeordnete Größen, welche die mhF 
gleicher Weise repräsentieren. Es muß jedoch darauf hingewiesen werden, daß diese Funktion 
elbst in ihrer ganzen Erstreckung nicht durchrirgendeine Summen- oder Integralformel nach 
2 ‚der Fo urier-Enwicklung dargestellt werden kann, da ja, wie schon bemerkt, der Ausdruck 


- = in der Grenze nicht existiert. Das hat, wie später noch deutlicher werden wird, seine Ur- 


N in, daß diese Funktioner i im allgemeinen einen statistischen Charakter tragen. Un- 


zw. K (®) können dann ihrem Wesen nach als Grenzwerte von Erwartungswerten auf- 
ßt werden. Für die Phasenfunktion gibt es aber prinzipiell keinen Erwartungswert. Alle 
hen Werte sind Bis wahrscheinlich, was as mit der Gleichwertigkeit aller Abszissenpunkte 
en hängt A - 


a 
)ieser . te ist nur im Limes mit dem vorhergehenden identisch, nicht aber für endliches 7, 


Die: 

man t/2 in den Integrationsgrenzen ist nur dann von Bedeutung, wenn der Grenzübergang 
ich vollzogen wird. Er bewirkt in diesem . Falle, daß X(r) auch für endliches Integrationsintervall 

isch bleibt. 

gilt natürlich nicht für die periodischen Bestandteile von F(t), da diese ja keinen 

sus Thnen: kommt, stets eine wohldefinierte Phase zu, 


a ra. u #794 
N et u Mn, ‚ 
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$3. Spektrum und Produktintegral quadratisch integrierbarer 
Funktionen 

Zum Schluß dieses Kapitels seien noch einige Transformationsgrößen quadratisch integrier- ; 

barer Funktionen eingeführt, welche in den folgenden Überlegungen eine wichtige Rolle spielen k 

werden. 


Wenn f(t) eine im Intervall (— 0, 00) definierte und dort quadratisch integrierbare Funk- 

tion 8) bedeutet, so sei } 

2 A e eiet 

fo) = [ro eben... Se. Snack) 

In | 

ihr Periodogramm und 4 3 

kloyset Wlan... 00 Pe 

ihr Spektrum. Analog zum Produktdurchschnitt existiert hier das „Produktintegral“ 

© oo I 

e eior 1 A . 

k (T) 1 (©) nn do A 7 (t) fit we T) dt OT ee,» (21). | 

—&© y2r y2r —o& > . P 

Diese Größe hat die dem Produktdurchschnitt ganz entsprechenden Eigenschaften: | 

FR k(r) =k(— Tr), 

ER | | 

2. k (0) = Te = fr (t) di= Max{k ()} a (22), 
an, | 


$ Sunao- 41 Io - us EEE ER IE)? 


Die Beweise gehen ganz genau so wie bei den analogen Sätzen für den Produktionsdurchschnitt 
der mhF. | 


ce ec ee 


III. Morphologisch homogene Funktionen als statistische Überlagerung quadratisch . 
integrierbarer Elementfunktionen | 
$1. Die Zerlegung der mhF in quadratisch integrierbare 

Elementfunktionen B 
Aus der Forderung, daß die mhF in jedem endlichen Intervall quadratisch integrierbar | 
sein muß, folgt sofort, daß man sie in quadratisch integrierbare Funktionen zerlegen kann. Wenn | 

wir die rein periodische Funktion ausnehmen, gibt es indessen unter allen möglichen Zerlegungen 
keine einzige, welche die Funktion in regelmäßiger Anordnung in gleichgestaltige®) Elementfunk- 
tionen aufteilt. Nun ist zwar eine solche Zerlegung stets auf unendlich viele Arten möglich, den- 
noch gibt uns dieser Gedanke einen Gesichtspunkt, wie wir eine große Unterklasse der mhFin ° | 

ihrer Struktur beschreiben können, indem wir nämlich für diese die Existenz einer besonders ein- 
fachen „natürlichen“ Zerlegung annehmen. Wir wollen nun allerdings nicht von einer solchen f 
mhF ausgehen und diese Zerlegung aufsuchen (das wird im allgemeinen nicht. durchführbar sein), 
sondern wir wollen den umgekehrten Weg einschlagen: Wir geben die Elementfunktionen vor 
und erzeugen daraus durch statistische Superposition eine mhF. Wie schon bemerkt, erfassen wir 
mit dieser Darstellung nicht alle möglichen mhF, z. B. nicht solche, die sich durch stückweises 
Aneinandersetzen quadratisch integrierbarer Funktionen in natürlicher Weise aufbauen lassen 10), 
Allerdings dürften diese Funktionen in der Anwendung von wesentlich geringerer Bedeutung sein 
als die von uns untersuchte Klasse. r 2 


$2.Aus lauter gleichgestaltigen Elementfwnktionen aufgebaute mhF 
Wir beginnen mit der Untersuchung der wohl einfachsten überhaupt möglichen mhF 


Fo) = I ft—i) RE Re SENSE Ce (24), 


x 
=. 


worin f{f) eine quadratisch integrierbare Funktion bedeutet und die 2; mit konstanter Dichte sta- 


®) Wir bezeichnen quadratisch integrierbare Funktionen im folgenden stets mit kleinen, mhF dagege n 


mit großen Buchstaben, ebenso die zugehörigen Transformationsgrößen, auch wenn diese einen andern 
rakter ‚ragen als die von ihnen repräsentierten Funktionen. 

m d. h.: gleich bis auf eine Verschiebung in der Abszissenachse. 0 ie Fule 

....0) Eine einfache Funktion dieser Art wird in $4 kurz behandelt, die periodische Funktion mit sta- 

tistisch verteilten Phasensprüngen. “ Ar 
\ E% 

De 
4 eu 


a 4 
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tistisch auf der Abszisse verteilt sind. Für die Rechnung werde diese Funktion durch die im In- 
tervall (— 7/2, 7/23 damit wesentlich gleiche Funktion 


N 
Fi) = ) AP DE 1 E; (24a) 
= 
angenähert, wobei ; 
Bet im Intervall (— 7/2, 7/2) rn 
s ww ZA Ro a 74 SS e 
die Wahrscheinlichkeit für die 2; bedeutet. 
Wir bilden nun den Produktdurchschnitt): 
% n oo N ? 
K()=7 | > ft) ft — 9) di = => Beer) 6.2 (26), 
Ber el ER k=1 


wo k(r) das Produktintegral von Ft bedeutet. Damit nimmt der Erwartungswert von K’(r) die 
Form an: 


B{K'(d} — f-ww--- Wii) K'ddty:--diy 


—o 


Ta Tj2 : 
N 1.5 1 WER 07); 
ce —Tja —Tj2 - 


Wir wollen das Integralglied umformen und führen deshalb eine neue Variable 9 — r— Ben 
ein: 
ra 7/2 ma *-y+T/2 
1 
kr—bt+h)dydy = 


ra —Tj2 —Tj2 7--T]2 


k(d) dddı,, 


TR 


So wie dieser Ausdruck dasteht, können wir ihn nicht weiter vereinfachen, wohl aber können wir 
seinen Wert im Grenzfall T— » sofort angeben. Denn dann wird 


2-1, + 7/2 © © a 
k(d) Ad [ra — [oa =: 82, 
7-17] . ee = 


also von 4; unabhängig, und somit nimmt das ganze Integral diesen Wert an. Wir erhalten 
asymptoticsh, denn die "Doppelsumme gibt den Faktor N Fe 1) 


65 E{K (N “ k(r) (9) dd 

Y Führen wir hierin die „Konzentration“ (d.h. die mittlere Anzahl der in die Intervalleinheit fal- 

2 lenden Elementfunktionen) 

; N 

2 ee er a 28 
I | eG (28) 

# ein und lassen dann N und T beide so gegen unendlich streben, daß dieser Quotient erhalten 

e. bleibt, dann geht der Erwartungswert von K’(r) in K(r) selbst über und wir erhalten dafür: 


Ko)=c-kn)+e: [io dd =: k( ae [roaı) Be (29). 


Der: zweite Summand stellt aber gerade das Quadrat des Mittelwerts von Ft) dar, so daß wir das 
Ergebnis auch schreiben können: 


Bene... ur (29a). 


WARE 


ar IR Im folgenden wird der Normiorungsfaktor Na Produktdurchschnitt und Produktintegral stets 
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Wir erkennen hier sofort die beiden Anteile: 


Hieraus können ce und k(r) dann auch einzeln bestimmt werden: 


Für K(o) können wir auf Grund der Transformationsgleichungen die entsprechenden Formeln x 52 
sofort angeben. Wenn k(w) das Spektrum von f(t) bedeutet, so ist “= ae 


und also umgekehrt: 


5°, 


Doch wollen wir uns damit nicht zufrieden geben, sondern die direkte o Ableitung dieser a 
chungen kurz skizzieren: - 


Die Ausgangsgl. (24a) wird mit rm 


ost wi ur [- az 
RL Pt 


"Das wird mittels der Substitution 9 = t— 1, ana zu 


Fo) = 


-Damit ergibt sich für das Spektrum sofort: 
Ko) = |F’(w 


Gefragt ist nun nach dem Erwartungswert von |S]? — P> eriey Das ist 


anschauliche Größe. Die Summe stellt nämlich eine Folge ER N Vektoren N Lä 
wegen der ‚Voraussetzung, daß die t; statistisch verteilt sind, ganz unregelmäßig ı 
tungen weisen. Der Erwartungswert für das Abstandsquadrat zwischen Anfa 
dieses. N-seitigen offenen Polygonzuges ist aber bekanntlich Eis! ‚Be 
schon die oben gegebene, Formel gewonnen, da der Grenzübergang sof: 
dings haben wir nur K“ '(o) erhalten, erkennen aber sofort, wa 
haben, denn die obige Überlegung galt natürlich nur für endliches 
Vektoren alle gleiche Richtung und folglich ist S( (0) = N. Dept 
ausführbar, vielmehr muß man setzen! 


Morphologisch homogene Fımkti 


KOagaıı...2r.. Ve : Eye: R . (80) 


KAid)=P=e: 


— 


riet. 
Re. 
= 
Zu 
S 
bi 
2 
S 
2 


2 ee f 
ee... Ve 
[KO (0) ao | 


[Ko )d®d Se 
= a EOS; en 


. N 


KO (0) & ce Klo wre 2 (32) 
= — © RR er 
Ko) =-P = ( [sw dr) — 2702: k(0) . (32a), 
fi a, s 
Pr 
= x . (83 
27 K”® (0) j = ; 
R Karo) Zar | Sr ae 
ku) = gr 9a) ee 


—iot 


i multipliziert und integriert: 


T a 
co - 4 a 5 


—o© 


io) Se RE 


ae IS: ee 


me jo Be 


A ie] 


je 


Flo Ist ea 2 


Fo) = 


SrZ es Math: Mech, 
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2 Wir haben beide Rechnungen nebeneinandergestll, um se Moment klarer 


 hervortreten zu lassen. Es besteht darin, daß in K“ ’(7) bzw. Kr '(&), die doch in ihrem Wesen 
BE quadratische Größen sind, ce linear und nicht quadratisch auftritt. Dies ist ganz allgemein das 
= Zeichen der statistischen Verteilung. Sie bewirkt in der ersten Rechnung, daß aus den N? 


Gliedern der Doppelsumme in Gl. (26) nur die N Glieder mit ,—t; einen Beitrag zu K'® (r) liefern’ 


und entsprechend in der zweiten Rechnung, daß die geometrische Summe der N Vektoren doch 
nur die Größenordnung YN ergibt. Die zweite Rechnung zeigt uns noch ein weiteres: Sie recht- 
fertigt nämlich unsere Annahme, daß im Grenzfall N—oo E{K’(r)} in K(r) übergeht, d.h. daß 
die Streuung dann verschwindet. Denn ‚diese ist wesentlich en durch die bekannte Be- 


ziehung 
Str{f 3 De ......070.% 86). 
2 Die relative Streuung wird damit: = 
Are Str{]2]} 1 
= rleP 02 


Diese Borel erlaubt auch, in praktischen Rechnungen, wo die Funktion ja stets nur in 
einem endlichen Intervall bekannt ist, die Unsicherheit A des dabei gewonnenen Ergebnisses 
abzuschätzen. Sie beträgt für den Einzelwert: 

| rel 
VN Ver 
wo T das zur Analyse benutzte Intervall kennzeichnet!?). Durch geschickte Mittelung über 


‘den Gesamtverlauf von K(r) bzw. K(o) wird die Genauigkeit jedoch meist noch vergrößert 
werden können. 


...(37), 


$3. Aus verschiedengestaltigen Elementfunktionen aufgebaute mhF 
Nachdem wir im vorigen $ den einfachsten Typ der mhF eingehend behandelt haben, ist 
es nun verhältnismäßig leicht, von da aus zu komplizierteren Funktionen fortzuschreiten. Wir 
betrachten eine endliche Folge mhF F,(t) mit den zugeordneten Konzentrationen c;, welche 
alle vom oben behandelten Typ sind, jedoch so, daß die Elementfunktionen der zu verschiedenem 7 
gehörigen Funktionen verschiedengestaltig sind. Daraus bilden wir nun durch Summation eine 
neue mhF 


N EG EN RA ELIE 
die also insofern von einem allgemeineren Typ ist, als sie N verschiedene Arten von Element- 


funktionen enthält. Um den Produktdurchschnitt zu berechnen, zerlegen wir die einzelnen 
Summanden in folgender Weise: 


Fe) =4AF()+F,. 2 
Damit wird 
= # . 1 T/2 | 
5 I = ers ren | 
! (35) HER) + 3 Im [| ARrWArı— di 
ek %7T >00 mp2 
kr je k 


Das letzte Glied verschwindet aber in der Grenze, da es eine Summe von Produktdurchschnitten 
‚statistisch unabhängiger Funktionen vom Mittelwert Null darstellt. Es folgt das Ergebnis, 


55 (K® (7 2 ist Produktdurchschnitt von AF,, also 26) von er 


r a = IE ()+P2= Zah > Er | 
ee 
N ne PX k;(?) 
E 6; und - ka-.— 2 (40) 
I=1 2 } PK 
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oe KO Er.. | 
Aus der Kenntnis von K®(r) und K(P)(r) =F” können wir hier allerdings c und 1 damik“ ke) 
nicht mehr exakt bestimmen. Mit der Be s; [s. Gl. (23)] folgt ne 


Diesen können wir genau so in der abgekürzten FOrE der-El. Se schreiben, nur daß hier die. 
Größen C und k(r) formal anders definiert sind: 


> GS; ne | 
il . . / m ne I 
die Beziehung’ Me 2 b 
. Fr 5, j $ 
C —(1 TESTS Be ana ee Teer ee (41), > 
[K®9 (8) dd E 
wo o die relative Streuung der s; bedeutet: ; \ 
Ng\e \ Fa 
(8) | = 
Da dies stets eine positive Zahl ist, gilt die Abschätzung r nr 
(6) > Se ee (6 ee e El jene 0% (42). a 
[EOO)ad en | a N 
—o N, i 
C* heiße die „effektive Konzentration“. 2 
Br Denken wir uns in Gl. (39) die Zahl der Summenglieder zunehmend und die c; entsprechend 2: 
uf kleiner werdend, dann kommen wir, wenn wir die f;(t) durch eine in willkürlicher Weise von n Ber 
le Parametern x, abhängige Enukon ersetzen, zu dem n- -fachen Integralausdruck: E 
® = 18% ;  Jet@ N day "da,+P REEUE a BE 
Br Ey Fu Et 
hr; 


= f--- Felar-- 2) day:--dam ee ee a er a (44). 


7 I [ea kn m dar: -den REIIER . (4a). 


ur ER PER. Por abe mark 
TE ALT ia 9 © Lech yaäe ıpr 


Die Schreibweise (39a) behält auch dann ihren Sinn, wenn wir statt endlich vieler Parameter, 

die wir anschaulich als n-dimensionalen Vektor auffassen können, deren unendlich vile-be- 
‚ trachten (d.h. zu Vektoren des Hilbert-Raumes übergehen). Die Menge der Vektoren ds 

Hilbert-Raumes ist aber der Menge der quadratisch integrierbaren ! aan 

so daß wir diese damit vollständig erfaßt haben. 


Als ein Beispiel, an dem der Mechanismus der Mittelung besonders dere wird, sei as . { = 
jenige. betrachtet, wo die Elementfunktionen nur von den beiden Parametern ‚Amplitude‘ BREI, 
und „Lebensdauer“ n abhängen: ag 3 EEE he FE IT 


k(&, y=&n (=) = ee 


54 


Damit erhält Gl. (43) die Form: 


5 ie Kia) Bi. el 


—o 00 


Für kleine x werden also hier besonders die Elementfunktio ' 
spielen, während ja ins Integral über KO (?), Herget Ba i 
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Die Ableitung der den Gin. (39a) und (42) entsprechenden Formeln für das Spektrum ist evident, 
nnd es sollen deshalb nur die Ergebnisse angeschrieben werden: 


KO%o)=Cko,; KOW)=-RM ........ (48; 48a). 
F? 
a . kehrt: en 2, 2 ee LA 
27 K® (0) 


$4. Elementfunktionen mit eingeprägten Perioden 
Wir kommen zu einem weiteren Typ mhF, wenn wir Elementfunktionen betrachten, 
denen eine periodische Funktion mit der Periodenlänge © aufgeprägt ist, deren Phasen o un- 
abhängig von den t,; der Elementfunktionen statistisch verteilt sind: 
Bo 0 FB... 2-2 (50) 
Ihr Pröduktintegral hat dann im Einzelfall die Form: 
= Sf Ftt— 9) Pit—g) Pi—p— 1) dt. 


—ox 


Wenn wir dies über alle Phasen mitteln, wobei wir die Substitution 9 = t— o einführen, er- 
halten wir: 


0/2 x s 
1 a 
ha)=7 [rw r@-nar- | ro Ma VE ER ERDG UL RY 
— 9/2 ee 


Der Produktdurchschnitt setzt sich also in derselben Weise wie die Elementfunktionen multipli- 
kativ aus den beiden Bestandteilen zusammen, nur, daß die Phase der periodischen Funktion 


selbstverständlich nicht mehr auftritt. 
Wir wollen die analoge Beziehung für das Spektrum aus der obigen kurz herleiten: 


Sr 2 [A Fe FEN U ; 5 
- 2.65 ") | enteo —0)do 
= ee Pe kr 2° 
kr(®) = > A n)R, (0 f2) „) FE EEE ER ET ee (32) 


y { o m £ 2n. j j 3 
Wenn die Halbwertsbreite von 4,(o) klein gegen = ist, so ist das ein ‚‚virtuelles‘‘ Linien- 


spektrum, d.h. wir haben keine exakt scharfen Linien, wohl aber an den entsprechenden Stellen 
eine sehr starke, in einem schmalen &-Intervall konzentrierte Kontinuumsintensität, während 
sie überall sonst nahezu verschwindet22). Die Ausdehnung obiger Überlegung auf verschieden- 
gestaltige Elementfunktionen ist trivial und führt nur zur Ersetzung obiger Größen durch ent- 
. sprechende Mittelwerte. 

Es gibt noch einen ganz anderen Typ von mhF, der auch ein virtuelles Linienspektrum hat, 
nämlich eine periodische Funktion mit unregelmäßigen Phasenstörungen, (d.h. die Stellen oder 
die Größen der Phasensprünge — oder alle beide — müssen statistisch verteilt sein). Eine sehr 
große Mannigfaltigkeit ist hier denkbar: Man kann sich die Phasenänderung unstetig oder 
auch kontinuierlich vorstellen, die Häufigkeit der Sprungstellen kann eine Funktion ihrer Größe 
sein, usw. Wir. wollen indessen hier diese Funktionen nicht weiter betrachten. Sie werden in 
der Theorie der Stoßdämpfung in der Optik eingehend behandelt. 

Noch einige Bemerkungen zu mhF mit echten Spektrallinien: Die einfachste Form ist die, 
daß zu einer mhF der bisher behandelten Art additiv eine periodische Funktion hinzugefügt 
wird. Dann ‚überlagern sich die Spektren beider Funktionen auch additiv. Ähnlich liegt der 

Fall, wenn ein periodischer und ein morphologisch homogener Bestandteil sich multiplikativ 


en 
18) Für fastperiodische Funktionen erhalten wir genau das gleiche Ergebnis, nur daß, wie immer, @n 


VCH de: 
an.die Stelle von 9" tritt. 2 


(tionen und ihre Erzengung zu.: »0 A: 


> Fi“ = _ 1 7Z 
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zusammensetzen. Allerdings hängen dann beide Spektren komplizierter zusammen. 3 der er 
Praxis wird es bei solchen Funktionstypen stets vorteilhaft sein, die periodischen Anteile?#) erst” an 
abzutrennen und dann hinterher das Spektrum bzw. den Produktdurchsehnitt zu bilden. a 
Eine weitere Art gemischter Spek werden, wir im nächsten Abschnitt gleich ‚kennen- Er 
lernen. E $ 
$5. 'MhF mit veränderlicher Dichte der Element 


Die Klasse der bisher betrachteten mhF kann man zusammenfassend dadurch kenn- de. 
zeichnen, daß für sie die Dichte der einzelnen Elementfunktionen eine Konstante (nämlich die R 
Konzentration c bzw. C) ist. Wir kommen zu einer größeren Mannigfaltigkeit, wenn wir Funk- at 
tionen zulassen, deren Dichte veränderlich ist. Natürlich kommt nur eine solche Veränderic- 
keit in Frage, die den Charakter der mhF erhält, d.h. die Dichte D(t;) der einzelnen Element- » 
funktionen muß selbst auch eine mhF sein. Außerdem muß sie selbstverständlich der Bedingung 
genügen, daß ihr Mittelwert mit der Konzentration übereinstimmt. Wir behandeln zuerst den 
einfachsten Fall, nämlich daß die Dichte eine periodische Funktion ist, unter der weiteren Voraus- 
setzung der Gleichgestaltigkeit aller Elemente. Dann ergibt die Rechnung des $2 für das Inter- 


vall (—T/2, T2): a 
% 1X BE a SE 
K cl SZ & \ N 
ja ee : 
Aber die Wahrscheinlichkeit der Elementabszissen ist nun: 4 
be} 1. Pie) KR 
Wii) — immer. (53), 
a & lo außerhalb IE 3 
} wo P(t,) eine periodische Funktion vom Mittelwert e bedeutet. Damit folgt: 
; Tea 72 } 4 
E{K (?)} =(C* ky(?) 2 j [ri (1) Pt.) ky(t a t; + tz.) di; di;. a i 1 3 
; ar — 2 —Tj2 s ern i 
so Das geht in der Grenze für große 7 über in ei | 3, | | 
Be, Kd=ckkd)+f Kr —OMlöde N. Be 


Das letzte Glied stellt im wesentlichen den Produktdurchschnilt Kr) En; ee: x 


ie EEE, 
: Di) = Pia RER: 
5 dar, jedoch mittels der Gewichtsfunktion %y(®) ausgeglichen. Mit gr Abkürzung X | 
; J[Krr—9)k (0) dd = [Kp(t re ERE oo 
Br - haben wir so das Ergebnis: E 8 r Er 
; r h K(r 3) =6(° hut 7) + u: al. ne .. ER ne SER .64n).. 
Fi Der einzige Unterschied gegenüber der Gl. (29a) besteht darin, daß K) x nun ı eine echte 
periodische Funktion und nicht wie dort eine Konstante ist. "F- Hr. „2 St 
Be» Die Erweiterung der Überlegung auf auf Funktionen mit verschiedenges altige Elemer en 
D- ist trivial: Man muß nur %,(z) durch &,(r) ersetzen. Ist auch die Dichtefunktion für die ver 2 
schiedenen Elementarten verschieden, so muß Kp(r) ebenfalls durch N Mittel $ 
Ei ersetzt werden. Dabei kann unter Umständen der periodische ‚Charakter yerl oren ‚schen. BT: 
Be Wir wollen noch die Gl. (54a) entsprechenden Gleichungen für das mans | 
Ko) =c- ky(o) SE 
Ki (0) = [Kr (@n)]y: = 2% Klon R: 
Man erkennt, wie die Spektrallinien K®(o,) durch den Faktor 
t dem & gegen Null geht, entsprechend geschwächt werden, doc 


kommenden Fällen für alle Linien mit wesentlicher‘ Intensitä on 
breite von Ty(o ) sein, so daß dieser Be nur eine gering! 9 = 2 


j tg en 
u; n \ ix we N 4? 
w L 4° Pr e # 


- 
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N alle no tie sind und A auch die Elementfunktionen 9(t;) der Dichtefunktion. Für 


: das Näherungsintervall (— 7/2, 7/2) ist dann die Wahrscheinlichkeit der t;: 
I; j R 
br: “ . 1 
EN 3 Wi, = — gl te) SE 0 (57), 
\ e-1 b 
Pt 1 
he | vwW=|T im Intervall (— 7/2, 72) (57a), 
Ss 0 außerhalb 
ü mit der Nebenbedingung: ’ 
St | T/2 
: T | a a a (57b). 
Y 4 | 
r Damit erhälten wir für den Erwartungswert: 
as N j N R Tj2 Ba A 
Fa Ti 153 aa | j [soo 00 — 1) kir—t; +1) di, di.dt, 
Bart ee e=1 _Ty2 —Tj2 —Tj2 
\ N: ze Te ma Ta 1 
ae > 33 | j j gl; — to) gie +to) Kir— y-+ tb) diz dig dt, dte. 
pe ee —Tj2 —Tj2 —7j2 —'Lj2 
Das ergibt mit y 
R f 
’ Mr T en Na A EN ee ar Re (98), 
wie immer im Grenzfall großer T': 
| K(ı) =c kr) + [ho EHEN NN EN . (59). 
| Das zweite Glied stellt wieder einen mittels der %,(0) ausgeglichenen Produktdurchschnitt dar: 
er SR 7 — 0) k,(d) dd — -[ruto) 909 dh. 2. (60) 
1 | so daß wir die Endformel auch schreiben können: 
vr \ ’ 
ern 7° Kid. te) hr Urt Ay + BE SEE NR RAN EEE NEN NEE 
e ER die en (7b) wird in Analogie zu Gl. (31) y Berueledt: 
Ri 752 
hr Ä nen 
= [.K®ı0) do 


wo K%(d) den zweiten Summanden von K(r) bedeutet. Für den Übergang zu Funktionen 
mit, verschiedengestaltigen Elementfunktionen f(t) und g(t;) gilt das gleiche, was oben für Funk- 
tionen mit periodischer Dichtefunktion gesagt 


Ken hatten in der vorhergehenden Rechnung für die Dichtefunktion 


EAN: ent EI KR DLG)=e, 
"% 


0=—&x© 


H) 


ö 2 also eine , mhF vom einfachsten Typ gesetzt. Wir können aber auch ebenso eine solche vom 
ige ‚oben IE Aue Typ einsetzen, d.h. 


’ an Kr 2(t) = 3 h(t«— to), D;(te) =t, 
Dr ns jr N 0=—»o 

önnen Kieo Re noch beliebig oft (allerdings nur endlich oft) wiederholen. 
für die Dichtefunktion höchster Ordnung Dy (t;y_,) noch eine periodische Funktion 
ehmen, so erhalten wir für den Produktdurchschnitt der allgemeinsten mhF der bisher , 
ten Klasse: 2 wi ; 


un 2 uor + NN 62). 


IH a de TRENNT IN % 
- i Be N a. Eat 
\ N . x ’ r4 
f # 


_ 
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Die Abkürzungen bedeuten hier: “ } 
ve? 0 [k(r)]o = k(?); 


#0: [k(r))» = Sin (T— 9) [k(d))-ı dd, (62a), 


Kr)lv= J Krtr— 9) [ka 40 \ 


wobei die Größen mit dem Index 0 die eigentlichen Elemente f(t), diejenigen mit dem Index 
v=1...(N—1) die Elemente f,(t;,_,) der Dichtefunktionen »-ter Ordnung kennzeichnen. 


Für das Spektrum gilt entsprechend: 


4 AR N—1 
KONVERNCHROERD En N (0 
l yv=0 
K») (©) ii \Np(on’in 5 t. ... . . . . . . . . . (63a). 
Für die „effektiven Konzentrationen‘ C% erhalten wir in voller Analogie zu den Gl. (42) und (49): 
F Pr 


O* — . (64). 


eo N a re 70) 
[E®"oas 2% 0) 


Die Mannigfaitigkeit der mit diesen Formeln umfaßten Funktionen ist so groß, daß wir uns 
fragen wollen, was diese Mannigfaltigkeit allgemein kennzeichnet. Dazu gehen wir noch einmal 
auf die einfachsten in $2 behandelten mhF zurück. Für sie war die Wahrscheinlichkeit aller 
Elementfunktionen im Intervall (— 7/2, T/2) die gleiche Größe W(t,) =1/T. Darin liegt ein 
doppeltes: Jede Stelle i der Abszissenachse ist gleichberechtigt, und die einzelnen Element- 
funktionen sind völlig unabhängig voneinander. Bei den Funktionen von $ 3 blieben diese beiden 
Eigenschaften für die einzelnen Arten von Elementfunktionen erhalten. In diesem $ schließlich 
wurde”die Gleichberechtigung aller Abszissenstellen aufgegeben, nicht aber die zweite der eben 
genannten Eigenschaften, so daß wir zusammenfassend sagen können: Die von uns bisher er- 
faßte Funktionenmannigfaltigkeit ist der Klasse aller mhF, die aus statistisch unab- 
hängigen Elementfunktionen aufgebaut sind, isomorph. 


$6. MhF mit nicht statistisch unabhängigen Elemenffunktionen 


Die letzte Erweiterung, die wir vornehmen wollen, besteht darin, daß wir nun auch noch 
die Bedingung der ‚„‚Unabhängigkeit‘‘ der Elementfunktionen fallen lassen. Zuerst müssen wir 
diesen Begriff mathematisch präzisieren und gehen dazu nochmaıs auf die in $2 behandelten 
Funktionen zurück. Numerieren wir ihre Elemente derart, daß t;_, <t4y<tj,, (was wir in 
der Rechnung dort nicht getan hatten!), so folgt aus der da gestellten Bedingung (25) für die 
Wahrscheinlichkeit, daß der Abstand zweier aufeinanderfolgender Elementfunktionen zwischen 
Ö und ®-+dd liegt, die bekannte Formel: 


WIO)dO = er dlod)  ee eles) 


Aus ihr leiten wir leicht auch die Wahrscheinlichkeit für den Abstand der j-ten und (+. n)-ten 
Elementfunktion ab: 


0-5 07] : : 
» 0-0, 1 n- 
WB Tr We) WO) Wed) W0— 3 9)d0,-- Add) (66). 
0.0 ö { el 
Einsetzen von Gl. (65) ergibt dann sofort: ; i a 
EN (cd)rA1 Re Sir 
W„(9) dd = et EACH). 5 


Diese Verteilung ist das Kriterium der statistischen Unabhängigkeit. Wir bezeichnen also de . 
Elementfunktionen unserer Funktion als nicht statistisch unabhängig, wenn ihre Verteilung 
von der hierdurch gegebenen ‚„‚Normalverteilung‘ abweicht. Wir wollen nun den Produktdurch- 
schnitt einer mhF berechnen, die aus lauter gleichgestalligen Elementen aufgebaut ist, deren 
W„(ö) jedoch willkürlich vorgegeben sind, abgesehen von einer einzigen, später zu formulierenden 


> 


1°) Die Formel kommt folgendermaßen zustande: Das Produkt bezeichnet die Wahrscheinlichkeit, daB & ine 
der Abstand zwischen der j-ten und (j + 1)-ten Elementfunktion ®,, zwischen der (j + 1)-tenund  +2)-ten 
Elementfunktion 9, beträgt, usw. Über alle möglichen Lagen der Zwischenpunkte wird dann integriert. 


3» a) ee ee a er 1 7 ah Ein au Are ne um 2 
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Bedingung, die aber für jede vernünftige Verteilung von selbst erfüllt ist. Wir gehen von Gl. (26) 
in $2 aus und ordnen die Glieder der Doppelsumme folgendermaßen um: 
1 Non 


Eid) =ec:ktt DPA, + imn m) Hklr— man tim) 0.» (68), 


=l ml 


Wir können hiervon die mathematische Erwartung bilden, indem wir mit W,„(®) multiplizieren 
und über alle positiven ® integrieren, da ja die i-Werte nur in der Verbindung 9 = tyın — Ims 
‚nicht aber unabhängig voneinander auftreten: 


. N—1 Ion 
BR Sek + > 5 mo RO) +kT— dd . ... (68a). 


> an 2, m 


Man kann die Summation über », welches ja unter dem Integral nicht mehr vorkommt, aus- 


en und erhält in der Grenze großer T: - 
Kir) =c-k(r) ef B3 Wr) iktt +9) +k(r— 9) dd 
0 n= = 


—=c:k(r) +e[ > Ww.\d)) k(t— ®) dd. 


—o n=1 
Ist W„(9) die Normalverteilung, so bekommt die Summe unter dem Integral den Wert cund dieses 
nimmt dann in selbstverständlicher Übereinstimmung mit Gl. (29) von $2 die Gestalt 


c?- Sr k(d Ma = BE 
an. Für den allgemeinen-Fäll führen wir die „statistische Anomalie‘ 


SR a oe Ra ae ee (60) 


nl 
ein. Diese Größe soll eine integrierbare Funktion sein!®). Das ist gleichbedeutend damit, daß 


“ 


oo “ 
& W„'0) für alle positiven © existiert und in der Grenze großer ® gegen e (das ist die Summe 
N £ 
der Normalvertleilungen (67)) strebt. 


"Dann nimmt unsere Gleichung die endgültige Gestalt an: 


: Kin) = elktr )+ Ju) k (a0) +7 ER OFT NT 


Wir wollen das hier neu auftretende Integralglied diskulieren. Es bewirkt eine Verzerrung von 
K®(r) derart, daß k(r) und K®(r) nun nicht mehr von gleicher Form sind. Ist die Verteilung 
„übernormal‘“ (d.h. im Bild: ‚‚Die Elementfunktionen ziehen sich an‘), so ist «(®) positiv und 
K'®(r) wird vergrößert, bei „‚unternormaler‘‘ Verteilung (‚Die Elementfunktionen stoßen sich 
ab“) ist #(9) negativ, K'®(r) wird verkleinert. (9) kann aber auch abwechselnd positiv und | 
negativ sein. Wir können diesen Tatbestand auch so beschreiben : Die Anomalie bewirkt, daß sich 


der durch # gegebenen Geraden symmetrisch zum Nullpunkt der r-Achse ein „Hügel“ (wenn 

u(d) positiv ist), eine „Mulde“ (wenn «(d) negativ ist) oder gar eine Folge von Erhöhungen und 
Vertiefungen bildet, (wenn «(d) um Null herum schwankt), so daß «+ k(r) nicht mehr wie früher 

einer Geraden, sondern dieser unter Umständen sehr komplizierten Kurve überlagert ist. 


Wenn wir das ‚„‚Anomalieperiodogramm‘‘ 1?) 


oo 


22 # 2 u (w) [io 2 N RE ie en zZ) 
24 k % 2 % Ei Y2a 

einführen, so nimmt die der Gl. (70) analoge Beziehung für das N die einfache Gestalt 

| eh KO (w )=c{1 +27 Ro)! kim ER 2 rn nen. (72), 

$ ONE 3 RE at. (Ta) 


an. Führen wir noch” zur Abkürzung das „reduzierte Pröduklintegral 5; 


In ac. pie) Nicht, zmendig quadratisch sehän: a(d) kann durchaus auch ö-Funktions-Bestandteile ent- 
a: Ist, 
e B Wald) = On —9), also. u(d) = Ss (On — dl) — ec, 

Luk n=1 

ten wir eine rein periodische Funktion den Periode ©. 

Wir dürfen nicht ‚Spektrum‘ sagen, da dieses seinem Wesen nach eine nichtnegative Denn ist. 


2 1 
nn aber sehr wohl negätive Werte annehmen, wenn auch nie kleiner als — 5 werden, da sonst K(w (®) 
werden apaßte,; was aber einen Perssaeruöh u. v2 


12* 
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bzw. das „reduzierte Spektrum‘ 
z(o) =-{1+NRnaw)}klo)..:..: 2... a . (73a) 


ein, so können wir damit formal unsere Gln. (70) und (72) auf dieselbe Form bringen wie die 
Gleichungen (29a) und (32) von $2 und können dann auch die Erweiterungen des $3 und 5 
genau so vornehmen. An diese letztgenannte Erweiterung allerdings wollen wir noch einige 
Bemerkungen anknüpfen: GR 
Ersetzen wir c durch die variable Dichtefunktion D(t), so gilt für solche Funktionen de 
Normalverteilung nicht mehr in der einfachen Gestalt der Gl. (65), sondern wir haben statt dessen a 
die allgemeinere Verteilung für den Abstand des Nachbarelements von einem Eu an der 


Stelle £: e 
-/ par RE h \ 
| We)dd=Di+We RE N ©...) “ 
und entsprechend für den Abstand des n-ten Elements: | 2 
. t+9 
Ba Fa oh Em t+9 r 
Er, W„(9) dd = Dit 4 d) e [ J fd D(dn-1) 29. ddna dd 
I [ ‚t 
f {+8 n1 BR (75). 
r (J D(t) ai) -7 par 18) 
R Fur Ü 
| er An e 


Diese Formeln sind allerdings mit den Formeln (65) bzw. (67) solange praktisch gleich, wie n/e 
klein gegen die Größe des t-Interyalls ist, in denen sich D(t) merklich ändert. _ 

Wir können damit nun leicht die Anteile des Produktdurchschnitts, die von den D(t)- 
Schwankungen herrühren, als Wirkung einer statistischen Anomalie ’auffassen. Allerdings 
müssen wir diese dann noch allgemeiner definieren, denn, indem wir Wo ) als Funktion von t 
auffassen, müssen wir in Gl. (68a) eine weitere Mittelung durchführen, und dann ergibt sich 
an Stelle von Gl. (69): z 


72 E FR iR 


1 u ‘ ERpA: 
u) = im 5 [DW Zwalha—e. 2.2.2... 0). 
2 Ne Ta en } HA 
Hier ist es nicht mehr vernünftig, die Forderung, daß «(®) integrierbar sein muß, aufrecht zu 
erhalten, sondern man wird zulassen, daß diese Funktion auch periodisch um Null schwanken A 
kann. Das Anomalieperiodogramm wird in diesem Fall auch Linien enthalten, sonst hat dies 


‚aber keine weiteren Folgen. In diese Formel setzen wir nun unsere Sn (75) ein A eT- Biol 
halten dann, da ; TAB RER 2 ET AR 


re 3 Dur a 
n=1 j Up. EEE N Kr RR Me 
für die statistische Anomalie RA a AB N er 1 Mta Ate ee Ra 
7/2 EX N RE ER AN Dre ar 
I As.” 
u(9) = — lim — | Dit) Du+di-e= Koi. u % 
> C T1T>% “b —7]2 KLraee ’ ih 
Das RN aber, in Gl. (70) eingesetzt, URN RR Mi 
Be“. "u #2 LA 
a ; Kd)=c.kd)+ | Ko(0) ke) 0. Be 


| also, i in eine etwas allgemeinere Schreibweise gefaßt, genau das Erg 
Diese Erkenntnis bedeutet zwar an sich, abgesehen von der allgem 
Anomalie, keinen besonderen Fortschritt, aber sie zeigt doch den. Be Z 
jener einfacheren Betrachtungsweise in $5 und der hier verwendeten n viel 
haben wir aber nun die größte Klasse von Funktionen, die Obeihaune pt 
Beschreibung zugänglich sind, grundsätzlich erfaßt, BE 


& 


| 
R 
5 


12) Die Ableitung geht hier genau so wie im Fall konstanter "Wahr 


des n-fachen Integrals in die n-te Potenz eines gewöhnlichen Integrals 
Integration. ur 


Eingegangen am 30. Juli 1949. 
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Fließgesetze normal-zäher Stoffe im Rohr 


Ein Beitrag zur Rheologie 
Von L. Prandtl, unter Mitwirkung von Fr. Vandrey!) 


Unter ‚‚Normalzähen Stoffen‘ werden hier Stoffe verstanden, deren Fließgesetz y = f(r) durch die 
Formel f(r) = C Sint/A angenähert werden kann, mit » = Deformationsgeschwindigkeit und x = Schub- 
spannung. Aufgabe der Abhandlung ist die Berechnung der Ahbängigkeit der Durchflußmenge eines solchen 
Stoffes in einem geraden kreiszylindrischen Rohr vom Druckgefälle, ferner die Errechnung der Konstanten 
A und C vorgelegter Stoffe aus den Messungsergebnissen. 


For materials called here,,normal-viscous‘‘ the law of flow, y = f(r) is described approximately by the 
formula f(r) =C sinh z/A (y = rate of strain, t = shearing stress). Considering materials of such kind 
ihe paper aims to calculate the rate of flow ihrough a straight circular pipe for given pressure gradients. 
Aiso a-method is indicated how to compute the constants A and C for given materials [rom measuremenis. 


Des substances, dont la loi de flux y = f(r) peut ötre approwimee par f(r) = C Sin r/A, sont denommees 
„de viscosite normale‘‘; y vitesse de la deformation, 7 tension tangentielle. Le pr£sent essai a pour but de 
calculer pour une telle substance la quantite d’Ecoulement a travers une tube droite eylindrique dans sa de- 
pendance du gradient de la pression. Une methode pour determiner les constantes A ei C des mesures 
pour des substances quelconques est demoniree. 

Ilon ‚‚HOPMalbHo-BASKUMH MATePmanamu‘ 31eCh IIOHHMAITCH TAKHE, YCHOBHe TEKYUecTH 
KOTOPEHX 9 =(T) MOxeT ÖBTb IPHÖNNKeHHO BbIpaskeHo dopmyaoi [(T)=Cshr/A, rne 
Yy = CKOPOCTB ehopmanun, a 7 = cnBuramımee Hampsekenme. lammoe uecmenoBanne cTaBuT 
cebe 3anayeii OIPEeNeIHTb 3ABUCHMOCTB TIPOTERAWMINETO KOAMYECTBA TAKOTO BEIMECTBA OT 
rpanueHTa NaBleHuA IpHu IpoTekaHunn B IPAMOÄ NMAIHHNPHYecKofi kpyrıoi TpyDbe. ApoMme 
TOTO BEIYHCHAOTCH MIA MAHHEIX MATEPHAJIOB HOCTOAHHBE A U C HC pesyIbTaToB n3Mepenni. 


1. Einleitung 


Die theoretischen Betrachtungen des Verfassers von 1928 [1] über die Häufigkeit von 
Platzwechseln der Molekeln in einem unter Spannung stehenden festen oder zähflüssigen Körper 
haben zu Beziehungen geführt, die sich für konstante Fließgeschwindigkeit y, die im folgenden 
mit c bezeichnet werden soll, durch die Formel 


BE SEN TEA EN EN RE: BI De Vet en ei hate 


annähern lassen. C und A sind dabei Stoffkonstanten. Soweit Temperaturunterschiede in 

Frage kommen, sind beide als Temperaturfunktionen zu betrachten. Formel (1) ist von A. Nädai 

1939 an denlangsamen Kriechbewegungen 

von: belasteten Metallstäben bei hohen 

Temperaturen rechtbefriedigend bestätigt 

worden [2] 2). 

Es ist zu erwarten, daß dasselbe 

Verhalten auch bei stark zähen Flüssig- 

keiten auftritt, soweit die Struktur dieser 
$ Substanzen genügend stabil ist, um nicht 
2 durch zeitliche Vorgänge oder auch durch 
den Meßvorgang selbst verändert zu 
werden. (Bei den Untersuchungen der 
„Rheologen‘ laufen häufig derartige Vor-, 
gänge nebenher). Stoffe der gekenn- 
zeichneten Art sollen hier mit dem 
22 - Namen „‚normal-zäh‘ gekennzeichnet 
e werden. Die gewöhnliche zähe Flüssig- 
keit (,Newtonsche Flüssigkeit‘ mit 
c = t/n) entspricht in Gleichung (1) dem 
4 - Grenzfall 7 < A. Die gewöhnliche Zähig- Bild1. Schubspannung r abhängig von der Gleitgeschwindigkeit e 
keit’ nist hier = A/C. | 
Be - Bild 1 zeigt ein Schaubild für Gleichung (1). Gemäß der in der Festigkeitslehre üblichen 
_ Darstellungsform ist die Umkehrfunktion. von Gl. (1), also 7 = A X Ein c/O = Aln [e/C 


+Y1 + (e/0)?] wiedergegeben, und zwar in.drei Maßstäben von c/C, um den ausgedehnten Zahlen- 


er 1) Der Verfasser verdankt Herrn Vandrey, der sich ihm zwischen der Beendigung seiner Göttinger 
- Tätigkeit und seiner Übersiedlung nach England auf zwei Wochen zur Verfügung gestellt hat, die Reihen- 

_ entwicklungen (4) und (5) und die numerischen Rechnungen hierzu, besonders aber die Auffindung der Lösung 

der Aufgabe des 4. Abschnitts und die graphischen Darstellungen der Abbildungen 4 und 7. 

2) Über die hauptsächlichsten Gedanken seiner Arbeit von 1928 hat der Verfasser kürzlich in den „Physi- 
lischen Blättern‘‘ berichtet und dabei besonders rheologische Zusammenhänge berührt [3]. Er hat sich auch 

ort ausdrücklich auf stabile Stoffe beschränkt. | 

SE Gh 
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bereich dieser Größe zur Anschauung zu N Die gestrichelt u; Anfangstangente 
stellt das Newtonsche Verhalten dar. 

Die gegenwärtige Abhandlung ist der speziellen Aufgabe N. die Formeln zu ent- 
wickeln, die an Stelle des Hagen-Poiseuilleschen Gesetzes zu treten haben, wenn es sich 
um solche stabile ‚„‚Nicht-New tonsche Flüssigkeiten‘ handelt, für deren Fließgesetz die Gl. (1) 
eine brauchbare Näherung liefert (Es ist zu erwarten, daß amorphe Stoffe, wie Pech, geschmol- 
zenes Glas und ähnliche, "besonders typische Vertreter dieser Gruppe sind). Um einen kurzen 
Namen für solche Stoffe zu haben, sollen sie hier „normalzähe Stoffe“ genannt werden. 

Die erste Aufgabe wird in der Ermittlung der Geschwindigkeitsverteilung bestehen, die 
bei gegebenem konstantem Druckgefälle für die Längeneinheit, DR = = PG2 einem Rohr vom 
Innenradius R auftritt. Im Zusammenhang damit ergibt sich leicht die Durchflußmenge Q für 
die Zeiteinheit, die beim Versuch gemessen werden kann. Da diese hier nicht mehr proportional 

' mit dem Druckgelälle ist, ergibt sich keine konstante „‚Zähigkeit 7°“ mehr. Wir wollen die Zähig- 
keit nach der Hagen-Poiseuillesche Formel, 7 = R? p'/8 Q, trotzdem als bequeme Rechen- 
größe beibehalten. Je nachdem man diese Pseudozähigkeit abhängig von der Durchflußmenge 
Baer abhängig vom Druckgefälle aufträgt, ergeben sich dann zwei etwas verschiedene Kurven, 
vgl. Bild 5 und 6. r 


2. Gesehwindigkeitsverteilung und Durchflußmenge 
Der Ansatz erfolgt ganz entsprechend dem für Newtonsche Flüssigkeiten: 


Das Kräftegleichgewicht in der unbeschleunigten Strömung verlangt in DE Aue Weise 
eine Schubspannung im Radius 


r-nr][2P) 
Die Geschwindigkeit u ist = 0 für r = R (Haften an der Wand!); sie sei positiv für alle Radien 
r < _R. Dann ist gemäß Gl. (1) 


| ce=—duldr = (0 &in Ba 
Mit der Abkürzung p’/2 A = a ergibt dies integriert 


= 
\ u(r) — [du = Ca (Cs & R— Eos ar) ERRRR I SR IE NE 


e 
Damit kann die Durchflußmenge in der Zeiteinheit Q =2r 1 ur dr gebildet werden. Durch 
partielle Integration ergibt sich u 

(&vs aR FArN 1) . . . . . (3). 


Q=2n Era |, Gos aR— SinaR + ——_ 


“ 

2a&R a 
Für die praktischen Rechnungen ei es sich, Reihenentwicklungen nach der Dimensions“ 

losen «R zu verwenden, Potenzreihen für kleine Werte von «R, Ausdrücke mit e*# und € 

für die größeren Werte, Setzt man der Kürze halber a«R =p'R/24A = a, so ergeben sich nach E 
den V ne yschen Ermittlungen für den Ausdruck in den eckigen Klammern in Gl. (3), der 

"jetzt f(a) heißen soll, die beiden Ausdrücke 


\ n2n+1), , 1 19% 1 1 Ä 3 
a) = van nn et U a aan EN ee LE ea 
fa) Sn } ae (4) 
und Fast, 1 1 1 ANTT f 
— u, 4b R f + ‚ , 
fa) = e E a ee) = +e (7 Kr + EN RE U Sn | 
Für großes a reduziert sich Gl. (5) praktisch auf? Ken Burn: . 
1 1 1 i win BR  . 
1), e® NER NR UN ER HRE  L7, USE 
fa) ze © ap eh REN R (6) Br 
Die Näherungslormel (6) ist bei a =5 schon auf 2 Dezimalen, bei «a = 7 bereits auf 4 Dezimalen BEN 
genau. Für, sehr große Werte von a kann f(a) = e“/4 gesetzt werden Nu 2. Ne 


Bei der numerischen Arbeit empfiehlt es sich, statt‘der Funktion f(a), die mit dem Glied 8 
beginnt, eine neue Funktion (a) = 8 f(a )/a? einzuführen,die beia =0 ‚mit dem Wert 1 beginnt, 


und die daher bei praktischen Aufgaben r Ai: auftritt. „Aus der 


Pl =145 + en 


2. . Math. h. : er a 
Bi. 30 Nr. 5/6 ultfsuni 1950  Prandtl-Vandrey, Fließgesetze normal-zäher Stoffe im Rohr 171 


und aus der asymptotischen Näherungsformel (6) 
| 2 2) 
ne ee N... ; 
#0) =2eieli + OF 


asymptotisch ist 
OS A. hen lea: 


Zahlentafel I gibt die Werte von o(a) für « = p’R/2 A zwischen «a = O und a = 10 wieder. 
Bild 2 zeigt den Verlauf von log p(a) abhängig von log a nach der Tabelle und nach der Nähe- 
rungsformel (8) und (8a). Die letztere ist jedoch wenig genau, sie gibt bei « = 10 noch 18 v.H., 
bei a=20 noch 9,5v.H. zu viel. 


Tafell. o(a) nach Gl. (7) und (8) 


“. \.vetaoelare | .Pta) | 6 Y(a) a A) a ?(a) 
0,00 1,000 1,5 1,272 3,0 2,402 4,5 5,800 7,0 33,81 
0,10 1,001 1,6 1,313 3,1 2,529 4,6 6,185 7,2 39,31 
0,20 1,004 1,7 1,358 32 2,671 4,7 6,599 7,4 45,80 
0,30 1,010 1,8 1,407 3,3 2,821 4,8 7,056 7,6 53,35 
0,40 1,018 1,9 1,459 3,4 2,982 4,9 7,540 7,8 62,31 
0,50 1,028 2,0 1,517 3,5 3,155 5,0 8,061 8,0 72,78 
0,60 1,041 1 1,578 3,6 3,341 5,2 9,234 8,2 85,10 
0,70 1,055 2,2 1,645 3,7 3,541 5,4 10,60 8,4 99,64 
0,80 1,073 2,3 1,716 3,8 3,757 5,6 12,18 8,6 116,7 
0,90 1,093 2,4 1,794 3,9 3,988 5,8 14,02 8,8 136,9 
1,0%" 1,115 2,5 1,873 4,0 4,238 6,0 16,18 9,0 160,6 
1,1 1,141 2,6 1,968 4,1 4,507 6,2 18,69 92 188,6 
1,2 1,169 27 2,065 4,2 4,794 6,4 21,64 9,4 221,6 
1,3 1,200 2,8 2,168 4,3 5,101 6,6 25,09 9,6 260,6 
1,4 1,234 2,9 2,281 4,4 5,438 6,8 29,10 9,8 306,7 
1,5 1,272 3,0 2,402 4,5 5,800 7,0 33,81 10,0 361,3 


3. Versehiedene Einzelausführungen 
Nach Gl. (3) war die Durchflußmenge 
Q:= 2% CR2]a-f(a). 
Mit «=aR ist 1/a = Rla, also 
Q=nCR°®.2fla) 
Bei Einführung von o(a) : a?/8 ei BE wird 
Q=nCR:-ald:p(a) =nCR3-Yyla). . (9) 
oder mit «= p’R/2A auch 
QS3T Re)... . (0). 


Im Grenzfall a—0, wo 9(0) = 1 wird, und A/O =n, 
mit der Newtonschen Zähigkeit identisch wird, geht 


. also Gl. (10), wie es sein muß, in die FormelQ =rzRt 05 KO HIER ZER N 3N 451 72510 


» '/Snoder Hagen-Poiseuilleschen Strömung über. Bild 2. Doppeltlogarithmische Darstellung von (a). Die 
Wenn die Stoffkonstanten unverändert bleiben beiden Näherungen GI. Ed ee SOSE Alad, Bertzlonele 
D eigefüg 
der Druckabfall aber alle Werte von Null an durch- 
‚läuft, so läuft auch a durch alle Werte. Der Verlauf der Durchflußmenge ist dabei aus Gl. (9) 


am besten zu erkennen. Er ist aus Bild 3 zu ersehen, die die Funktion y(a)=ay (a)/4 abhängig 


. von a wiedergibt. Auf die qualitative Ähnlichkeit mit Bild 1, die das Elementargesetz darstellt, 


sei hingewiesen. 
Bild 4 zeigt den Gang des Geschwindigkeitsprofils ion mit a nach Gl. (2). (Die mittlere 


Geschwindigkeit w ist =@/aR?.) Mit wachsendem a werden die bei a— 0 paraboloidförmigen 
Profile immer mehr abgeflacht. Die Grenze @—> & bedeutet den idealplastischen Körper. Der 
Ausfluß von Zahnpasta aus der Tube ist ein Beispiel für sehr hohe « bei sehr kleinem A. 

Die Schubspannung an der Rohrwand, 7,, ergibt sich nach früherem zu c, = p'R/2. Das 
Argument in Gl. (1), z,/4, also das Verhältnis der Wandschubspannung r, zu der für den Stoff 
 charakteristischen Spannung A, ist demnach = p’ R/2 A, also = a. Man erkennt hieraus, warum 
diese Größe im Mittelpunkt unserer Rechnungen steht. 


- *4) Diese Formel mit der Reihe (7) für p(a) findet sich auch bei Vandrey. 


I 
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Wie schon in der Einleitung erwähnt, ist die Zähigkeit, die eine Newtonsche Flüssigkeit 
haben müßte, um bei gleicher Durchflußmenge Q in demselben Rohr denselben Druckabfall 
und nach obigem auch dieselbe 
Wandschubspannung zu liefern wie 
ein konkreter Fall einer Nicht- 
Newtonschen Flüssigkeit, eben- 
falls eine lehrreiche Hilfsgröße. 
Sind C und A die Stoffgrößen 


Bild 3, Verlauf von y (a) (Durchflußmenge abhängig vom Druckgefälle) Bild 4, Geschwindigkeitsprofij u(r)/% 
abhängig von « 


dieser Flüssigkeit, ferner Q die Durchflußmenge und 9’ der Druckabfall, so ist die ‚„‚rechnerische 
Zähigkeit‘“ 
N Fa SEO NET Re RL) 


h Wenn wieder wie oben die Stoffkonstanten A und ( fest bleiben sollen, dann läßt sich 7, 
2m abhängig von a angeben. Mit Q nach Gl. (9) und 9’ =a.2 A/R sowie A/C = n, findet man 


NEENODÄN RATEN TEL ING ED ER 12 


7, nimmt also mit wachsendem a, d.h. wachsendem 9’ zuerst wenig, später immer stärker ab, 
s. Bild 2. Abhängig von @, also von y(a), ergibt sich ein Verlauf nach Bild 5. 
Man sieht hieraus, daß für Veränderlichkeit von n nicht immer eine ‚‚Strukturviskosität‘ 
| vorzuliegen braucht, sondern auch normalzähes Verhalten eine bestimmte Form einer solchen 
| Veränderlichkeit hervorbringen kann. 


Bild 5, n,/n, abhängig von y (a), d.h. von Bild 6, n„/n. abhängig von a, d,h, vom 
der Durchflußmenge Druckgefälle ’ i 


4. Ermittlung der Stoffkonstanten A und © aus Meßergebnissen 


"Die Lösung dieser Aufgabe verdanke ich völlig Herrn Vandrey. Der hier abgedruckte RN 

Text enthält gegenüber den Vandreyschen Rechnungen einige leichte, meist formale Alände- 

rungen, die sich aus inzwischen erreichten Resultaten des Verfassers ergeben haben. a 

y Die Aufgabe sei in der Form gestellt, daß mit dem zu untersuchenden Stoff zwei Messungen 

a gemacht werden, bei denen für zwei Druckabfälle für die Längeneinheit des Rohres pyundp 
j zwei Durchflußmengen in der Zeiteinheit Q, und @, festgestellt werden. Hieraus sollen die zwei 

Konstanten A und ( der Gl. (1) ermittelt werden. A ist eine Spannung (Kraft auf der Flächen- 


einheit, Einheit z. B. dyn/cm? oder kp/em? °)), © eine reziproke Zeit (Einheit z. B. 1/sec oder 2. 
5) „Kilopond‘, technische Krafteinheit, früher Kilogramm genannt. BIN SR sm Dr 
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1/Std., wenn Q in cm?/sece oder cm?/Std. gemessen wird). Beide Versuche mögen mit demselben 
Rohr vorgenommen werden oder aber auch mit zwei verschiedenen Rohren vom Halbmesser R, 
und R,. Dieser allgemeinere Fall bringt keine Erschwerung der Aufgabe und mag deshalb hier 
als vorliegend angenommen werden. an 

Die in den theoretischen Formeln in den Vordergrund gestellte Dimensionslose a = p'R/2A 
ist hier zunächst nicht greifbar, da sie die Unbekannte A enthält; dagegen ist p/ R, und p,R; 
bekannt. Da a, und a, gesucht werden, empfiehlt es 
sich aus den Meßwerten auch Dimensionslose zu bilden. *9 
Vandrey wählte p/R,/p, Rz = € und ny/n, =, wo- 
bei n, und n, die rechnerischen Zähigkeiten nach 08 
Formel (11) sind. Wenn p,R, > p/R, ist, dann ist 
0<&E<1lundO <{<I1T, wie sich aus dem Vergleich 
von (9) und (11) ergibt. Wegen der Konstanz von A 
ist a, = &a,. Wir wollen a, als unsere Unbekannte 
ansehen; nach (11) ist, dan, = A/C ein Festwert ist, 


orlepla e. (13). 


Diese Gleichung kann gelöst werden, indem man an ı 
Hand von Tafel 1 die rechte Seite von (11) probeweise . 97 
für verschiedene Zahlwerte von a, ermittelt und sich 

so systematisch dem Zahlwert & annähert. Herr 
Vandrey hat für eine Reihe von runden Werten o 


von a und von & aus Gl. (13) die zugehörigen Werte 0 02 07 08 2 Na 
von ermittelt (die später noch etwas ergänzt wurden) > r 2 
und er hat hiermit ein Nomogramm hergestellt, mit Be De nenn ie 


dem sehr rasch Näherungswerte für a, gewonnen 
werden können, vgl. Bild 7. Das obige Probierverfahren braucht dann erst angewendet werden, 
wenn genauere Werte verlangt werden„als das Nomogramm sie liefert. Die Zahlentafel 2, 


Tafel2. Z abhängig von a und &®) 


a N) EN Be ERTER VERA SE N TAL 08 EN: 
a=1 0,898 0,900 0,915 0,933 0,945 0,960 | 0,980 1,000 
2 0,660 0,674 0,708 0,772 0,815 0,867 0,927 1,000 
3 0,417 0,433 0,487 0,586 0,658 0,748 0,859 1,000 
4 0,236 0,253 0,310 0,424 .| 0,512 0,630 0,788 1,000 
5 0,124 0,138 0,188 0,298 0,392 0,526 0,720 1,000 
. 6 0,062 0,072 0,111 0,207 0,296 0,436 0,655 1,000 
7 0,030 0,037 0,064 0,142 0,223 0,360 0,641 1,000 
8 0,014 0,018 0,037 0,097 0,167 0,298 0,541 1,000 
9 0,006 0,014 0,021 0,066 (0,125) 0,245 (0,492) 1,000 
10 0,003 0,004 0,012 0,045 0,094 0,201 0,445 1,000 
12 0,000 0,001 0,004 0,021 0,052 0,136 (0,363) 1,000 
14 0,000 0,000 0,001 0,009 _ | 0,029 (0,095) (0,306) 1,000 


die £ abhängig von «a, und £ wiedergibt, mag zur Herstellung des Nomogramms in größerem 
Maßstab dienen. 
Sobald a, ermittelt ist, findet man 


A=pR,/2a, . 


C = QylaR; ' Y(a,) . (15) 


mit y gemäß Gl. (9) und Bild 3. Für kleine Werte von a, wird das Verfahren ungenau, da hier 
wegen Annäherung an den Newtonschen Zustand A und C ihren Sinn verlieren. Sie nehmen 
bei « = 0 beide die Form 0/0 an und nur A/C liefert einen gut definierten Grenzwert. Für große 


. (14) 
und / 


‘a, dagegen, wo das Nomogramm versagt, kann aus der hier recht genauen Formel (8) für @(a) 


Q=nCR?:o(a) ad zum OAR?]/p, ; er'Rti2A . (16) 
Nach Bildung von Q,/Q, kann nach A aufgelöst werden. Es wird 
f k 5 2 ” } ; i ARE Q p. R: E 
deu R, ur) Fan, 7... see (17). 
2 Pi, — pi)! AmR: (17) 


Er Bed, Eingeklammerte Werte sind interpoliert. 


A 


% 
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Mit A ergibt sich aus (16)-auch ©. 
Für mittelgroße a kann jetzt aus der genaueren Gl. (8), gegebenenfalls unter mehrmaliger 
Iteration, ein verbesserter Wert von A gefunden werden. 
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Zur Bestimmung der Eigenwerte einer Matrix durch Iteration 
Von Alfred Vogel in Stuttgart 


Es werden verschiedene Verfahren zur Bestimmung der Eigenwerte ung Eigenvektoren einer Matrix 
angegeben und auf ihre praktische Brauchbarkeit für numerische Berechnungen verglichen. 


For the determination of the eigen-values and Ihe characteristie vectors of a matrix, several methods are 
given, and their usefulness for numerical computation is compared. 


Plusieurs methodes sont donne&es pour la determination des valeurs propres et des vecleurs propres d’une 
malrice; les methodes sont compar£es a l’Egard de leur utilit&E pour une Evaluation nume£rique. 

lamrcea paskmynsie MeTONBI ONpenesteHnst COÖCTBEHHBIX BHAyYeHHÜH MH BEKTOPOB MATPHILBI 
MH CPaBHUBAIOTCH OTHOCHTENIBHO HX HPAHTHYeCKOH IIPUTOAHOCTU IIPU BBIYHCJICHUAKX. 


1. Einleitung 


Nach der teils graphischen, teils numerischen Bestimmung der Eigenwerte gewisser Rand- 
wertprobleme durch L. Vianello!), J. J. Koch?) und F. H. van den Dungen?) haben 
R. von Mises und H. Pollaczek-Geiringer?) einen zusammenfassenden Bericht über prak- 
tische Verfahren zur Eigenwertberechnung gegeben. Sie legen das lineare Gleichungssystem 
r = AXr mit reeller symmetrischer Matrix WA zugrunde und bestimmen zunächst den kleinsten 
Eigenwert durch Iteration eines beliebigen Anfangsvektors ;„und Quotientenbildung entsprechen- 
der Komponenten von aufeinanderfolgenden Iterationsvektoren. Zur Berechnung des nächst- 
höheren Eigenwerts gehen sie von zwei Anfangsvektoren 3, und 3, aus, die hinreichend oft iteriert 
und dann linear kombiniert werden. Mit welchen Schwierigkeiten bei der numerischen Durch- 
führung man zu rechnen hat, ersieht man aus dem von den Verfassern vorgeführten einfachen 
Beispiel einer dreireihigen Matrix: Um den zweiten Eigenwert auf drei Dezimalen genau zu 
bestimmen, müssen sie zweireihige Determinanten aus zwanzigstelligen Zahlen bilden, von denen 
sie jeweils 16 richtige Stellen berücksichtigen. Seit der Arbeit von von Mises und Pollaczek- 
Geiringer sind wiederholt Untersuchungen über das gleiche Problem durchgeführt worden. 
In England haben W. J. Duncan und A. R. Collar?) Schwingungsprobleme mit Hilfe von 
Matrizen gelöst. Ihr lineares Gleichungssystem ist vom Typ Axt = Ar. Sie bestimmen zunächst 
durch Iteration eines beliebigen Anfangsvektors 3 den absolut größten Eigenwert A, (Maximal- 
wurzel) und den zugehörigen Eigenvektor r,. Mit Hilfe von r, bilden sie aus Neine neue Matrix X, 
deren Maximalwurzel mit dem absolut zweitgrößten Eigenwert }, von X übereinstimmt. Später 
untersucht A. C. Aitken®) den Fall einer beliebigen, nicht mehr notwendig symmetrischen, 
reellen Matrix W, insbesondere die Fälle von mehrfachen reellen und komplexen Eigenwerten 
mit linearen und nichtlinearen Elementarteilern. Er gibt Formeln zur Berechnung, von Aug 
e Bis Vianello, Graphische Untersuchung der Knickfestigkeit gerader Stäbe. Z. d. VDI, Bd. 42 (1898), 
D«. I» N 

®) J. J. Koch, Bestimmung höherer kritischer Drehzahlen schnellaufender Wellen. Verh. d. 2. int. gi 


Kongr. f. techn. Mech. Zürich (1926), 8.213. — Eenige toepassingen van de leer der eigenfunkties op 
vraagstukken mit de toegepaste mechanika. Diss. Delft (1929). 
») F.H.van den Dungen, Les Problemes gönsraux de la technique des-vibrations. M&m. d. sciences R 
phys. Bd.4 (1928). — Über die Biegeschwingungen einer Welle. Z. angew. Math. Mech. Bd. 8(1928),8.25. 
‘)R. von Mises und H. Pollaczek-Geiringer, Praktische Verfahren der Gleichungsauf- 
lösung II. Z. angew. Math. Mech. Bd. 9 (1929), S. 152. ; ; s SEA, 


5) W.J. Duncan und A.R. Collar', A method for the solution of oscillation problems by matrices. Bar 
Phil. Mag. Ser. 7, Bd. 17 (1934), 8.865. — R.A.Frazer, W.J. Duncan, A.R. Collar,,Elemen- 
tary Matrices. Cambridge (1938), 8. 143. | KR a ET 

)A.C. Aitken, Studies in practical mathematies II. Proc. Roy. Soc. Edinb, Bd. 57 (1937), 8. 269. . re: 
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A, A, A, usw. an. In der letzten Zeit haben sich N. J. Lehmann?) und W. Schmeidler') mit 
der Berechnung der Eigenwerte und Eigenvektoren durch Iteration befaßt. Lehmann schließt 
sich an die Untersuchung von von Mises und Pollaczek-Geiringer über die Bestimmung 
des absolut kleinsten und zweitkleinsten Eigenwertes von reellen, symmetrischen Matrizen an. 
Sein Verfahren für die höheren Eigenwerte ist eine Ausnutzung des Ritzschen Gedankens. 
Schmeidler bestimmt für Hermitesche Matrizen den größten Eigenwert und den zugehörigen 
Eigenvektor durch Iteration eines beliebigen Anfangsvektors und den nächsten Eigenwert 
ähnlich wie Duncan und Collar durch Abspalten; außerdem gibt er eine Fehlerabschätzung 
für den größten Eigenwert an. 

Es gibt Schwingungsprobleme, die sich unmittelbar in der Form Ar = Ar darbieten, 
z.B. die Berechnung der Eigenfrequenzen von Torsionsschwingungen. Die Iterationsverfahren 
liefern hierfür nur die größten Eigenwerte. Praktisch interessieren jedoch nur die kleinsten 
Eigenwerte. Durch Reziprokbildung der Matrix X könnte man auf das System r = AW-!r über- 
gehen und hierfür die kleinsten Eigenwerte bestimmen. Die Bildung einer Reziprokmatrix hoher 
Ordnung ist aber immer mit umständlicher Rechenarbeit verbunden. Es wird im folgenden 
unter anderem gezeigt, wie man auch für Systeme des Typs Ar = Ar durch zweimalige Berech-, 
nung einer Maximalwurzel zur Minimalwuürzel gelangen kann, ohne daß man reziproke Matrizen 
bilden muß. 


2. Bereehnung der Maximalwurzel und des zugehörigen \Eigenvektors 
Es sei Xeine beliebige reelle Matrix n-ter Ordnung. Ihre Eigenwerte genügen der charakteri- 
stischen Gleichung |AE — A| = 0. Zu jedem einfachen Eigenwert 4; gibt es einen linear unab- 
hängigen Eigenvektor rt, (vom Spaltentypus) als Lösung des linearen Gleichungssystems 
Ar; = A;r;, und einen linear unabhängigen Eigenvektor 9; (vom Zeilentypus) als Lösung des 
linearen Gleichungssystems 1; WM = A;yi. Im Falle einer symmetrischen Matrix A ist 1; — y.. 
Die Eigenvektoren r; und y; kann man aus der adjungierten Matrix A* (A) von AA) =, E— A 
erhalten. Diese Adjungierte A*(A;) ist eine singuläre Matrix vom Rang 1 mit proportionalen 
Reihen und läßt sich stets als Produkt A*(A,) = ©; 8; aus einer Spalten- und einer Zeilenmatrix 
darstellen. Es ist 1; = t6&; und y; = r&i. Die Proportionalitätsfaktoren it und können noch 
so bestimmt werden, daß r; und 9; der Normierungsvorschrift yit; =1 genügen, Gehören 
die Eigenvektoren r; und 9; zu verschiedenen Eigenwerten, so gilt die Orthogonalitätsrelation 
Ye =0 (ÜZ£ k). Hat WA lauter einfache Eigenwerte, so gibt es n linear unabhängige Eigen- 
vektoren rt; bzw. y;. Jeder Vektor 3 mit » Komponenten läßt sich nach den linear unabhängigen 
Eigenvektoren von A entwickeln. Ist z.B. 5 ein Vektor vom Spaltentypus, so läßt sich 5 in 
der Form 
x BE A Th RER Fa 8 1) 
mit skalaren Komponenten a; darstellen, die sich unmittelbar aus y; 3 = a,4ir; = a; ergeben. 
Man nennt den Eigenwert mit größtem absolutem Betrag die Maximalwurzel von %; 
sie braucht nicht unbedingt die algebraisch größte charakteristische Wurzel von X zu sein. Hat 
die reelle Matrix X nur positive Elemente, so hat sie einen reellen, positiven, einfachen Eigen- 
wert A, von absolut größerem Betrag als irgendein anderer Eigenwert von AP). Hier ist also die 
Maximalwurzel selbst positiv. Wir wollen von unserer Matrix X annehmen, daß ihre Maximal- 
ji wurzel einfach sei und daß es auch keine zwei verschiedene Maximalwurzeln mit gleichem 
absolutem Betrag gebe. Es sei also 
en, 
Zunächst soll die Maximalwurzel A, berechnet werden. Zu diesem Zweck multipliziere 
* _ man Wrechts mit einem beliebigen reellen‘ Vektor (Spaltenmatrix) 5 # vo, der nur nicht zu W3 
proportional sein soll. Dann bilde man nacheinander die Folge der Vektoren W3, (Az) = W’y, *-, 
allgemein AA” 37) = W"+1;,---, indem man X jedesmal rechts mit einem neuen Vektor multi- 
; pliziert. Dividiert man dann zwei entsprechend gelegene Elemente von W"*1 3; und W”3, so 
=. konvergiert für m—- oo die Folge dieser Quotienten gegen die Maximalwurzel },. Denn aus (1) 
ns folgt 
En Umz—aAry, tazältg hr" + Ann In: 
Für irgendein Element des Vektors A”3, etwa für die ö-te Komponente, hat man dann 
2 7“ lm) = aa Fir 4. 0, mt a0 0.2.0.6.) 
ET  ERENE ST, Lehmann ‚ Eingrenzung der Eigenwerte und Berechnung der Eigenvektoren homogener 
linearer Gleichungssysteme durch Iteration. Dresden (1947). Nicht veröffentlicht. Manuskript beim Lehr- 


stuhl für angewandte Mathematik der Technischen Hochschule Dresden. ! 
- 8) W.Schmeidler, Vorträge über Determinanten und Matrizen. Berlin (1949), 8. 99. 


ken 


__90.Perron, Zur Theorie der Matrizes. Math. Ann. Bd. 64 (1907), 8.248. und G. Frobenius, 


a RN Über Matrizen aus positiven Elementen, Sitz.-Ber. d. preuß. Akad, d. Wiss. (1908), 8. 471. 


ei 
F 
- 


176 


und als Quotient entsprechender Elemente von A”+1; und A”z entsteht so, falls!%) c, #0, 


film +1) Cy A,\m Ay\m fAsiem. K. > Er 
BD, (m) = Fum) = A, + % (A, —A,) © +0 e) + (0) e a . (4). 
Aus (4) folgt wegen (2) lim ®,(m) = A.. 


Bei der Auswahl der Vektorkomponenten hat man noch Freiheit. Man kann auch eine 
Mittelbildung vornehmen oder innere Produkte aufeinanderfolgender Vektoren bilden. Ist 4 
symmetrisch, so entsteht durch Multiplikation der Zeilenmatrix (A?3)’ mit der Spalten- 
matrix X”; (inneres Produkt des Vektors Y?3 mit A”3z), der Skalar i 


UA Wr; — Ss (pP +m) 
insbesondere für p= m ’ D 
| (Am ;) an; = UP sam) 2... u 222 10) a 
Die Quotienten R a 
s, (2m +2) I Se 
——— ., SE 
sam+ı) B 
konvergieren gegen A,. Mit Verwendung der inneren Produkte aus ‚gleichen oder benachbarten 
Vektoren der Vektorfolge W”3 ergeben sich also bei einer symmetrischen Matrix X nach m 
Iterationsschritten Skalare s, (2 m — 1) bzw.s; (2 m), während sonst nach m Schritten nur Skalare 
der Ordnung m zur Verfügung stehen, z. B. (3). Man braucht bei der Aufstellung der Vektor- 


folge für eine symmetrische Matrix X nur halb so weit zu gehen wie bei einer unsymmetrischen 
' Matrix, um die gleiche Güte der Näherung zu erzielen. 


2,am1) = 


@.. zii. na 


Wie man aus (4) erkennt, wird die Folge ®,(m) um so rascher gegen , konvergieren, je 


> du 


ei. 


ve | unter 1 bleibt!!). Man kann auch in den Fällen, wo das Verhältnis , ungünstig aus- 


fällt, die Konvergenz gegen A, dadurch beschleunigen, daß man aus der Folge ® ‚(m ) die neue 


Folge rk 
DS, $ | ® er 
Y,(m) = Fam) mit D&®,(m) = TE 1). et. | 
und RER, Kum. 
A®®, (m) : —= dB, (m +1) —2 B,(m) + DB, (m —1) 2% : 2 
bildet. Wegen (4) ist re | 3 n ” 


IN ).\m .\? m 
nm = +0) +0) 
\ \ 1 al 


N . A j 
: Der Fehler, den man begeht, wenn BER von einer Nummer m, an 9, ( m) : == A, setzt, ist also nur nee 


- 


“ nz n; 


u 


‚noch von der Ordnung Aa)" bzw. | 


v 


h, 
Fehler bei Benützung der Folge ®,(m) von. ‚der Ordnung i | 


h 

"je nachdem, welcher Wert der Brähere ist, während der Erol. 4 
1% ‘ 

N 

' 


Die Verbesserung Raun Er 


beträchtlich sein. Aus der Folge (6) kann a in gleicher Weise eine BASED Pegen E8 konver- Pr = 
gierende Folge \ du 


aq 
» 
S 
6 
| 
hs 


ER bilden. | | \ ee Re S REM Es 
| Im Falleiner positiv definiten symmetrischen Matrix Xist die Folg el 
» steigend, ihre Elemente sind untere Schranken fun größten Eigenwert A,. j 


trische Matrix A nicht positiv definit, so sind die Quotienten Kur au ur 


für A}. Nach oben gilt die An Le Ee u 
2 < s,(2m 19), 


F Fin) 


sı (2m) en a 
IRRE ER 
2%) Wenn zufällig = ,=...=a-1=0, c%7=0 ist, kon ars 


4%) In dem Beispiel vonSchmeidler, 2.2. 0,, das er dem Beric] 


Geiringer entnommen hat, ist dieses Verhältnis” gleich 0,06 


wöhnlich günstig. Re Ch 
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 Eigenwert 
FG. 


% 


N a ae Er a he Fi ve A an 
> rn, N um i 


nr 


mit | 32m +2) _ (2m) 
| 5,2 m — 2) a 
F(m) = en und ee. RL) 


un ae ET RENT HE BEER 3 1 
ks,(2m —2) — s, (2 m) 
Also hat man die Näherung®®) 


a ea N TO) 


mit einem Fehler, der höchstens + > - F(m) ist. 


Will man zum Eigenwert A, auch noch den zugehörigen Eigenvektor x, bestimmen, so kann 
man dazu die Vektoren der Folge X” ; verwenden, die gegen r, konvergieren. Für die praktische 
Rechnung wirkt sich jedoch nachteilig aus, daß man bei fortwährendem Iterieren schließlich zu 
recht großen Komponenten der Vektoren W”3 gelangt. Für die Rechnung bequemer ist es, 
wenn man nach jedem Schritt die Komponenten ähnlich verkleinert, also die Länge der Vektoren 
ändert, indem man sie auf eine Standardform bringt, z.B. so, daß stets die 1. Komponente 
gleich 1 wird, Dann ergibt sich folgender Rechnungsgang: 


pP, Ur 
Az = = (Pi); re: 
s=Pı Piı) 31 en 
p Ur 3 
AYı=r =); I a 8 
Pı2 Mia HL RE STD), 
” ® m+1 
Ugm = Int = (Pi,m-+1) 3 öm-1 =, = ee s 
; Aa, Pu 


Die Folge der Skalare 9,,, konvergiert gegen die Maximalwurzel A, und die Folge der Vektoren zm 
gegen den zugehörigen Eigenvektor x, mit der 1. Komponente 1. Für die praktische Durch- 
führung wird das so aussehen: wenn sich zwei aufeinanderfolgende Vektoren 3, und 3m-+1 Inner- 


- halb der durch die Rechenmaschine erreichbaren Genauigkeit nicht mehr unterscheiden, ist die 


Iteration abzubrechen. Man hat dann Azm = Pı,m+1 Zmpı mit praktisch gleichen Vektoren 
3m und 3m+1, so daß P,,m+ı für den Eigenwert A, und 3m = äm+ı für den zugehörigen Eigen- 
vektor y, praktisch erreichbare Näherungen darstellen. 


3. Bereehnung der Minimalwurzel 

Wenn man eine einzelne Matrix rationalen Operationen unterwirft, so machen ihre Eigen- 
werte entsprechende Operationen mit. So hat z. B. die reziproke Matrix A-! die Eigenwerte A7!, 
At, An, wenn A, Ay, **‘, An die Eigenwerte der nichtsingulären Matrix A sind. Für 
mehrere Matrizen hat G. Frobenius ") einen entsprechenden sehr allgemeinen Satz auf- 
gestellt: Es seien &, 9, 8, ... vertauschbare Matrizen n-ter Ordnung. f(z, y, 2,...) sei eine 
rationale skalare Funktion. Dann bilde man die Matrix R—= f(%,9, 8,...). Bezeichnet man 
die Eigenwerte von & mit &;, die von Q) mit %;, die von 3 mit 2, ..., während die Eigenwerte 
von X mit o(@=1,...,n) bezeichnet werden mögen, so kann man die Eigenwerte x, Yi, 2 » 


so ordnen, daß sich die Eigenwerte 0; von R durch ,; = f(, y, 2...) =1,...,n) dar- 


stellen lassen. 

Dieser Satz erlaubt eine für unsere Zwecke nützliche Anwendung. Die Matrix X habe 
lauter reelle Eigenwerte &,,..., &n. Mit einem von ihnen, z.B. mit x;, bilde man die Skalar- 
matrix 9) = %;€, die mit jeder Matrix gleicher Ordnung vertauschbar ist. Nach der rationalen 
Vorschrift f(®, y) = y— x stelle man die neue Matrix R= — = x,&— X her. Dann sind 
= %p-:, On = %— %n die Eigenwerte von RW. Sind A, >A, >... > 4n die Eigen- 
werte von W, so sind 0, = 0, 09, = A, — Ay... ., On = A, — An die Eigenwerte von # = A, E— N. 


R ‚Dabei ist on > On-ı >... >0,=0. A, ist die Maximalwurzel von X, on ist die Maximal- 
_  wurzel von ®. Mit Hilfe dieser beiden Maximalwurzeln erhält man die Minimalwurzel 
von X zu in=A,— On. Interessieren also z. B. die beiden kleinsten Eigenwerte von X, so 


RR) Für den größten Eigenwert des Eigenwertproblems Ar = Ar bei Schmeidler,a. &. 0. 8.103 in 
nlehnung an eine Darstellung von L. Collatz über Integralgleichungen; für sämtliche Eigenwerte des 

roblems g=AWXr bei Lehmann,a.a.O,., 8.21. 

robenius, Über vertauschbare Matrizen. $.-B. preuß. Akad. Wiss. (1896) I, 8. 601. 
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berechne man zunächst den größten Eigenwert 4, von X und hernach die beiden größten Eigen- 
werte 0, und 02-, von R. Die gesuchten zwei kleinsten Eigenwerte von X sind dann 
A -h,— | 
® Renz. oe 


P\ An-ı 4A, — On-ı Ä . 


4. Berechnung der nächsthöheren Eigenwerte 
A sei eine Matrix mit lauter reellen und verschiedenen Eigenwerten |A,| > A| > >1An. 
Die Maximalwurzel A, sei bereits berechnet. Gesucht werde jetzt der nächsthöhere Eigenwert ,. 
Bildet man aus der Folge (3) die neue Folge 


A a A En ee 

so streben die Quotienten j i x 
ea ge at =) RT re kai 
D,(m) = m) Aıds SR IBREr (A, As —- Ar 5) 1, Ho 1, +0 i ; 


alls c, # Oist, für m— oo gegen IE Die Konvergenz gegen A, A, kann man noch beschleunigen, 
D®, (m) 


wenn man zur Folge i 
9, (m) = — — mit‘. Dd, (m) = | en) 2 vn 
AFP,(m) | D, (m) B,(m —1) 
und 
A2?d,(mM)=d,(m +#+1)—29, (m) +9, (m —1) 


übergeht. Es ist 
; m—1 211—2 , 
gem to) Hol)” © 


Der Fehler, den man begeht, wenn man von einer Nummer m, an 95 (m) = A, A, setzt, ist also nur 


1, m—1 1.2 m—2 
Aa bzw. |® ‚je nachdem, Sa Wert der größere ist, w ährend. 


ki ser “ 
der Fehler bei Benützung der Folge ®, (m) von der Ordnung war. Trotzdem ist dieses Ver- 


noch von der Ordnung 


3 
A. 
fahren zur Berechnung des absolut zweitgrößten Eigenwertes nicht zu empfehlen. Bei der prak- 
tischen Durchführung der Rechnung treten Schwierigkeiten auf. Nach (13) ist 

fı(m+1) fı (m) | 


Yan) =Fulmg im Ur mn 


Nach (4) gehen sowohl A isn = & wie 7 „a n) mit wachsendem mgegen A,, so daß für Fe (m 
ı( Az: 
die Differenz von zwei wenig ver ARE, Zahlen zu bilden ist, was beisehr großen Zahlen steisz zu 


einem Verlust von wertvollen Ziffern führt. Größte Rechengenauigkeit ist erforderlich, um auch 
nur einige Stellen von A, zu sichern. Kleine Rechenmaschinen reichen dazu nicht aus. 
- Brauchbarer ist folgendes Verfahren: Die Maximalwurzel A, und der zugehörige Eigen- 
vektor g, seien nach (11) hinreichend genau bestimmt. Dann bilde man die neue Matrix 
A AT L (dm Ga ne Se en 
Mit Ausnahme von A, hat X, die gleichen Eigenwerte A,, A3, . . ., An wie die ursprüngliche Matrix X. 
In der 1. Zeile von %, stehen lauter Nullen. Deshalb ist Nullein Eigenwert von X,. Erist an die 
Stelle des Eigenwerts A, von A getreten. Der nach A, nächste Eigenwert A, von X kann deshalb 
auf gleiche Weise wie ” als Maximalwurzel der Matrix %, bestimmt werden. Man kann so As 
miterträglichem R öchenenewand sehr genau ermitteln. Nurist dazu erforderlich, daß r, genügend: 
genau bestimmt ist. Von der Güte von r, kann man sich dadurch SD daß man Ar, mit 
At, vergleicht. 
Wäre tr, genau bekannt, so könnte man den Freien Eigenwert A), und den zugehörigen 
Eigenvektor r, einfach wie folgt bestimmen: v sei ein beliebiger Anagst Seine Kompaz Be 


nente in r,-Richtung ist nen, t,. Der neue Vektor FA A Ne: 
1&1 . NE fi; £ Ei 
v’ 5 3 ai BEE 5% 
gern Fe BROEN 2 Se 
; Di Ear ea GR 
hat dann keine Komponente in r,-Richtung mehr, ist also zu r, orthognd, d. he Er 8- I 
wickelt man jetzt ® nach den een age von m so RR * FRE 


Sn Sri 
‘Es + 


„angew. Math. Mech. AN 


1 AR a a a a EL AR A rt 
Lu 


ERNG 
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und damit 
AR = a,an 2 + HF ande In 
Irgend eine Komponente von A” sei 
g2(m) = oz1y 4°, + 0ndn. 
Der Quotient entsprechender Komponenten von aufeinanderfolgenden Vektoren Ar+18 und 
RB wird dann 
9 (m +1) 
Jam) 
und strebt, falls ce, # 0, für m— oo gegen A,. Die Vektoren der Folge A” ® streben gegen 18 
Ist XY symmetrisch, so empfiehlt sich wieder die Bildung innerer Produkte 


(WS WB = VWI"® =s,(p+ m) 


X. (m) = 


insbesondere für p=m 
Ar]? = s, (2m). 
Die Quotienten 


Im Ba 
ee non yuaaemtd 


ss (amy N r (2m +1) 
konvergieren gegen Ay. e 

Nun ist aber r, leider nicht exakt bekannt, sondern nur näherungsweise, etwa in Gestalt 
des Vektors "3 von Nr. 2.,. Deshalb kann man den Ausgangsvektor ® nie so wählen, daß seine 
Komponente in r,-Richtung exakt verschwindet. Bei der Iteration X” ® wird sich die nicht ganz 
getilgte Restkomponente in t,-Richtung wieder auswirken. Das Kochsche Verfahren sieht des- 
halb vor, nach jedem Iterationsschritt erneut die noch verbliebene Komponente in r,-Richtung 
zu beseitigen. 

Die beste Näherung für r, in der Vektorfolge X” 3 werde jetzt der Kürze halber einfach mit 
31 bezeichnet. v, sei ein beliebiger Anfangsvektor. Wir nehmen ihm seine 3,-Komponente und 
schreiben 


\ voA 
1 
\ ! BD) =] > 4 S 
7 > 
ö1 091 


| . Der erste iterierte Vektor sei v1, = A 2,. Subtrahiert man von ihm seine 31-Komponente, so 


entsteht 
v13 
Bed a 1> 
11 
mit dem man nun wieder einen Iterationsschritt 0, = A ®, ausführen kann, usw. Allgemein wird 
Dim? 
Bam Au eo und — AB, ER ee (LO): 
5181 
E ; ar : r On? 
Im Verlauf der Rechnung wird sich nach einigen Schritten zeigen, daß der Anteil u 3ı keinen 
| 1ö1 
wesentlichen Einfluß mehr hat. Dann ist ®,, = d„. Von da an kann man iterieren, ohne nach 
E jedem Schritt die „Reinigung“ vornehmen zu müssen. Schließlich bilde man, falls X symmetrisch 
% ist, die inneren Produkte 
2 VB) Bun = t(2 m), R,, Vnti = t(2 Mm +1) a ABS HIRTEN (16). 
\ - Die Quotienten 
' {2 m-+1) {2 m +2) 
EL 2m +1) = — 17), 
T (2 m) 1m) T( +1) mr)! (17) 
'konvergieren gegen A,. Noch rascher konvergiert die Folge 
DT (m) 
ar 18 
| > m A®T' (m) (19) 
gegen Ay. 
% 5. Beispiel '*) 
3 "Die Bestimmung der Drehschwingungszahlen einer mit fünf Massen arbeitenden Kolbenmaschine führt 
auf nalgende Frequenzmatrix 
| 1 0,866025 0 0 
DE gr — [ 9866025 3° 1,095445 0 
BE", 0 1,095445 1 0,447214 
RN FE 0 0% 0,447214 2 


PAR der numerischen Durchführung des Beispiels hat mir Herr Dipl.-Ing. K.E. Hafner geholfen, 


Se] we ihm Sue: an Ger Stelle meinen herzlichen Dank aussprechen möchte. 
en 
N ee 3 
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1 
A ist symmetrisch und positiv definit. Mit dem Anfangsvektor 5 = ; erhält man beim 6.,7. und 8. Itera- Kl 
tionsschritt die Vektoren ‘ 1 zn? b DEN 
Us z DIKZF ; Us; £ | 3 a e 
1218,3373 4563,0911 17066,288 en 4 
3862,1908 14437,455 53936,195 - | 
1639,3086 ‚ 6090,7480 22675,292 I HD, 
493,32488 1719,7715 6163,411 
Nach (5) bilde man ?5) 

> s, (12) = |A° 5]? = 193315666 ; Sn — 13,931449 
s, (14) = |W 3]? — 26931673 (9); i in ' a RE 

5 116) = |U8 3]? = 37525278 (10); n a)” 13,933512. | | | a R 
Damit wird nach (9) und (10) ER; EN 
N) s| (14) | ’ N 
EN a NL ea 1. 
7, (1) m 5 0, (14) > — 13,933512 + 0,000646 | L Ba. 


ks,(12) —8,(14) 


x . w 
mit einem Fehler, der höchstens gleich 0,000646 ist. Also ist mit drei Stellen genau A? = 13,934 und A, = 3,733. 


Genauere Werte erhält man, wenn man die Folge (6) benützt und das Verfahren (8) zur Beschleunigung 
der Konvergenz anwendet. (Siehe Tafel I) _ 


K AR a fell 
m s,(m) | 2, (m) | An RE TFDE, | a o,(m) 3 
3 1 - 
Ri EURO 12 | 19331566 | » 3. 
U 97028873 i > 
M 13 | 72152612 - 225 . j 5 
® 598 (0,569) | 97 86806 895 | 3,732 895 - 
14 | 26931673 (?) 128 : 4 
726 (0,563) 56 49872 891 sg , Je 
15 | 10052856 (19) 72 72 RR | i 
16 | 37525278 (19) ’ rar 4, 
Also ist mit 5 Stellen genau A, — 3,73289. (Der genaue Wert ist 3 ‚1328913.) | \ | ie: 1 DR 7 
Das Verhältnis der ersten Differenzen J1, in Tafel T ist 0,569 bzw. 0,563, also nahen konstant; ei. 7 
konvergiert gegen 2: Mit dem mittleren Wert 0, 566 erhält man eine Näherung für A, und zwar x = 2, 11. , 
(Der genaue Wert ist 2,10690, also beträgt der Fehler der durch einfachste Rechnung erhaltenen Näherung® 4 2 
weniger als 0,15%.) Für die genauere Bestimmung von 2, reicht die Anzahl der Iterationsschritte noch nicht _ Bi 
aus. XP, W735, A®3 geben einen viel zu ungenauen Wert für den zu A, gehörigen Eigenvektor r,. Man müßte ER 
entweder die Anzahl der’ Iterationsschritte verdoppeln oder das Ko c h sche Verfahren PUSEBT Wir, RS ii 3 
_ führen diese 2. Methode vor. A A u 
Aus den drei Vektoren a | ERBE Dre Mn 
. er 2. SE 2) De NUDE. 
145 Az ME URS a ee 
4547,9203 4563,0911\ - ,./ 4571,8688\ er At be 
14417,137 14437,455 14448,907 5.) Wer Te SE 
| 6119,3603 6090,7480 ‚1 ..6074,4586 | VE Re 
u ’ 1841,5280 1719,7715/. 1651, 1092/ H q 2% # 
x T. bikle man durch Anwendung des Verfahrens zur Beschleunigung der; Krranan, wi Voktor nn 
E75 
4582,6281 | i 


14463,699 ER, 
. 6052.9260 oder 4582,6281 *3, mit 3, —g 


1562,3140/ N 


Er der als beste Näherung für x, genommen werden kann. = 


% 


Rechnung wechseln Base nach (15) ‚Beseitigung d. ER Ä d TI 
Tafel I.) 1m) Era a r 


15) Der obere Index (9) zeigt an, daß es sich um eine ahl it *2 
8 ersten angeschrieben wurden. - % ur 


er do z vd — AD, | 1P3 a AB, D5 =, AB; 
ar £ N .0,432070 0,023547 0,568038 
En 1 a —0,384192 —0,683040 —| 349013 
A ö 1 WISE 0,305419 —- 0,620652 1 "594295 : 
En = 1 So: 1.860821 3.850392 7.978329 
N ie el! , 
| WERT: v13ı 0341 did, 
Ian SH | He id" 
Be 0, 4537 32 0,020304 0,000015 —0,0000003 
Ar 1 ‚432070 y -0,064085 0,000047 —0,0000009 
0, "599308 - ... 0,026819 0,000020 —0,0000004 
$ 0; 154687 0,006922 0,000005 —0,0000001 
A 7 Do3ı N Dj 31 2 0331 i D33 
> Bud z V=-— oo: %,=-u,—- 7 2,0, 
Nez 3 a EEE Fa RE ae ERTL 
..0,546268 0,411766 0,023532 —0,568038 
7% N —0,432070 Sr ‘448277 —0,683087 —1,349012 
Be > = 0,400692 0,278600 .0,620632 1,594295 
E ini, s - - 0,845313 i 1.862899 3,850387 7,978329 
= | | = 
B. Yon 3. Schritt an ist YO = ms die 3,- -Kömponente‘ von by fällt nicht mehr ins Gewicht. Man kann jetzt 
nacheinander iterieren. Die nächsten drei Iterationsschritte liefern & E 
j Seh 
Br er, B, %, 
Be. 3 R —1,736318 —4,154706  . —-9,216642 
Be  . ee —12,792515 —5,845023 —12,281516 
> 2% e : 3,684545 8,080399 17,324209 Y% Er 
F n ENT : ; - 16, 669649 34,987078 73 ‚587824 % 


Nam kann man entweder eine Komponente herausgreifen und die Quotienten dieser Komponenten in auf- 
5% einanderfolgenden Vektoren bilden oder nach (16) "die inneren Produkte der % und ihre Quotienten (17) 
Bee herstellen und nach (18) diese RuospEnt eh gleich verbessern. (Sale Tafel IIT.) 


4 3 4 \ 2 l 


: : 2 Tafel IM. 
2 N H N hg e A 3 nl DT 2 
Sg | 17 2 
x >: ee ker} | al (m) | A X A Re DT 17 (m) 
ee er oz EL“ 
a a a: et N | 2,10 5864 
92636; a =1.r.8575 5 
2 AR bee 6439 | (0,44) | 318 | 2195377 | 6904 | 2,10 6904 
DE Pr ee TO, 1340,814 | u 2522 32 
FE ‚6696 (0,45) | 144 | 993265 | 6898 - 6898 < 
Para e 4 je 11 5 
we 11 B 2824,€88 6x0) SE . 
rn Du RT ‚5951, ‚078 e Ex Ö ER 


Stellen genau 2, = 8, 10690. (Der genaue \ Wert ist 2,106900.) Die Verhältnisse der 1. Differenzen 
ch ungefähr bei 0,445. Sie streben gegen = Damit erhält man einen Näherungswert A; = 0,937. 


hy 
ist 0, 93609.) s 
en wir uns noch für die RE kleinsten Eigenwerte A, und A, von Ninteressieren. Dies 
ßten Eigenwerte go, und Bs der Matrix R: =, &—X. Wenn wir mit dem letzten in Tafel I 
EV 132891 Teolnen, wird { 


- 0,866025, Te ) 


ER 0,732891 ee > 095445 ) 
en 1095445 en 0,447214 


o 447214 1,732891 


g Eeliorigei. Zigenvektor tr, bestimmen und benützen 


13 


« . 


lösenkann. Inder Tat muß jedoch solch ein schie- 
fer Stoß eine Änderung der Rotation zur Folge haben, 
wie stets die Wirklichkeit zeigt und wie im folgenden 
auch mathematisch nachgewiesen wird. 

Aus dieser mathematischen Untersuchung folgt 


weiter, daß das Reflektionsaxiom, wonach z.B. b 


eim, 


schiefen Stoß elastischer Kugeln streng mathematisch 
physikalisch der Einfallwinkel gleich dem Ausfall- 
winkel sein soll, eine unbewiesene, falsche Behauptung 
ist. Abgesehen davon, daß jeder baliebige schiefe Stoß 
eine Rotation des Stoßkörpers in Erscheinung treten 
läßt, drängt z. B. ganz offensichtlich folgender einfache 
Versuch zur exakten Klärung der Zusammenhänge. 
Wählt man als Stoßkörper gemäß Bildl einen an 
langem Pendelfaden aufgehängten, mit erheblicher 
Einschnürung (R >2r) versehenen Vollzylinder und 
bildet die Stoßfläche als in die Einschnürung ragende 
S:hneide aus, kann sich abgesehen vom Auftreten der 
Rotation zeigen, daß der Reflektionswinkel B.in völli- 
gem Widerspruch mit dem Reflexionsaxiom größer 
ist als z/2. Demnach ist die Voraussetzung berechtigt, 
daß bei schiefem Stoß eines den Abmessungen nach 
. bekannten Körpers auf eine feste Stoßfläche bei ge- 
gebenem Einfallwinkel und gegebener Stoßgeschwin- 


7 


--. Ebene, 
-  Körperschwerpunkt bewegt, 


fläche steht und sie inihrer 
"Mittellinie schneidet. Sonst 
würden noch räumliche 
 Kreiselbewegungen ausge- 
löst und die drei een er 
 tenwürdennichtausreichen.. 
Zur Ermittlung der dei 
Z Unbekannten stehen 
völlig elastischem schie 


Für nur drei Unbekannte 


ist dabei vereinfachende 


Voraussetzung, daß die 
in der sich der 


senkrecht auf der Stoß- 


Bob Luz SVertOgnnzE 
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Tafel IV. 
ETF ' Vektorfolge m 
Vektorfolge 3m 
Pım Pam Pam Pım 
\ 
a 1 3,6 2,7 4,3 2,2 
1) 0,7 19 0,61 3,38 2,72 4,35 1,59: 
ı| 0,80 29 0,47 343 2,87 4,61 1,39 
1| 0,84 34 0,40 - s 
3,508773 | 3,143550 | 5,191896 | 1,307455 
1 0,895911 1,479690 0,372625 3,503773 3,143550 5,191896 |: 1,307455 - 
Beim letzten Schritt hat sich keine Komponente mehr verändert. Also ist 0, = 3,508773 der größte Eigenwert 5 
i | | | \ 
-von R undt, = ER der zugehörige Eigenvektor. Mit seiner Hilfe bildet man nach (14) die neue Matrix 3, 
’ 4 “ 
0,372625 j 
N=R— 1, (2,732891  0,86605 0 0) 3 
0 N 0 0 4 
—1,582403 —0,042991 1,095445 0 a 
— 1 —4,043831 —0,186004 2,732891  0,447214 a 
—1,018343 05322702 0,447214 1,732891 
Jetzt wird genau so der größte Eigenwert von N, und der zugehörige Eigenverktor bestimmt. Wählt man % 
als Anfangsvektor 5 = ‚so erhält man als Maximalwurzel 03 = 2,796801 und als Eigenvektor 
RES Er 
: 2,592363 
0,786380/ " 
Nach (12) hat man mit o, und o, zugleich die beiden kleinsten Eigenwerte von X, nämlich > - 
)) = 1 — 04 = 0,224118 E 
Az > )ı 037 0,936090 . 
Eingegangen am 19. 9. 1949. ;. 
Untersuchungen des für den schiefen Stoß ela- . digkeit entgegen bisher Fi 
stischer Kugeln gültigen Reflektionsaxioms | üblicher Vorstellungsweise i 
und einige Folgerungen daraus. | mindestensdrei Unbekannte 4 
Im S:hrifttum findet man den vorsichtigen Hinweis, zu ermitteln sind: Der Aus- 3 
daß der schiefe Stoß in der Regel neben der fallwinkel, die translato- R 
translatorischen Reflektionsbewegung nach Stoßbeen- tische und die Botetions- y 
* digung auch Rotationen der stoßenden Körper aus- geschwindigkeit auf ihm. N 
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fangsgeschwindigkeit auf dem Trägheitskreis des Stoß- 
körpers; der zu jedem Glied der Energiegleichung ge- 
hörende Faktor m/2 fällt heraus; m ist die Masse des 
Stoßkörpers. Neben Einfallwinkel x und Reflektions- 
winkel ß werden zweckmäßig ihre Komplementwinkel 
A und B eingeführt. In Bild 2 ist 0 der Berührungs- 
punkt an der Stoßstelle, $ die Lage des Stoßkörper- 
schwerpunktes im Augenblick der Stoßbeendigung. 
In weiterer Vereinfachung sei der Stoßkörper eine ho- 
mogene Kugel. Strecke 0 — $ entspricht also ihrem 
Radius r. 


UCOSsA 


« 


Der Drall des Stoßkörpers für 0 als Bezugspunkt 
kann sich durch den Stoß nicht ändern, da die durch 
0 gehende Raaktionskraft als einzige äußere Kraft des 
Systems bezogen auf 0 kein Drehmoment hervorruft. 
Demnach gilt die Drallgleichung 


mrusin A= morsin B-+ mwre. 


"Hierin bedeutet cr den Trägheitsradius des Stoßkör- 


pers. Gekürzt ergibt sich: 
2. Drallgleichung: 


usinA=vsinB-+ we. 


Beim, Zentralstoß hat die Stoßzahl k den Wert 1 bei 
5 völlig elastischem, 0 bei völlig unelastischem und 
e 0< k<1 bei beliebig elastischem, Stoß. Bei schiefem 

Stoß entspricht die Normalkomponente des Stoßes 

dem Zentralstoß mit der Stoßzahl%k, die Tangential- 

-  komponente dem Tangentialstoß mit der hiermit neu 

- eingeführten, gleichfalls empirisch feststellbaren Stoß- 

=> zahl g, die noch na:h Einfallwinkelgrößen abgestuft 
BE werden kann. Gleitende Reibung tritt beim schiefen 
Stoß elastischer Körper im allgemeinen kaum in Er- 
scheinung, wie selbst bei ganz schiefen Stößen das 
Billardspiel lehrt und wie auch durch die plötzliche, das 

- Gleiten hemmende Verformung der Stoßflächen durch- 

- aus einleuchtend ist. Eine um so größere Rolle spielt 
> aber beim schiefen Stoß die Haftreibung, durch die 
"der Stoßkörper an der Stoßstelle plötzlich eine ruck- 
artige Abwälzbewegung ausführen muß. Erfahrungs- 
‚gemäß ist der Haftreibungsverlust bei solcher rollenden 
Bewegung sehr klein. Der Tangentialstoß wird durch 
die bei der Verformung auftretende Schubspannung 
in ähnlicher Weise federnd abgefangen wie der Zentral- 
stoß durch die Druckspannung. Soweit beide Span- 
7% nungsarten dem Hook schen Gesetz folgen, ist daher 
in erster Annäherung die Annahme berechtigt, daß 


“ 
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k=q. Fürk=qg=1, also im Idealfall völlig elasti- 
schen, gleitfreien schiefen Stoßes, haben sich bei Stoß- 
beendigung die beiden Geschwindigkeitskomponenten 
an der Stoßstelle genau umgekehrt und ihre Ursprungs- 
größen wieder erreicht. Das gleiche gilt für ihre Resul- 
tierende |wl. Führt man nach Föppl, ‚Dynamik‘, 
R. Oldenburg 1942, S. 246ff. die reduzierten Massen m} 
und m, ein und verlegt m, in O und m, in (, so folgt 
für $ als Kugelschwerpunkt sofort: 


2 
en Et ec? 1 
08=r7 80 =re; m =m 


Auf alle Fälle hat S genau wie m, und m, nach Stoß- 
beendisung die Normalgeschwindigkeit u cos A. Aus 
Bild 2kann man daher zunächst ohne Einschaltung der 
reduzierten Massen m, und m, unmittelbar ablesen: 


3. Geschwindigkeitsgleichung: 
uwcos AFvcosB. 
Aus 1 bis 3 folgt: 


\ 


26 ‚ 
W=U - sin A 
1.+ 


und nach Einführung von 


;5 11.77 
Ir 
BD. v=U V eos?A + n?sin?A 
6. tg B=ntgA, ferner 
n. sin Kr 


sinß u 


\ 
Diese Gleichung erinnert im Aufbau an das Brechungs- 


gesetz. Durch 6. ist der Nachweis erbracht, daß für 
den völlig elastischen Stoß das Beflexionsaxiom falsch 
ist; denn Körper mit ce = 0 gibt es nicht. Nach Ein- 
führung der Stoßzahl k = g bleibt 2. gültig, nicht aber 
1. und 3. Nach Stoßbeendigung hat jetzt die Tangen- 
tialgeschwindigkeitskomponente von m; den Wert uk 
sin A und die Normalkomponente von m} und: m, so- 
wohl wie von S den Wört uk cos A. Nach Bild 2 folgt 
dann geometrisch bei Einführung von 


1— ke? 1— ke? 
en Beer undz=r TS 
8. ukcos A=vcos B, 
9. tgB=ptgA, 
10. =Uu EL sin 4A, 
‚uk, v=uk ) co8® A + psim A 


Bei der Stoßzahl für Stahl (%k = 5/9) wird bei homo- 
genen Kugeln (c?—= 2/5) durch ein seltsames Zahlen- 
spiel der Natur p—=1. Die gleiche mittlere Stoßzahl 
dürfte für Eilfenbeinbälle auf Billardbande gelten. 
Diesem, den Augenschein täuschenden Zusammenhang 
verdankt das Reflexionsaxiom in Sonderfällen seine 
Gültigkeit. Nach 10. ist bei völlig unelastischem Stoß 
die durch den Stoß aufgeprägte Botationsgeschwindig- 
keit halb so groß wie die bei völlig elastischem Stoß. 

Ergibt‘ sich empirisch ein Unterschied zwischen k 
und g, so ist inBild 2 lediglich uk sin A gegen ug sin A 
auszuwechseln, um die entsprechenden Gleichungen 
anzuschreiben, und es folgt für diesen verallgemeinerten 
Fall bei unverändert bleibender Drallgleichung: 


1— ge? 
12. fe An lag 
201.4): 
13. a 
14. v= uk Yco® A+ n?sin® A. 
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Hat der Stoßkörper schon auf dem Einfallstrahl Ro- 
tation, liegen die Verhältnisse gleichfalls einfach, so- 
fern die Rotationsachse senkrecht auf der Strahlebene 
steht. Die Strahlebene sei die Ebene, in der der Ein- 
fallstrahl liest und die senkrecht auf der Stoßfläche 
steht; d sei genau wie w die auf den Trägheitskreis be- 
zorene Umfangsgeschwindigkeit, jedoch vor dem Stoß. 
Durch die Vorzeichen von d und w regelt sich ihr Dreh- 
sinn. Als Drallgleichung folgt dann: 


15. usind-+de=vsinB-+ we. 


Nach Stoßbeendigung ergibt sich als Tangential- 
geschwindigkeit für m, der Wert u sin A — d/c und für 
m, der Wert u sin A + de. Nach gleichem Verfahren 
wie zuvor gilt dann z.B. für den völlig elastischen 
Stoß: - 


16. 


Der Index 0 kennzeichnet hierbei die entsprechenden 
Werte des rotationsfreien Stoßes gemäß 4. und 6. 
Ferner gilt 


17.tg B=nstg A+(2dc) | (wcos AX1+ cc), 
wodurch auch » gegeben ist. 
Als Energiegleichung gilt die Kontrollgleichung 


18. wı+rR—=vr-LW. 


Die zugehörigen Gleichungen nach Einführung von 
q und %k regeln sich entsprechend den früheren Ablei- 
tungen. 

Hat eine stoßende Kugel vor dem Stoß Eigen- 
rotation um eine beliebige Schwerpunktsachse, 
so ist dieräumliche Lage der Geschwindigkeits- 
vektoren zu ermitteln. Dies ist eindeutig durchführbar. 
Insbesondere sind dabei die Kreiselgesetze zu berück- 
sichtigen, ferner die Tatsache, daß auch die räum- 
liche Richtung des Dralles vor und nach dem Stoß 
unverändert bleibt. Der Reflektionsatrahl liegt dabei 
im allgemeinen nicht mehr in der Strahlebene., 

Stoßen zwei Kugeln zusammen, so gilt sinngemäß 
alles Vorausgesagte. Zur Vereinfachung des Rechnungs- 
ganges nimmt man dabei den gemeinsamen Schwer- 
punkt des Systems als ruhend an und wählt ihn zum 
‚Koordinatenursprung. Zusätzlich zu den bisherigen 
"Buchstaben treten für die zweite Kugel die großen 
Buchstaben mit gleicher Bedeutung. Nach dem Schwer- 
punktsatz gilt dann: 


19. 


w=w, — An.. 


MU= mu, 
20. MV =mı, 
21. ‚em=/M, 
22. gm =hM. - 


@und f sind die lotrechten Abstände der Einfallstrahlen 
vom gemeinsamen Schwerpunkt und g und A die der 
Roflektionsstrahlen. Nach dem Scehwerpunktsatz haben 
beide Körper die gleichgroßen Einfallwinkel A. Dies 
gilt auch für B. Der Gesamtdrall des Systems, bezogen 
auf den gemeinsamen Schwerpunkt, wird durch den 
Stoß nicht geändert: 


23. mue+MUf=mvg-+- mwre+MVh+-MWRC. 


Das gleiche, gilt beim völlig elastischen, gleitfreien 
Stoß für die mechanische Energie: 

24. MU? + mu? = MV?+ MW: -+ mv? + mu. 
Bezeichnet man den Abstand zwischen Berührungs- 
punkt und Systemschwerpunkt mit H, so gilt noch 
geometrisch: - 

25. 9=(r+H)sinB. 

Durch 19. bis 25. sind die 7 Unbekannten P, », W, w, 
B, g, h, zu ermitteln. Bei Einbeziehung von k, q, d 
und D sind die reduzierten Massen einzuführen. M} 
und m, sind in den Berührungspunkt zu verlegen; 
alles Weitere regelt sich dann wie vorher beschrieben. 


’ 


Mitteilungen 


fd 


das Sonnensystem manche rätselhafte Spuren, deren 


‘wie in der Natur Umwandlungen von Hyperbel- in 


Hiermit sind wesentliche Hilfsmittel zur Klärung- : 
der Frage geschaffen, wie sich der elastische, schiefe Zn 
Zusammenstoß von zwei Himmelskörpern regelt. Ein 
solcher Zusammenstoß muß ein gewisses Maß von 
Elastizität (Federkraft) besitzen, da völlig unelastische 
Körper in der Naturlehre unbekannt sind. Auch hängt 
diese Fragenicht gänzlichin.der Luft, da z.B. Hertz die 
Stoßzeit zweier im Zentralstoß aufeinandertreffender 
Körper von der Größe der Erde in seine Rechnungen 
einbezogen hat. Sofern überhaupt Himmelskörper zu- 
sammenstoßen, ist der schiefe Stoß die Regel, die 
Wahrscheinlichkeit des reinen Zentralstoßes aber un- 
endlich klein. Die Wahrscheinlichkeit des Zusammen- 
stoßens überhaupt ist bei zwei Himmelskörpern, zu- 
sammen als Zweikörperproblem betrachtet, lediglich 
von der Größe ihres unveränderlichen Systemdralles 
abhängig, soweit er die translatorischen Bewegungen E 
der Körper auf ihren Kegelschnittbahnen betrifft. 
Durch ihn wird entschieden, ob sie, sich auf Hyperbel- 
bahnen aus weitester Ferne nähernd, miteinander in 
Berührung treten. Sofern dies geschieht, verursacht < 
der Stoß ruckartig eine mehr oder weniger radikale 
Änderung von Form und Lage der beiden momentanen E 
Kegelschnittbahnen, jedoch bei unverändert bleiben- PER 
dem Gesamtwert des Systemdralls. Einigermaßen ‘2: 
klare Voraussagen über den Verlauf eines solchen Vor- 
ganges sind denkbar, sofern insbesondere über die - 
Stoßzahlen k und g brauchbare Aussagen vorliegen. 
Je größer die Massen der Körper und damit ihrer Ent- 
weichungsgeschwindigkeiten sind, um so hypothetischer 
sind aber solche Aussagen; denn k und g sind von der 
Stoßgeschwindigkeit abhängig. Es handelt sich also 
dabei um Werte der Himmelsmechanik, die allenfalls 
der Rechnung, kaum .aber experimenteller Ermitt- 
lung zugänglich sein dürften. Gleichwohl sind z.B. 
bei hinlänglich klein und als Stahlkugeln gedachten 
Himmelskörpern unter voller Einbeziehung. der Schwer- 
kraft solche Werte exakt bestimmbar. Dann ist aber 
auch die Möglichkeit gegeben, den Verlauf der Kegel- 
schnittbahnen und Rotationen genau anzuschreiben, 
um unter Berücksichtigung durchaus möglicher Stoß- 
serien, vielleicht auch unter theoretischer Zuhilfe- 
nahme von Gezeitenwirkungen nachzuprüfen, ob und _ 


Ellipsenbahnen unter gleichzeitiger Anderung der 
Eigenrotationen und Minderung der mechanischen Ge- 
samtenergie durch den schiefen Stoß vor sich gehen 
könnten. Für den schiefen Stoß auf flüssige _ 
Oberflächen geben zudem z.B. die Untersuchungen 
„Über den Ricochetschuß“ (C. Ramsauer, Dis. 
Kiel 1903) gewisse Anhaltspunkte. Im übrigen zeigt 


Zusammenhang mitsschiefen elastischen Stößendenk- 
bar ist. - z Ss a - a 


q=—1 werden, sodaß da 


fürk=+1 nach 12. & 
wird. ° DE 
Rheydt. _ ey 


ne 


Ermittlung einer Reihendars 
In r in elliptischen Koc 
- a e 


E 


Für speziell Proble 
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23. - ern ellipt  Koordi- Schreiben wir 
ursen Ansatz einer emreihe. 


In 


a Bene 


n .n 

[02 „ o ‚ 

= [a I m —— ne (, 
(e —0,)(0 — 05) 


x 


ee Ä Er ER“ ) so liefert f,, keinen Beitrag für das Residuum im Null- 
in der Form dargestellt: Ba 
urn VERF En 


‚ punkt. Die Entwicklung von f_„ dagegen liefert 


ER: en ro). ES" S ei S" 3 
Bet = Realln (£ ae =i-=i — Real pe “gl j 5 u 0,05 2 PURE e 
_ Führen wir elliptische Bosinaten ein N die Schar “ © © utvon 
der konfokalen Bllipsen und Hyperbeln 2 B> | ex nn 
“ u gu +Lo’+1 
= -Cofw SET =u-iv; 0<v<2n u=u v0 1 & 
3 2 E mo SiW; Ww U Be V); 0sV/< Be und man findet nach Umordnung als Koeffizienten 
I Er s e 2 3), 1 
Shefarr ze) | 9  vont 
5 Se DEE {23 
= ST EI U = 0 R—1 1 
mit den umkehrbar eindeutigen Beziehungen 2 ZERERE 
— — Gofucosu; » = — bj U -cosV ae R 
7=+ Einusinv; y=-+ Sin UsinV » und wegen 0,0, = 1 
so ist wegen der Periode 27 von » folgender Ansatz | ; 
nahegelegt: De a" 0, PX 2) 
E ; 1 1 U ee : e 
Re u:  umw- Sw 2 Das Residuum im Nullpunkt besitzt demnach den Wert 
= 7 co Ugk: (5), ) i 77 
x 72 (a „cosnv + b„sin nv) Bes. f_„(=0)= m (8). 
1 1 
mit Ar \ z Als Beitrag im Punkte o, ergibt sich 
} ) 2 n 
& 2. f cosn v’ dv. EDER TERN E 
| nen) GW- | In(e= 91) 01— 0,’ : 
ER A=0,12... ni Bu = 
} e o=0,) = —— 
an  sinnv’ do’ en .) 0103 
se W — Gojw’ Demnach wird 
= .n -nW 
’ EN NEN 
= gutin —osiodor= Be A leere 2&SinW 
= so Sin die Fourierkoeffizienten durch folgende Inte- . so daß wir erhalten 
„Er Bet 4% N 123 & 
{ E > Lahr -nW 
a IR: 8, Soinw) . . . (10 
der © wri“ SE $ "und nach Integration 2 
SEE Fat dz SinW aW i 
ZI ER de. —» ee Han ra 
ze en: je | = 
u I : ER BE, Sen { 
SEE = Mr ET: ni SB ei W -+ const. — 2 = Cojnw. (11). 


san? nw 13 ur © ww. x 6 I "Betrachten wir nun den Fall, daß der Aufpunkt z 

= ’ innerhalb der durch den Punkt £ erklärten Ellipse liegt, 

EEE EEE BO liegen wegen U< u o, wie auch o, außerhalb des 

=. Integrationskreises ‚el =e", so daß nur die Anteile 

R 3 des Nullpanktsresi uums in Betracht kommen. Es 
) en sich mit 


. Rn SinnW 
n=0; Dr 9 = un 
2 er 2, — 0 Sin W 
Entwicklung 
lern 2 nw 


= < am“ SinnW 7.2) 
GE. u 


—&fn uw EE- 118): 
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Durch Zerlegung in Real- und Imaginärteil ergibt EAN ”: eererbind findet man fü ür die Intogrodifferetit. | 


für U >u DEE E 
In | -21= In VE Zee Em N? | | re 
is: : nlrwane-n non : 
-U- w2 a3, = 
RE 1, e(N), 020 
(Con ucosnvcosnV + EinnusinnvsinnV) (14) die Lösungen RER 
er. Are (di —2) = D. - a u us> u 3 
| 3 I REND ze rg F ; 
=2 1 2% ig 90) ! Er 2 U va a 
Ai 2 90) =cosnvV er —2n N Cojnu 3 R ES 
«(Sof nu cos nv sinn V — SinnusinnvcosnTV) p,()=sinnn = —2n er! Be : 
und für U<u U 2, 
In|&—2] =nY(—- 2? +m-y% a } 
0.24) 


kl 


mad 1, =2ne"" SinnU 
= ne ""Cojn U 


.(Cof n U cosnV cosnv+ Sinn U sin nVsin nv ar 
as e (15), Weitere Integralformeln ergeben sich z.B. auch da- 


% , Arc —2)= durch, daß man die vorangehenden nach den Para- 
S Se) metern U,u oder V differentiiert. 82 Er ne man aus 
“ R en (17) 
Ir m n - ge: On a ER 
; . 22 S 4 2 
ei ; ? TE 5 n . 
FR -(Cojn U cosnV cosnv — Sinn sinnV sinnv), duouU i elolar (U, u; V,») dv Be: 
&) ® v i x 2 
DE wobei die Integrationskonstanten mit U— x und A (U,u) VER x ee 
durch Vergleich bestimmt worden sind. = FU p(P). BE Pa Te g 
Aus diesen Reihenentwicklungen folgen einige spe- $ N BE ER N 
zielle Bilinearformeln, von denen die für In Is —ı| Mtu=U=0 und Er a ee a A 
die bekannteste ist!): ; 2: Ne 
In ]E— z| = In |cosv — cos Be“ u hör — Sin noosnz; + Coluninnz, Per e: 
o erhält man aus 9 und (19). die en Integral 2 Ass EN; x 
| = —_1n2—2 > cosnncosnV , (16). gleichungen | BSN; Er 
} n Bl Na BER, ne SE E- 


; €) ? N 4 i 
; IV. Anwendungen. ee e|dE = urn 3 b 

' Mit Hilfe der gegebenen Reihenentwicklungen sind TArre BEL: ee 
wir in der Lage folgende a | zu be- = cos 4 


handeln 


a er 


2n j 
S yWw)lar(U,u; V,v)dv = 
f) 


= DE 
2 -  AU,u)p(V), U20 EB 


Nach (14) und (15) wird sie gelöst durch die folgenden 
Eigenfunktionen und Eigenwerte 


U>u \ 
YM)=l A = 2r(U—1n2) 
R c gnU vu 
5 n=e—%r Coinu 


BE de B grU Ar 
Mm) =sinnv, A, =—?%ı Sinnu 
et (17a). 
Er u 
Betr 3 Er ; 


4 $ ar 
An SU -2n- 


Ein nu Ra 6753 


Wien man Be ar RR des Deebnt 
\ = u ch in ‚die ‚des Problems U a ee DER 


' Z.angew. Math, Mech. 
Bd. 30 Nr. 5/6 Mai/Juni 1950 


wegen der Vollständigkeitsrelation von der Form 


5 =A(U-m2—-2. Er 


-(A,CojnucosnV +B, SinnusinnV) | 
Der Sonderfall einer Streckenbelegung auf der x-Achse 

„zwischen — 1 und -+ 1, für den ein cos-Reihenansatz 
genügt, liefert mit 


o(v) = N A, cosnv 0 <v S2n 
) 


7) Kr = grUT 
>= = As(U-m2) == - 


Die übrigen erwähnten Kernfunktionen mit ihren Ent- 
wicklungen erlauben es, auch andere potentialtheo- 
retische Probleme zu behandeln, wie z.B. die Frage 


A, cos nV. 


nach der Wirkung einer auf ‘Ellipse oder Strecke ge- -» 


legenen Wirbel- oder Quellschicht. 


Hannover. K. Jaeckel. 


Experimentelle Untersuchung zur Grenz- 
schichtentstehung an Zylindern bei Anfahrt 
aus der Ruhe. 

Die Blasiusschen Ansätze [1]!) zur Berechnung 
instationärer Grenzschichtströmungen aus dem Jahre 
1908 — noch bis heute der wesentliche Bestand unserer 
theoretischen Methoden zur Behandlung instationärer 

\ Grenzschichtprobleme von einiger Allgemeinheit — 

beziehen sich auf die Ausbildung der laminaren Grenz- 
schicht an Zylindern bei Anfahrt aus der Ruhe relativ 
zur umgebenden Flüssigkeit. Diese Ansätze sind in- 
zwischen einerseits analytisch und numerisch weiter 
ausgebaut und atch für rotationssymmetrische Pro- 
bleme durchgeführt, zum anderen auf allgemeinere Be- 
wegungsgesetze der Anfahrt aus der Ruhe übertragen 
worden [2, 3, 5, 6, 7,8, 9, 10]. (Der Falleiner Strömung 
an einer unendlich ausgedehnten ebenen Wand, der,. 
wie einige andere spezielle Strömungen auf eine lineare 
Differentialgleichung vom Wärmeleitungstyp führt, ist 
kürzlich unter neuartigen Bedingungen untersucht 
worden [11].) 

Da meines Wissens noch keine quantitative Nach- 
prüfung der Blasiusschen Theorie. vorliegt, er- 
scheint es mir gerechtfertigt, im folgenden einige ex- 
perimentelle Ergebnisse aus dem Jahre 1939 mit- 

R zuteilen, deren Veröffentlichung bisher unterblieb. Die 
etwas mühsame Untersuchung wurde mit den mir da- 
mals gerade zur Verfügung stehenden Mitteln am Kai- 
ser-Wilhelm-Institut für Strömungsforschung in Göt- 
tingen durchgeführt. 

“Man erhält nach Blasius in verhältnismäßig ein- 
facher Rechnung eine befriedigende Approximation des 
Grenzsshichtverlaufs vom Zeitpunkt des Bewegungs- 
beginns an für nicht zu große Fahrzeiten. Da die Theo- 

 rie von der Voraussetzung Gebrauch macht, daß die 
Potentialströmung außerhalb der körpernahen Grenz- 
schicht durch diese nicht wesentlich gegenübar der rei- 

 bungslosen Strömung modifiziert wird, kanm ohnehin 

- jedenfalls nıch einer stärkeren. Grenzschichtablösung 

_ eine brauchbare Approximation der wirklichen Ver- 

“  hältnisse nicht mehr erwartet werden. (Vgl. zu dieser 

Ag Frage [8].) <E 


‘ 


Die im folgenden geschilderten Beobachtungen be- 
ziehen sich auf den ruckartigen Übergang 
eines elliptischen Zylinders aus der Ruhe in die gerad- 
zu uftune Bewegung in Richtung einer Haupt- 

se des 


Zylinder). Die zugehörigen Ergebnisse der 
det man in meiner Arbeit [10]. Die Experi- 


{ kieen, Klammsra angegebenen Zitfern verweisen 
ı Sc) 


BER DEE ERBE 3 


ir - 
5 u 5 
EM FR 
ne ir Re E: 
FE 
Er 1 > 
Vera 
en 
. 


Kleine Mitteilungen 


gungen (Bodenerschütterungen, 


- Ellipsenquerschnitts (,„‚nicht-angestellter‘‘ 


n in einem- für Sshleppversuche eingerich- sungsstelle mußte darauf geachtet werden, daß solche 


se dieser Mitteilung genapnte Literatur. 
1f a ze ter > r 
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teten Wassertank desgenannten Institutsvorgenommen, 
derin der Literatur mehrfach ausführlich beschrieben 
worden ist, vgl. etwa [4]. Sie betrafen Ort und Zeit des 
Ablösungsbeginns und das erste Stadium der Wande- 
rung der Ablösungsstelle längs der Zylinderkontur für 
weiter wachsende Zeiten. Die Sichtbarmachung der 
Strömungsvorgänge erfolgte durch Aufstreuen feiner 
Teilchen (Al-Pulver) auf die Wasseroberfläche (Methode 
von Ahlborn). Die Bewegung der Teilchen beim 
Hindurchschleppen des Zylinders durch die Flüssigkeit 
wurde von einem mit dem Zylinder mitbewegten Auf- 
nahmeapparat kinematographisch festgehalten. Die 
ausnutzbare Fläche des Tanks betrug 340 x 50 em? 

Beim Einschalten der elektrischen Antriebsvorrich- 
tung bewegte sich der Zylinder in guter Näherung nach 
dem theoretisch gut erforschten Gesetz des ruckartigen 
Übergangs aus der Ruhe in die geradlinig-gleichför- 
mige Bewegung. Die Beschleunigungszeit (Zeit vom 
Bewegungsbeginn bis zur praktisch vollkommen gleich- 
förmigen Bewegung) betrug in den angestellten Ver- 
suchen nach sorgfältiger Ermittlung 0,2 bis 0,3 sec. 

Um die zeitliche Entwicklung des Bewegungsvor- 
gangsin der Flüssigkeit leichter zu erfassen, mußte an- 
gestrebt werden, die Zeit i = t4 vom Bewegungsbeginn 
bis zum Einsetzen der Grenzschichtablösung möglichst 
groß werden zu lassen. (Selbstverständlich sollte sie 
zumindest groß sein gegenüber der oben genannten Be- 
schleunigungszeit.) Nach der Theorie ist (vgl. etwa 
[10], Gl. (20)) 

m? a 
L--E:m.U7: 
Dabei ist m = b/a das Verhältnis der Ellipsenhalbachse 
b quer zur Bewegungsrichtung zur Ellipsenhalbachse a 
in Bewegugsrichtung und U ist die Fahrgeschwindig- 
keit des Zylinders. Ferner ist m < 2/3 vorausgesetzt, 
wonach die Ellipse so schlank sein soll, daß nach der 
Theorie Ablösung im hinteren Staupunkt einsetzt. 
Nach dieser Formel waren somit kleine Fahrgeschwin- 
digkeiten U und große Zylinderquerschnitte zu wäh- 
len. Zur Vermeidung eines mierklichen Einflusses der 


ta = 0,64 


. Kanalbreite (50cm) durften jedoch die Zylinderab- 


messungen nicht zu groß gewählt werden, und die Fahr- 
geschwindigkeit durfte wiederum nicht zu klein sein, 
damit die Geschwindigkeit der zu verfolgenden Flitter- 
teilchen noch im Mittel hinreichend groß blieb gegenüber 
den Geschwindigkeiten der unvermeidlichen Störbewe- 
thermische Bewe- 
gungen nach Einschalten der starken Beleuchtungskör- 
per), die nach unseren Feststellungen unter 1 mm/sec 
lagen. Die Fahrgeschwindigkeiten betrugen in den 
schließlich angestellten Versuchen etwa 8 mm/sec. Die 
unten ausführlicher wiedergegebenen Resultate ins- 
besondere beziehen sich auf einen Zylinder mit a = 
10cm und b=5cm, also m =!/,. Nach der Theorie 
war somit eine Ablösungszeit t4 von etwa 1,3 sec zu er- 
warten. Die Versuche wurden über das Vierfache dieser 
Zeit ausgedehnt. Es wurden Filmaufnahmen gedreht 
mit etwa 25 Aufnahmen pro Sekunde. Die Bewegung 
eines Uhrzeigers und eines an einem raumfesten Lineal 
mit dem Körper entlanglaufenden Zeigers wurden zur 
sorgfältigen Ermittlung von jeweiliger Laufzeit und 
-geschwindigkeit mitphotographiert. 

Die Auswertung der Filme erfolgte durch Projektion 
auf einen Wandschirm mit aufgezeichneter Ellipse, 
wobei das jeweilsins Auge gefaßte Al-Teilchen von Bild 


zu Bild leicht zu verfolgen war. In Schritten von im 
" allgemeinen fünf Bildern wurden Lage und ‚Uhrzeit auf 


der Wandschablone eingetragen (also in Zeitschritten 
von etwa 0,2 sec, womit die jeweils in der anströmenden 
Flüssigkeit zurückgelegten Wege etwa 1,6 mm, auf dem 
Projektionsschirm bei der gewählten dreifachen Ver- 


 größerung somit 4,8 mm betrugen). Bei der Verfolgung 


der Al-Teilchen in der Nähe der momentanen Ablö- 


Teilchen für die Beobachtung ausgeschieden wurden, 
die beim Überholen eines anderen Teilchens einen Zu- 
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sammenstoß mit diesen erlitten oder sich gar mit diesem 
durch Hängenbleiben zusammenballten, ferner daß ein- 
zelne an der mit einer Paraffinschicht versehenen Zy- 
linderoberfläche hängenbleibenden Teilchen nicht zu 
merklichen Störungen des verfolgten Teilchens Anlaß 
gaben. S:hließlich mußte selbstverständlich für eine 
völlig saubere Wasseroberfläche gesorgt werden, um 
auch durch Verschmutzung verursachte Zusatzbewe- 
gungen auf der Oberfläche auszuschließen. Dies alles 
erschwerte sehr die Auswertung, bei der pro Film je- 
weils nur bis zu drei der verfolgten, hinreichend wand- 
nahen Teilchen in Nähe der wandernden Ablösungs- 
stelle als einwandfrei akzeptiert wurden. Bei den Ver- 
suchen erfreute ich mich der reichen Erfahrung und des 
besonderen Geschicks von Herrn W.Müller, beider 
Auswertung der großen Geduld von Frl. J. Wolters. 

Eine gewisse unvermeidbare Ungenauigkeit der Be- 
obachtungsmethode mußte hingenommen werden: Die 
Teilchen wandern zunächst mit unerheblich abnehmen- 
der Geschwindigkeit bis fast unmittelbar an die ihnen 
entgegenwandernde (stromaufwärtsrückende) Ablö- 
sungslinie, bremsen dann verhältnismäßig rasch ab, 
rücken dabei sehr langsam etwas von der Wand «ab, 
bleiben relativ lange Zeit praktisch fast in Ruhe, um 
dann in eine rückläufige Bewegung überzugehen; es 
sind daher der momentane Ort der Ablösung und die 


zugehörige Zeit nicht mit der wünschenswerten Ge- - 


nauigkeit zu ermitteln. Ferner scheint das Verhalten 
der Grenzschichten an der Wasseroberfläche etwas von 
dem Verhalten in größeren Tiefen abzuweichen. 

Die Beobachtungsergebnisse für den genannten Zy- 
linder mit m =!/, sind im nebenstehenden Bild ein- 
getragen und mit der Theorie (ausgezogene Kurve, 
„zweite Näherung‘“ im Sinne meiner Arbeit [10]) ver- 
glichen. Dabei ist die Abszisse u der zur jeweiligen Ab- 
lösungsstelle gehörige Inkreiszentriwinkel der Ellipse, 
gemessen vom hinteren Staupunkt aus (vgl. [10]), und 
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Zusammenhang zwischen Ort der Ablösungsstelle und Fahr- 
weg des Zylinders (m=0,5) nach Theorie (ausgezogene Kurve) 
und Experiment. 


die Ordinate stellt den in Vielfachen von a aufgetra- 
genen, vom Zylinder zurückgelegten Fahrweg $ dar. 
Wenn auch die Ausbeute angesichts des Aufwandes be- 
scheiden bleiben mußte, so sind die Ergebnisse doch 
aufschlußreich. Nicht nur die Zeit des Ablösungs- 
beginns wird von der Theorie gut wiedergegeben, son- 
dern — und dies ist bemerkenswert — die Wanderung 
der Ablösungsstelle während der weiteren Fahrt des 
Zylinders stimmt noch bis zu Zeiten von etwa 44 be- 
friedigend mit den Ergebnissen der Theorie überein, 
einer Theorie, die eine Modifikation der äußeren Po- 
tentialströmung durch die anwachsende Grenzschicht 


und insbesondere durch die Ablösung derselben vom 
Körper unberücksichtigt läßt. Die Veränderung der Tr 


“ 


äußeren Strömung infolge der Verdrängungswirkung 
der Grenzschicht beeinflußt offenbar die zeitliche Wan- 
derung der Ablösungsstelle erst sehr spät. ig 
Entsprechend sorgfältige Messungen wurden auch 
bei Querstellung des Zylinders (also m = 2) und für 
‘einen Kreiszylinder (m = 1) vorgenommen. Die Er- nr 
gebnisse sind vernichtet worden. Es kann daher hier 
nur mitgeteilt werden, daß auch in diesen Fällen Ort 
und Zeit des Ablösungsbeginns (für m = 2 an zwei zum Pt 
hinteren Staupunkt symmetrischen Stellen der Kon- ehe 
tur einsetzend) in gleich guter Übereinstimmung mit x 
der Theorie beobachtet wurden. Für weiter wachsende 
Zeiten ergaben sich bei diesen, größeren Dickenverhält- 
nissen m etwas früher systematische Abweichungen 
zwischen Theorie und Experiment für den Zusammen- 
hang zwischen Ort und Zeit der wandernden Ablösungs- 
stelle. ER 
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gerungen vernachlässigbar sind. Der ‚Regler folgt da- 
mit der Gleichung sa 


. r 
Inn 9+r, | Pdt=np+ 2 |pdı (2). 


Dabei bedeutet NEO, = 
r, = Verstärkungsverhältnis des Reglers bei ‚starr‘ 
gemachter Rückführung (d.h. 7, — »). 
T, = Abklingzeitkonstante der nachgiebigen Rück- 
führung. 2 - 
pr = vom Regler eingestellte Stellgröße. 


Die beiden Konstanten r, und 7; sind bequem aus der 
Übergangsfunktion des Reglers zu ermitteln, ebenso 
sind die Bestimmungsgrößen s, und 7’, (bzw. s,) der 
Regelstrecke aus deren Übergangsfunktion abzulesen. 

Sucht man nun die Grenzeinstellungen des Reglers 
auf, bei denen der Regelvorgang gerade instabil (an- 
schwellend) zu schwingen beginnt, dann können dieser 
Grenzeinstellung Dauerschwingungen konstanter Am- 
plitude @, zugeordnet werden: 9 =9@,sinwt. Für 
solche Schwingungen kann das aus der Elektrotechnik 
her bekannte Zeigerbild gezeichnet werden, in dem alle 
Glieder der Gleichungen (1) und. (2) als Zeiger (rotie- 
rend gedachte Vektoren) dargestellt sind. Man erhält 
so Bild 1, worin die beiden Zeiger 3 und ur infolge der 


Bild 1. Zeigerbild des Regelvorganges, 


Laufzeit 7} um den Winkel »T7 gegeneinander ver- 


dreht sind. Aus diesem Bild kann man folgende beiden 
Gleichungen ablesen: 


A a r 12 +2) - 2 +82 w* 
SE | ze °o i se rninkie 
123 alien (3b). 
Mit r_, = ro/Tr und ru/s, = V =Verstärkungs- 


z 


kreises folgt aus da) 0» 


verhältnis des aufgeschnitten gedachten Regel- 


f 


7 - > (oT)? = 3m — 1) = 


> #y tr -1P + rm | 
unter Berücksichtigung von %—t g oT, 
ro: En - Ne . En 


- 


; Regelkreises auf- 
und 7’, noch gerad 
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BEER SDORIX- Dabei ist r,, das Verstärkungsverhältnis des Reglers 
bei ausgeschalteter Rückführung. Da im allgemeinen 


(starre Rückführung) seinen Größtwert. Ändert man 
T , dann schreitet man auf den strichpunktierten Li- 
nien fort; auch mit wachsendem 7, wird das zulässige V 
größer und bei 7, = 0 (d.h. nur 77 und T},) bleibt ein 
endlicher Wert für V zwischen V=0 und V=1 
übrig. Andert man schließlich 77, dann bewegt man 
sich parallel zur Ordinate; bei 7; = 0 wäre immer Sta- 
bilität vorhanden. 


Neben der Teilung für V ist noch eine Teilung für den 


-Regelfaktor R=1/(1-+ V) angebracht. Daran 


kann man ablesen, auf welchen Betrag die Auswirkung 
einer Störung durch den Regler im ersten Augenblick 
verringert wird; nach einer längeren Zeit wird die Aus- 
wirkung der Störung infolge der Nachgiebigkeit der 
Rückführung ja sowieso beseitigt. 

Fir 


N 


Bild 2. Einstelldiagramm. 


In Gl. (2) ist angenommen, daß der Regler beim Aus- 
schalten der Rückführung astatisch wird und. somit 
dann die Gleichung ur = r-, [ p dt annimmt; gültig 
beispielswetse für nicht federbelastete Kolbenstell- 
motoren. In der Verfahrenstechnik werden jedoch 


häufig federbelastete Stellmotore, beispielsweise Mem- 


branventile, benutzt, die den Regler auch beiausgeschal- 


'teter Rückführung statisch machen. Statt Gl. (2) gilt 
- jetzt: 


el 2 ai lt (2a). 
EuRET m uRdt nP+7, pdt (2a) 


rg, >r, gewählt wird, um die beabsichtigte ausglei- 


chende Wirkung der nachgiebigen Rückführung gegen 
dauernde äußere Störungen des Gleichgewichtes zu er- 


eichen, wirdin Gl.(2a) das Glied mit [ur dt klein und 
(2a) kann dann mit genügender Genauigkeit durch 
rsetzt werden.- 
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Das Ausschlag-Zeit-Diagramm einer „ein- 
fachen Schwebung“. 2 


In dem (im übrigen ausgezeichneten) Werk von K. 
W. Wagner!finden sich über die sogenannte ‚‚ein- 
fache Schwebung‘‘ Ausführungen, die es angezeigt er- 
scheinen lassen, die Frage nocheinmal kurz zu streifen. 

Unter einer Sshwebung versteht man die Überlage- 
rung zweier harmonischer Schwingungen, 


x = A,sinozt+ Assin(wt+0)...(D) 


die benachbarte Frequenzen &, und &, aufweisen. ES 
ist bekannt, daß ein solcher Vorgang sich darstellen 
läßt als eine sinusförmig amplituden-modulierte Schwin- 
gung. Wenn die Amplituden der beiden Summanden 
gleich sind, A} = A, = A, spricht man von einer ein- 
fachen S:hwebung. Sie läßt sich darstellen als 


u) a re X) & 
AR, 1 2 Sin (Sr ee 
© 2 Acos| RER. ‘) sn 5) +3) @ 


Esist ferner bekannt, wie das Ausschlag-Zeit-Diagramm 
durch Projektion des ‚‚erzeugenden Diagrammvektors‘“ 
(bei Wagner: ‚‚Speers“) zustande kommt: Die er- 
zeugenden Vektoren a, und a, der beiden Teilschwin- 
gungen addieren sich; es gibt im Abstand 


TrT= 2rllwı = @5) 


Augenblicke [bei den Nullstellen des ersten trigono- 
metrischen Faktors in (2)], wo der erzeugende 
Vektor der Summenschwingung verschwindet. Soweit 
herrscht Übereinstimmung in allen Darstellungen. 


An die zuletzt beschriebene Erscheinung, das Ver- 
schwinden des resultierenden erzeugenden Vektors, 
schließt Wagner nun folgende Bemerkung ?) an: 


„Beim Durchgang des resultierenden Speeres durch 
Null ändert er zugleich seine Richtung um 180°; die 
Schwingung erleidet an dieser Stelle einen entsprechen- 
den Phasensprung, derin Abb.18 gleichfalls er- 
kennbar ist.‘ 


Bild 1. 


Die Abb. 18 ist hier als Abb. 1 wiedergegeben; sie 
zeigt bei den „Schwebungsknoten‘“ eine Unstetigkeit 
der Richtung der r—t-Linie, also eine Unstetigkeit der 
Ableitung (Geschwindigkeit) #. Eine solche Unstetig- 
keit kann es aber nicht geben. Man erkennt dies schon 


‘einfach aus der Darstellung des Vorgangs nach Gl. (2). 


Die Ableitung enthält wieder nur Produkte aus Sinus- 


“und Cosinus-Funktionen; Sprünge im Sinne von Un- 


stetigkeiten können nicht zustande kommen. 
Welchen Sinn kann man nun dem Ausdruck ‚‚Pha- 
sensprung‘‘ unterlegen ? " 
Die meisten der Nullstellen der Funktion x(t) werden 
durch die Nullstellen des zweiten harmonischen Fak- 


tors in (2) bezeichnet. Die Funktion z(t) verschwindet 


1) K. W. Wıgner, Die L>hre von den Schwingungen und Wellen, 
Wiesbaden 1947, Dieterichsche Verlagsbuchhandlung. 
®) a.a.0. S. 15, letzter Absatz. S 
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dort deshalb, weil der erzeugende Vektor (Speer) in die 
Projektionsrichtung fällt; es verschwindet also die Pro- 
jektion eines von Null verschiedenen Vektors (Speers). 
Wir wollen diese Nullstellen ‚reguläre Nullstellen‘ 
nennen. Daneben gibt es noch ‚‚überzählige‘‘ Null- 
stellen. Sie sind die Nullstellen des ersten (harmoni- 
schen) Faktors. Dort verschwindet x(t) deshalb, weil 
der erzeugende Vektor selbst verschwindet. 


Nun unterscheiden wir zwei Fälle je nachdem, ob 


die überzähligen Nullstellen mit einer der regulären 
Nullstellen zusammenfallen oder nicht. Wir beginnen 
mit der zweiten Möglichkeit: Bei den Mullsialen von 


in (4) 


wird die Zeitachse abwechselnd im steigenden und im 
fallenden Sinn durchlaufen. Sobald eine „überzählige‘“ 
Nullstelle auftritt, kehrt sich (weil der Faktor 


N ar 
cos ( 5 t = 


sein Vorzeichen ändert) diese Zuordnung um. Man kann 
diese Umkehrung der Zuordnung von Nummer der Null- 
stelle und Richtung des Durchgangs durch die Achse 
einen ‚„‚Phasensprung‘“ nennen. Der erzeugende Vektor 
der Summenschwingung hat seinen Phasenwinkel plötz- 
lich um 180° geändert. Das Bild der Funktion x(t) 
sieht dann so aus, wie Abb. 2 angibt. 


- Bild 2. 


“ 


5 . \ 
Fällt die überzählige Nullstelle mit einer der regu- 
-Jären Nullstellen zusammen, so sieht die Funktion so 
aus wie Abb. 3 angibt. Auch hier tritt ein „Phasen- 


„sprung“ beim erzeugenden Vektor auf. Er zeigt sich 


wieder in der Umkehrung der Zuordnung von Nummer 
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Dr.se. techn. F. Schultz-Grunow (o. Prof. d. allg. 
Mechanik a.d. Techn. Hochschule Aachen), Ein- 
führungindieFestigkeitslehre. 2448. 
m. 239 Abb. Düsseldorf-Lohausen 1949. Werner- 
Verlag G.m.b.H. Preis brosch. 15,— DM. - 

Das umfangreiche seitherige Schrifttum über Festig- 
keitslehre und Elastizitätstheorie kann man im wesent- 
lichen nach zwei Gesichtspunkten ordnen: Einerseits 
handelt es sich um Bücher, welehe nur die ganz ele- 
mentare Festigkeitslehre nebst den vereinfachten Ge- 
brauchsformeln vermitteln; andererseits um Werke 
über höhere Festigkeitslehre und Elastititätstheorie, 
welche weitgehende mathematische Kenntnisse voraus- 
setzen und daher ihren Leserkreis vorwiegend in Krei- 
sen der. Hochschulstudierenden und der akademisch 
gebildeten Ingenieure suchen. Die Kluft zwischen dem 
wissenschaftlichen Niveau beider Bucharten ist infolge 
der starken Entwicklung der theoretischen Festigkeits- 
forschung ständig im Wachsen begriffen, so daß sich 
bei den Studierenden bereits ein gewisser Mangel an 
solchen Büchern herausgestellt hat, welche den Leser 
mit einfachen und mehr anschaulichen Methoden bei 
einem Mindestaufwand an mathematischem Rüstzeug 
von der elementaren Festigkeitslehre aus zu einzelnen 
für die Praxis besonders wichtigen Teilgebieten der 
höheren Festigkeitslehre und Elastizitätstheorieführen. 

Neben den heute vielfach vergriffenen Föpplschen 
Büchern gehört hierher auch das in dieser Zeitschrift 
besprochene Büchlein von ©. Weber. Das nunmehr 
vorliegende Buch des Aachener Mechanikprofessors 
Schultz-Grunow ist gleichfalls der Aufgabe 
gewidmet, eine Einführung in die (elementare und 
x höhere) Festigkeitslehre zu geben und dabei mit den 

einfachsten mathematischen Methoden auszukommen. 
Entsprechend dem Buchumfang konnte der Verf. 
seinen Stoff auf wesentlich breiterer Grundlage ent- 

5 wickeln, als es ©. Weber in seinem kurz gefaßten 

Buche möglich war. Dank der vielen klar ausgeführten 
Abbildungen und Erläuterungen führt der Verf. den 

Leser verhältnismäßig gründlich in die Vorstellungs- 

welt der Festigkeitsbegriffe ein. Einzelne Probleme 
sind‘weitgehend durchgerechnet und diskutiert, wenn 

auch die Stoffauswahl des als Einführung gedachten 
Buches für den Gebrauch des in der Praxis stehenden 
Ingenieurs natürlich nieht ausreichen kann; dieser 

. wird sich jedoch nach dem Studium des Buches mit 
um so größerer Leichtigkeit den im Rahmen der je- 

3 weiligen speziellen Aufgabe erforderlichen Spezial- 

{® werken zuwenden können, Schließlich wird auch der 

| . angewandte Mathematiker manche Anregung vorfinden 

- obwohl ihn naturgemäß manche mit der Anschauung 
geführten Herleitungen und Beweise nicht befriedigen 

+ werden, deren exakte Formulierung er jedoch leicht in 


; = der Spezialliteratur nachsehen kann. So hilft das neue 
gr Buch, die Kluft zwischen der einfachen Festigkeits- 
2a lehre und der höheren Theorie zu schließen, und es 


. wird sicher von weiten Kreisen begrüßt werden, 
Eu Dresden,, H. Neuber, 


SE Helmut Hasse (o. Prof. a. d. Univ. Berlin, Direktor 
am Forschungsinstitut für Math. d. deutschen- Aka- 
- demie d. Wiss. zu Berlin), Zahlentheorie — 
 Oktav bröschiert, XII + 4688. Berlin 1949. Aka- 
-Verlag. Ladenpreis geb. 44,— DM. 

as Werk ist der erste Band einer geplanten Gesamt- 
stellung der modernen Zahlentheorie. Es stellt die 


ndung dar, wie sie von Hensel, ‚dessen 


7 
1 


Igebraischen Zahlen in divisorentheore- 


lage wieder. Der Verfasser selbst kennzeichnet den we- 
sentlichen Unterschied dieser Begründung der algebra- 
ischen Zahlentheorie gegenüber der idealtheoretischen 
folgendermaßen: ‚Während dort eine zunächst willkür- 
lich erscheinende Definition des Begriffes der ganzen 
algebraischen Zahl an die Spitze gestellt wird, ergibt sich 
hier der Integritätsbereich I der ganzen algebraischen 
Zahlen aus K, der der Arithmetik in K zugrunde liegt, 
zwangsläufig aus dem bewertungstheoretischen Ansatz, 
und seine Kennzeichnung in der bisher üblichen Weise, 
durch den Koeffizienten der algebraischen Gleichungen, 
erscheint als ein Satz, nicht als Definition‘ (S. 221). 

Hinsichtlich der großen Fülle des verarbeiteten Stof- 
fes verdient besondere Hervorhebung, die für das Stu- 
dium des Werkes glückliche Art, in der die Analogie der 
algebraischen Funktionenkörper mit derjenigen der al- 
gebraischen Zahlkörper dargestellt wird. Allen Ab- 
schnitten, in denen das wegen der anzubringenden Ab- 
änderungen notwendig ist, wird der arithmetischen 
Theorie der algebraischen Zahlkörper ein Anhang mit 
derjenigen der algebraischen 'Funktionenkörper einer 
Unbestimmten hinzugefügt. \ 

Eine abschließende Bemerkung, die den Bereich der 
Betrachtungen des mathematischen Inhalts überschrei- 
tet und eine Betrachtung über das Wesen dieses mathe- 
matischen Inhaltes ist, erseheint mir zur Gesamtcharak- 
terisierung zweckmäßig: das Werk als Ganzes ist dank 
seiner überlegenen Gedankenführung (z. B. bezüglich 
der Einordnung der älteren elementaren Zahlentheorie) 
ein klarer Ausdruck für die tiefgehenden Änderungen, 
die sich in der Aufgabenstellung und logischen Grund- 
legung der Zahlentheorie vollzogen haben. In der Be- 
zeichnung ‚‚neuere‘‘ Zahlentheorie steckt wesentlich 
mehr, als dieses Wort vermuten läßt, steckt eine Struk- 
turänderung, die sich als Übergang von Größen zu 
Formen kennzeichnen läßt und von der ein anderer 
Autor einmal sagt: „Das Reich der mathematischen 
Größen kann den mächtigen Zug zur Form nicht ver- 
bergen; es ist, als bräche seine Sehnsucht nach Form 
um so unbezwinglicher hervor, je mehr seine unlogische 
Seele sich verleugnen will.‘‘ Wenn ein rein mathema- 
tisches Werk, in dem von mathematischer Philosophie 
explizite nyjt keinem Worte die Rede ist, dem aufmerk- 
samen Leser auch diese weitesten Perspektiven er- 
schließt, so ist das ein weiterer Vorzug, der überdies ein- 
drucksvoller wirktalsin einem großen Teilder Literatur 
zur Philosophie der Mathematik. 


“ Dresden. M. Draeger. 


Dr.P. Luckey (Oberstudienrat i. R.), Nomo- 
graphie. (Mathematisch-physikalische- Bibliothek, 
Reihe I Nr. 59/60) 6. Aufl. 107 S. m. 57 Abb. Leipzig 
1949. Verlagsgesellschaft B. G. Teubner. Preis kart. 

; T. 
Das treffliche kleine Bändchen, das im wesentlichen 
seit der in Bd. 8 (1928), S. 497 dieser Zeitschrift be- 


- sprochenen zweiten Auflage unverändert blieb, er- 


scheint jetzt in sechster, vom Verfasser nicht durch- 
esehener Auflage. Auch so bleibt es eine gute Ein- 
ührung in die Nomographie, die sich vor allem zum 


- Dresden. Willers. 


-Dr.K.L. Schecher und Studienrat Dr. K. Seebach 
(Privatdozent a.d. Techn. Hochschule München), 
EinführungindieMathematik. Ein Lehr- 
buch für Selbstunterricht und Wiederholung. Bd. T. 
288 S. m. 408 Abb. Legelshurst (Baden) 1950. Ver- 


Selbststudium eignet. 


lag Dr. Schecher & Co. KG. Preis geb. 13,— DM. 
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und daher schlecht einprägsame äußere Aufmachung, 
Die wichtigsten Ergebnisse usw. hätten drucktechnisch 
besser hervorgehoben werden müssen. Die Darstellung 
wirkt dadurch eintönig und ermüdend. An .dem 
geometrischen Teil ist zu bedauern, daß er durchweg 
in dem heute veralteten Stile euklidischer Darstellung 
gehalten ist und manchen entbehrlichen Ballast ent- 
hält. E ; 

Einige Einzelheiten, die bei einer Neuauflage der 
Korrektur bedürfen, seien noch herausgegriffen: S. 13: 
daß ‚‚3a + 2b nicht weiter zusammenfaßbar ist‘, an der 
Nichtaddierbarkeit von 3 Hasen und 2 Apfeln zu er- 
klären, ist nicht angängig. Die Buchstaben haben die 
3edeutung von Zahlen und nicht von Benennungen. — 
Auf 8.24 wirdnichtklar, was a° = list: eine Definition, 
ein Grundsatz oder ein beweisbarer Satz. Gerade 
darauf aber kommt es für den Fortgang und für die 
Einsicht.in das Wesen der Mathematik sehr an. — 
5.29: Dieimaginären Zahlen als „‚die nicht wirklichen, 


die nur in der Vorstellung existierenden‘ einzuführen, - 


bedeutet die Verewigung längst überholter Irrtümer, — 
Auf $S. 178 £ werden die Winkelfunktionen als Strecken 
definiert ; sie sind aber reine Zahlen und daher die Maß- 
zahlen jener. Strecken. 


Dresden. M. Draeger. 


Hans Wehner, Untersuchung mikro- 
barographischerWellen aufJanMayen 
(Veröffentlichungen des Geophysikalischen Institutes 
der Univ. Leipzig. Herausgegeben von W. Hesse. 
2. Serie: Spezialarbeiten aus dem Geophysikalischen 
Institut und Observatorium, Bd. XVI, Heft). 56 8. 
mit 12 Abb. u. 9 Tabellen. Berlin 1949. Akademie- 
Verlag. Brosch. 7,50 DM. 

Verf. weist nach, daß die auf der Insel Jan Mayen 
registrierten kurzperiodischen Drucksehwankungen dy- 
namischen Ursprungs sind. Zu ihrer Erklärung be- 
nutzt er eine Karmansche Wirbelstraße und zeigt, daß 
die Beobachtungen die Stabilitätskriterien einer solchen 
erfüllen. Die Schwankungen sind Auswirkungen lo- 
kaler Verhältnisse, Die von dieser Station gemeldeten 
Druck- und Windwerte brauchen daher nicht immer 
repräsentativ für die Großwetterlage zu sein. 

Dresden. Willers. 


Dr.-Ing. Gerhard Grüß (o. Prof. an der Bergakademie 
Freiberg), Differential- und Integral- 
rechnung. (Mathematik und ihre Anwendung in 
Physik und Technik. Herausgegeben von E. Kamke 
und A. Kratzer, Reihe A, Bd.21). XV u. 6428, 
m. 355 Abb. Leipzig 1949. Akademische Verlags- 
gesellschaft. Preis geb. 39,—-DM. - 

Ein gutes Buch über Differential- und Integral- 


„ rechnung können wir in Anbetracht des immer stärker 


werdenden Interesses für den Gegenstand immer noch 
gebrauchen. Eine gewisse Konkurrenz wird-auch hier 
ganz nützlich sein, Das vorliegende Buch hält Ref. 
für durchaus konkurrenzfähig. Es bringt den Stoff, den 
jeder Studierende der Mathematik, der Physik und der 
Technik zu allererst kennen lernen muß, in ansprechen- 
der, ja eindringlicher und anschaulicher Vortragsweise,. 
355 Abb. erleichtern das Verständnis. Numerische und 
instrumentelle Methoden werden angemessen vor- 
getragen. Auch sie sollte ja jeder Anfänger kennen 
lernen, auch der Studierende an der Universität. Was 
die Strenge anbetrifft, so erscheint sie Ref. für eine 
Anfängervorlesung durchaus genügend. Das Buch 
wurzelt, wie Verf. selbst in der Vorrede sagt, in der Ära 
der Cauchy-Weierstraßschen Mathematik. In dieser 
Beziehung hätte Ref. einen Vorschlag zu machen. Bei 
der Einführung -der Exponentialfunktion wird diese 
zwar durchaus richtig erklärt, dann aber wegen des 
Beweises des Additionstheorems auf $ 34.1 verwiesen, 
wo das einfacher geschehen werde ($. 65). An der ge- 
nannten Stelle aber findet man zwar eine rein arith- 
metische Theorie der trigonometrischen Funktionen 
vom Integral aus, hinsichtlich der Exponentialfunktion 
aber nur auf Seite 388 eine kleingedr uckte Bemerkung, 
a ee =, Br 
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daß man den Logarithmus und also auch die Ex- 
ponentialfunktion vom Integralaus begründen könne, 
ohne daß dies wirklich geschieht. Hier Klafft zweifellos 
eine Lücke, die Ref. zu schließen vorschlägt. Ami 
besten an beiden Stellen; es ist ja auch gar nicht so 
schwer. ns 

Das Buch gliedert sich in neun Kapitel: 1. Die 
Grundlagen: Zahl und Funktion. ı2. Die elementaren 
Funktionen. 3. Grenzwert und Stetigkeit. 4. Diffe- 
rentialrechnung. 5. Integralreehnung. 6. Weiterer 
Ausbau der Lehre vom bestimmten Integral. 7. Taylor- 
sche Formel und Taylorsche Reihe. 8. Unendliche 
Reihen. 9. Geometrische Anwendungen der In- 
finitesimalrechnung. Es handelt sich also im wesent- 
lichen um Funktionen einer Variabeln; von komplexen 
Größen wird nur einmal, bei der Partialbruchzerlegung, 
ein kurzer Gebrauch gemacht. : 


Landshut. Hamel, $ 
Wilhelm Scholle, Die mittlere Entfer- E 

nungeines Punktes voneiner Fläche, ® 

(Thünen-Archiv der Universität Rostock: Beiträge zur _ 

Thünen-Forschung. Herausgeber Prof. Dr. E. Schle- 

singer u. Prof. Dr. A. Petersen. Heft2.) 798. mit F 

33 Abb. Berlin 1949. Akademie-Verlag. Brosch. u 

12,50 DM. g ö 
Das für Nicht-Mathematiker geschriebene Heft stellt 

die von Thünen auf bestimmte Dreiecke angewandte 3 

Methode der Auswertung des Doppelintegrales zur Be- 

rechnung der mittleren Weglänge bezogen auf einen 

Eckpunkt ausführlich dar und erweitert sie auf gerad- 

linig begrenzte Flächenstücke und Kreise auch bei 

anderen Lagen des Bezugspunktes. Die Darstellung 

ist vielfach ungeschiekt und nicht fehlerfrei. Nurkurz 

wird erwähnt, daß das Integral, das heute in der Sied- ä 

lungsplanung viel verwendet wird einen besonderen H- 

Namen ‚Vial‘ bekommen hat, und man zu seiner Aus- h 

wertung für beliebig begrenzte Flächenstücke und be- 

liebige Lage des Bezugspunktes in der Praxis Um- 

fahrungsplanimeter mit Gleitkurventrieb verwendet, 

deren Genauigkeit durchaus dem praktischen Be- 

dürfnis entspricht. 


“Dresden, . -Willers, “3 
Dr.-Ing. habil. Otto Luetkens. Die Methoden _ 
derRahmenstatik,Aufbau,Zusammen- 
fassungundKritik. VII+ 2818. mit38 Abb, 
und 9 Zahlentafeln» Berlin— Göttingen— Heidelberg 
1949. Springer-Verlag. Preis brosch. 33,— DM, geb. NE 
36,— DM. ans FR f - 
Die Literatur zur Statik des Stabwerks hat teils 
durch die zahlreichen technischen Probleme, die dr ey 
Stahlbau und der Stahlbetonbau in den Jahrzehnten Schul } 
einer stürmischen bautechnischen Entwicklung ge- 


den Hilfsmiti 
zwischen ‚den iı 


 geometrisch unbestimmten Komponenten des Ver- 
: schiebungszustandes- aus linearen Gleichungen und 
ihre Vereinfachung bei einfacher, mehrfacher und zy- 
2 klischer Symmetrie des Tragwerks durch Umordnung 
der äußeren Kräfte. Hieraus haben sich die zahlreichen 
Verfahren der Baustatik ergeben, die unter dem Namen 
; Ausgleichsverfahren, Abklingungsverfahren usw. be- 
” kannt sind und oft nur in einer Approximation der 
Lösung bestehen, die den technischen Zweck erfüllt. 
Diese Methoden werden erörtert, auf viel verwendete 
" Tragwerke angewendet, durch gut gewählte Zahlen- 
- beispiele erläutert und vor allem aufihre Brauchbar- 
2 keit geprüft. Darin liegen Bedeutung und Eigenart des 
Buches, mit dem sich der Verfasser den Dank vieler 
ge, Ingenieure erwerben wird, die für die Baustatik und 
nn ihre Zusammenhänge Interesse besitzen und neue 
Erkenntnisse für ihre beruflichen Arbeiten suchen, 


Be Dresden, K.Beyer. 


- Dr.-Ing. habil. Ernst Eckert, Einführungin 
j den Wärme- und Stoffaustausch. VII 
+ 203 S. mit 125 Abb. Berlin-Göttingen-Heidelberg 

1949. Springer-Verlag. Preis brosch. 21,— DM, geb. 

. 24,— DM. 

R Das vorliegende Buch behandelt die Wärmeüber- 
| tragung durch Leitung, Konvektion und Strahlung. 
Dem Stoffaustausch ist nur am Schluß ein kurzer 
Abschnitt über Dampf-Luftgemische und den Diffu- 
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sionsvorgang bei der Verdunstung gewidmet. Das 
Kennzeichnende der Behandlung des Wärmeaus- 
tausches ist die Anwendung der Prandtl’schen Grenz- 
schichttheorie auf den konvektiven Wärmeübergang 
in Verbindung mit der Kärmänschen Näherungs- 
lösung im Gegensatz zu der sonst in der Literatur 
bevorzugten AÄhnlichkeitsmethode. Die Benutzung 
der Grenzschichttheorie ist daher besonders dankens- 
wert, zumal sie in außerordentlich klarer und über- 
zeugender Weise dargestellt ist. Weniger gelungen 
ist dagegen die anschließend allzu summarisch be- 
handelte Ähnliehkeitstheorie. Den gleichen Einwand 
kann man gegen die Darstellung der Gasstrahlung 
innerhalb des Abschnittes Wärmestrahlung erheben, 
zumal gerade diese zu dem engeren Forschungsgebiete 
des Verfassers gehört. Es wäre vielleicht im. Interesse 
des sonst vorzüglichen Buches besser gewesen, zu- 
gunsten einer ausführlicheren Darstellung der tech- 
nisch so wichtigen Gasstrahlung auf den im Hinklick 
auf den Buchtitel enttäuschenden Abschnitt Stoff- 
ausgleich zu verzichten. Der Stoff wird belebt und 
vertieft durch eine große Anzahl gut ausgewählter 
Rechnungsbeispiele. Eine umfangreiche Tafelsamm- 
lung von Wärmeleitzahlen und Zähigkeitswerten 
beschließt das den Fachkollegen empfehlenswerte 
Buch. 


Dresden. Bas 1tin: 


Die besprochenen und angezeigten Bücher sind durch den Buchhandel zu beziehen. 


Dr. Ernst Hölder (Prof. a. d. Univ. Leipzig), Über 
dieVariationsprinzipederMechanik 
der Kontinua. (Berichte über die Verhand- 
lungen der Sächsischen Akademie der Wissenschaften 

zu Leipzig, Math. Naturwiss. Klasse, Band 97, Heft 2) 
138. Berlin 1950. Akademie-Verlag. Preis brosch. 
- 2,75 DM. . 


- Dr. phil, habil. H. Gebelein, Zahlund Wirk- 
lichkeit. Grundzüge,einer mathematischen Sta- 
!  tistik. Zweite Auflage. ‘XII + 430 S. mit 52 Abb. 

. Heidelberg 1949. Quelle & Meyer. Preis 12,— DM. 


Prof. Dr. Friedrich Tölke, Mechanikdefor- 
mierbarerKörper. Erster Band: Der Punkt- 
förmige Körper. VILL-+ 388 S. mit 339 Abb. Berlin- 
.  Göttingen-Heidelberg 1949. Springer-Verlag. Preis 
Be geb.245,-= DM, 


_ Otto Toeplitz, Die EntwieklungderIn- 
 finitesimalrechnungnachder gene- 
tischen Methode. (Die Grundlehren der 

mathematischen Wissenschaften i. Einzeldarst. Bd. 
 — LVL) Erster Band: Aus dem Nachlaß herausgegeben 

von G.Köthe. IX -+181S. mit 148 Abb. Berlin- 
 ,  Göttingen-Heidelberg 1949. Springer-Verlag. Preis 


‘ _brosch. 19,60 DM., geb. 22,60 DM. 
Dr. 6. Hamel (em. o. Prof. a. d. Techn. Univ. Berlin), 
oretische Mechanik. Eine einheitliche 
ung in die gesamte Mechanik. (Die Grund- 
r mathematischen Wissenschaft in, Einzel- 
Band LVIL) XVI-+ 7968. mit 161 
öttingen-Heidelberg 1949. Springer- 
brosch. 63,— DM., geb. 66,— DM. 


ıd (o. Prof. a. d. Univ. Jena), Ein- 
ie Theoretische Physik. 
eorie der Wärme, (Meyers. 
ıer, Band 54/55.) ‘330 S. m. 46 Abb. 
B Diographisches Institut VEB. Preis. 
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EINGEGANGENE BÜCHER 


_ Bei der Schriftleitung sind folgende Bücher ‘eingegangen (ausführliche Besprechung bleibt vorbehalten) 


6. Truesdell, Anessaytoward a Unified 
TheoryofSpecialFunctions. (Annals of 
Mathematies Studies Number 18.) IV-+182S. 
Princeton 1948. Princeton University Press. Preis 


? 


Eleetromagnetie Theory (Proceedings 
of the second Symposium in Applied Mathematics 
of the American Mathematical Society). 91 8. 
mit 8 Abb. New York 1950. American Mathematical 
Society. Preis geb..3,— $ 


Dr. W. Weitzel, Lehrbuch der theoreti- 
schen Physik. Erster Band: Physik der Vor- 


_ gänge, Bewegung, Elektrizität, Licht, Wärme. XIV 


+ 7718. m. 270 Abb. Berlin-Göttingen-Heidelberg 
1949. Springer-Verlag. Preis geheftet 53,— DM, geb. 


56,90 DM. 


Veröffentlichungen des Institutes für Erdmessungen, 
Band VI, Teill u.2. Beiträge zur Lotab- 
weichungsausgleichung und Geoid- 
bestimmung. (Allgemeine Betrachtung und 
Darehführung im Bereich des Zentraleuropäischen 
Netzes.) 106S. m. 2Abb. und 25 Tafeln, Bamberg 
1949. Bamberger Verlagshaus Meisenbach & Co. Preis 
5,— DM. 1 


“Dr. se. techn. Edward R. van Driest, Dielineari- 
sierte Theoriederdreidimensionalen 
kompressiblen Unterschallströmung 
unddieexperimentelleUntersuchung 
von Rotationskörpern in einem ge- 
sehlossenen Windkanal. (Mitteilungen aus 


- dem Institut für Aerodynamik an der eidgen. Techn, 


Hochschule in Zürich. Herausgegeben von Prof, Dr. 
J. Ackeret. Nr. 16.) 31°S. m. 22 Abb. Zürich 1949. 
Verlag Leemann. Preis brosch. 7,— fr. 


- Dr. se. techn. William T. Sawyer, Experimen- 
talInvestigationofaStationaryas- 
cade of Aerodynamic Profifes. (Mit- 


 teilungen aus dem Institut für Aerodynamik an der 
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eidgen. Techn. Hochschule in Herausgegeben 
von Prof. Dr. J. Ackeret, .Nr..17.) 78 S. m. 45 Abb. 
Zürich 1949. Verlag Leemann, Präis brosch. 12,— fr. 


D:.-Ing. habil. Alf Pflüger (Prof. a. d. Techn. Hoch- 


schule Hannover), Stabilitätsprobleme 
der Elastostatik. VIIL + 339 S. m. 389 Abb. 
3erlin-Göttingen-Heidelberg 1950. Springer-Verlag, 


Preis geb. 34, 50 DM. 
Hans Jönek,Mechanik und Festi EEOERL, 
lehre., Fünfte Auflage. 2808. mit 322 Abb. 


Braunschweig 1950. Friedr. 
geb. 7,380 DM. 


Prof. Dr. R, Rothe, Höhere Mathematik 
für Mathematiker, Physiker, ’Inge- 
nieure, Teil V, Formelsammlung unter Mitwir- 
kung von Stud.-Rat Oscar Degosang und Dr.-Ing. 
Gerhard Dobbrack (Teubners Mathem. Leitfäden Bd. 
43). Zweite Auflage, 124 S. mit 74 Abb, Leipzig1949. 
B. G. Teubner Verlagsgesellschaft. Preis kart. 4,— DM. 


Festschrift zur 125-Jahrfeier des 
Physikalischen Vereins zu Frank- 
furt am Main 1949. Herausgegeben im Auftrage 


NACHRICHTEN 


Tagung der Gesellschaft für angewandte Mathema- 
tik und Mechanik vom 16. bis 19. April 1950 in 
Darmstadt 


In der Zeit vom 16. bis 19. April fand in Darmstadt 
. die diesjährige Tagung der GaMM statt. Der Besuch 
dieser Tagung war außerordentlich stark, Nicht we- 
niger als 240 "Mitglieder und- Nichtmitglieder der Ge- 
sellschaft aus dem In- und Auslande hatten sich ein- 
gefunden. Rund 60 Vorträge wurden gehalten, so daß 
es sich leider erforderlich machte, Parallelsitzungen 
einzurichten. Dieser Schönheitsfehler wurde dadurch 
gemildert, daß sich jeder an Hand der gedruckt vor- 
liege nden kurzen Vortragsauszüge die Vorträge aus- 
suchen konnte, die ihn besonders interessierten, und 
daß die Vortragszeiten so festgesetzt waren, daß ein 
Hinüberwechseln von der einen zur anderen Gruppe 
möglich war. 


Schon im Laufe des Sonntages hatten "sich eine 
Reihe von Tagungsteilnehmern eingefunden, die zu 
einem zwanglosen Begr üßungsabend in der Mensa der 
Hochschule zusammenkımen. In der Eröffnungssit- 

.zung am Montag morgen wurden die Erschienenen 
durch den Vorsitzenden der GaMM der amerika- 
nischen Besatzungszone Herr Grammel-Stutt- 

gart, Herrn Marguerre-Darmstadt in Vertretung 
des erkt ankten Rektors der Hoshschule und Herrn 
Metzger,dem Oberbürgermeister der Stadt Darm- 
stadt begrü üßt. Anschließend sprach Prof. Walther- 
Darmstadt übsr Stand und Probleme der angewandten 
Mathematik, und es wurde ein ausführlicher Bericht 
von Prof. Hartree-Cambridge, der leider am Er- 
scheinen verhindert war, über automatische Rechen- 
maschinen von H.-J. D r eyer- Darmstadt verlesen. 
Daran schlossen sich die Vorträge: 


C. Weber, Schlewecke: Randverformung einer 
rar mit Normalbelastung. 

Gross, Rom: Über eine neue Mntegrations- 
niet für lineare partielle Differ entialgleichungs- 
systeme und ihre Anwendung auf das Gleichgewichts- 
problem einer Brücke, 

C. Torre, Wien: Anwendung der zuyklogpapbı, 
schen Abbildung i in der Bruchtheorie. \ 


Am Nachmittag tagten nebeneinander die Sektion 
für Mechanik und die für Stochastik (Wahr- 
scheinlichkeitsrechnung und Mathematische Statistik). 
In der ersten wurden dıe folgenden Vorträge gehalten: 


Nachrichten 


Vieweg & Sohn. Preis. 


angewandte Mathematik und Mecha- 
“ deutsam, weilsich in ihr die beiden Gesellschaften der 3 3, 


2. es "Math. Ma 
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des Vorstandes des Praha Vereines von Dr. E 

Ludwig Protz. 478. z 4 

» Veröffentlichungen des Instituts für Erdmessung: N 

1ı— Das Zentraleuropäische Drei- . 
ecksnetz, Grundlagen. 1224.28, 

5 Tafeln. Bamberg 1949. - Bamberger Verlags- £ 

haus Meisenbach & Co. Preis brosch. 4,— DM. I 

3 — Professor Hans Boltz und sein & 

Werk. 77 8. Bamberg 1949. Bamberger Ver- £ 

lagshaus Meisenbach & Co. _Preis brosch. 3 

3,— DM. x 

4 — Beiträge zur Ausgleichungastro- \ 

nomisch- geodätischer Netze; 3 


unter : besonderer Berücksich.- 
tigung des Zentraleuropäischen 
Netzes. 143 S. mit 4 Abb. Bamberg 1949. 
nn Meisenbach & Co. Preis brosch. 
4,— DM 
Be zur Zentrierung: von ; 
Richtungen. 47 S.mit 6 Abb. Bamberg = _ = 
1949. Bamberger Verlagshaus Meisenbach & Co. 


Preis brosch. 3,— DM. _ - N 

F 

£ 

» ‚ > % b) 

W. Haacke, Braunschweig: Untersuchungen = 
über die Stabilisierung eines Dopp>lpendels mit perio- 2 
disch erschüttertem Aufhängepunkt. I 
H. Jung, Stuttgart: Über eine neue Iterations- N 
methode zur Berechnung der Knicklast, 3 
H. Münz, Braunschweig: Zur Berechnung 3 
schwach gewölbter, kreissymmetrischer Platten. R 
F. Weidenhammer, Clausthal: Dar einge-, 3 
spannte, axial-pulsierende, belastete Stab als Stabili- e 
tätsproblem. 3 
W. Lode, Leverkusen: Zur thermokinetischen u 
Theorie’ dar Elastomere., i -3 
M.Bergsträßer, Bingen: V-Verzahnung mit 2 
‚großen Profilverschiebungen. . 
O0. Heck, Wiesbaden: Beitrag zur analytischen 2 

= 


Behandlung des Kurbel _etriebes. 


In der Sektion Stochasti Krenfachen: 


ER Münzner, Göttingen: ‚Die Methoden der Ber 
athematischen Statistik. ag 

"ww. Ludwig, Heidelbarg: Biomathematik (au: f 
wendung ma shematischer, insbesondere stochastischer 
Methoden in der Biologie). 

K. Freudenberg, Berlin: Medizinische Sta- 
tistik. 

H. Kellerer, München: Ausg>wählte Probleme 
der Wirtschafts- und soziologischen Statistik. 

M. P.Geppert, Bad Nauheim: Mathematisch- 
statistische Fragen der Fabrikationskontrolle. . 


Am nächsten Vormittag fand zunächstdieH a u pt- 
versammlung der Gesellschaft für 


DO Le 


nik statt. Diese war deswegen außerordentlich be- 


englischen und amerikanischen Zone zu einer Gesell- 
schaft zusammenschlossen. Damit besteht. wieder die 
GaMM in ihrer alten Form. Da Prof. Prandtl, 

Göttingen, der langjährige Vorsitzende der Gesell- 
schaft, wegen seines hohen wre yernahme des. 
Vorsitzes ablehnte, wurde Pr ramm el, Stut 

gas auf drei Jahre zum Vorsitzenden und Pro 


randtl, Göttingen, z renv len au 
Leben gewählt, - Beide ins 
mig. 
Sodann wurder 
Grammel die 


die in allen wese 
a u 


 mungslehre: 


Be 
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genehmigt. Entsprechend diesen Statuten wurden 
Prof. Walther, Darmstadt, zum stellvertretenden 
Vorsitzenden, Prof. Görtler, Freiburg, zum 
Schriftführer einstimmig gewählt. Weiter erfolgte 
die Wahl des wissenschaftlichen. Beirates, der 5—7 
Mitglieder umfassen soll. Er setzt sich zusammen 
aus: Prof, Betz, Göttingen, Prof, L-Föppl, 
München, Prof. Pöschl, Karlsruhe, Prof. E. 
Sehmidt, Braunschweig, Prof. G. Schulz, 
Aachen, Prof. ©. Weber, Schlewecke und Prof. 
Fr.A. Willers, Dresden. Zu Kassenprüfern wur- 
den Prof. Sauer, München, und Prof. Vogel- 
pohl, Göttingen, gewählt. ; 

Die Mitgliederbeiträge wurden für das Jahr 1950 
und 1951 auf je 5,— DM festgesetzt; doch wurde dem 
Vorstand die Ermächtigung erteilt, diesen Beitrag für 
das Jahr 1951, falls es die Lage gestattet, herabzu- 
setzen und ihn in einzelnen Fällen zu ermäßigen oder 
erlassen. Nichtmitglieder der GaMM, die an dieser 
Tagung teilnahmen, wurden gebeten, einen entspre- 
chenden Betrag als Tagungsbeitrag einzuzahlen, So- 
dann wurde über den Ort der nächsten Tagung, die 
wieder im Frühjahr stattfinden soll, gesprochen, der 
Beschluß gefaßt, möglichst bald die Satzungen und 
ein Mitgliederverzeichnis drucken zu lassen ünd von 
den Preisermäßigungen Kenntnis gegeben, die den 
Mitgliedern der GaMM für den Bezug einzelner Zeit- 
schriften und Bücher zustehen. Die Gesellschaft 
tritt dem fachwissenschaftlichen Beirat des deutschen 
Forschungsrates bei. 

Nach Beendigung der Hauptversammlung nahm 
man wieder die wissenschaftliche Arbeit auf. Neben- 
einander tagten die Sektion für Elastokinetik, 
in der folgende Vorträge gehalten wurden: 

K. Ludwig, Hannover: Der Massenpunkt mit 
konstanter Horizontalgeschwindigkeit auf dem beider- 
seits unbegrenzten, elastisch gebetteten Geleis "mit 
Reibung. 

Th. Schade, Aachen: Neuartige Behandlung der 
Poissonschen und der inhomogenen Bipotentialglei- 
chung bei rechteckigen Bereichen mit Anwendung auf 
Probleme der Torsion in der Plattenbiegung. 

W. Zerna, Hannover: Allgemeine Theorie ela- 
stischer Flächenträger ohne Bernoullische Annahme. 

W. Swida, Karlsruhe: Plattengleichungen für 


„- den elastisch-plastischen Zustand. 


J. Weinhold, Holzgerlingen: Eine näherungs- 
weise Integration der Differentialgleichung y’’ = f (y) 
bei Druckbiegeproblemen. 

U. Wegner, Heidelberg: Variationsmethoden 
zur praktischen Lösung von ebenen und räumlichen 
Spannungsp-oblemen. 

0.Emersleben, Berlin: Über Berechnung der 
Gitterenergie endlicher Kristallstücke, 
und die Sektion für Mathematische Ma- 

schinen und Instrumente: 

W. Bader, Stuttgart: Auflösung von Polynom- 
gleichungen auf elektrischem Wege (mit Vorführung). 

H.-J. Dreyer, Darmstadt: Automatisches licht- 
elektrisches Kurvenabtasten bei Integrieranlagen. 

K. Zuse, Neukirchen: Programmgesteuerte Re- 

- chenmaschinen in Deutschland. 

Th. Fromme, Weil/Rhein: Die Schachmaschine. 
 K. Ramsayer, Stuttgart: Die Funktions- 
rechenmaschine. - 


Am Nachmittag sprachen in der Gruppe Strö- 


_W. Barth, Nonnenborn: Modellversuche mit 


- wasserdurchströmten Zyklonabscheidern zur Voraus- 
_ bestimmung der Abscheideleistung. 


M. Strscheletzky, Kreßbronn: Strömung 


im Übergangsraum der Wasserturbinen. 


_W. Richter, Frankenthal: Lösung einer Auf- 
gabe über Strömung in Leitungs-Systemen. 


Nachrichten - 
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M. Herbeck, Göttingen: Der Wärmeaustausch 
zwischen einem geheizten Band und einer Konvek- 
tionsströmung, & 

K. Nickel, Tübingen: Über ein Minimumpro- 
blem aus der Tragflügeltheorie. 

H. Calame, Karlsruhe: Der Einfluß von Staub- 
größe und Staubdichte auf die Druckverhältnisse bei 
beschleunigten Strömungen. 

N. Scholz, Braunschweig: Ein einfaches Sin- 
gularitätenverfahren zur Erzeugung von Schaufel- 
gittern. 

E. Pieruschka, Kassel: Die mathematischen 
Grundlagen zu einer Meßmethode des Schubmoduls 
zäher Flüssigkeiten. 

und in der Gruppe Mathematik: 

E. Ullrich, Gießen: Ausgezeichnete Betrag- 
flächen analytischer Funktionen. 

J.Heinhold, München: Schmiegungsverfahren 
der konformen Abbildung. 

F. Löbell, München: Eine reelle Deutung kom- 
p!exer Vektoren, 

H. Buchholz, Darmstadt: Die Entwicklung 
beschränkt willkürlicher Funktionen nach den Eigen- 
funktionen des Ringparabols. 

Anschließend wurde das Institut für Praktische 
Mathematik (Mathematische Instrumente, Unter- 
richtsmodelle, Lösungsmethoden linearer Gleichungs- 


systeme, Integrieranlage, Rechnen mit Dualzahlen) ° 


besichtigt. 

Am Mittwoch Vormittag setzte die Grupps für 
Strömungslehre ihre Sitzungen fort. Es 
sprachen: e a 

H. Görtler, Freiburg: Über die Lösungen nicht- 
linearer partieller Differentialgleichungen vom Rei- 
bungsschichttypus. 

K.Schröder, Berlin: Bemerkungen zum Diffe- 
renzenverfahren bei der Berechnung der laminaren 
Grenzschicht. ; 

H. Schubert, Rostock: Über ein gemischtes 
räumliches Randwertproblem der Potentialtheorie. 

H. Behrboom, Waldkirch: Zum Schalldurch- 
gang in zwei- und dreidimensionalen. Düsenströ- 
mungen. 

H. Richter, Haltingen: Die Stoßherzkurve bei 
allgemeiner Zustandsgleichung. 

Fr. Wecken, Haltingen: Expansion einer Gas- 
kugel hohen Druckes. 


M, Sehäfer, Göttingen: Über die Berechnung 


von kleinen stationären Störungen in einer gemischten 
Überschall-Unterschall-Strömung mittels einer ein- 
heitlichen Charakteristikenmethode. 


In der parallel stattfindenden Sitzung Prak- 
tische Analysis wurde zunächst ein von J. ©. 
P. Miller, London, angesagter Vortrag: Survey 
of Aids to Computation in deutscher Übersetzung von 
H. Unger, Darmstadt, verlesen. Es folgten die 
weiteren Vorträge: 

W.Schnell, Darmstadt: Interpolation im Kom- 
plexen, 

A.van Wijngaarden, Amsterdam: Grund- 
sätzliche Probleme der Abrundungsfehler. 

W., Quade, Karlsruhe: Grundsätzliches zur 
numerischen Integration von gewöhnlichen Differen- 
tialgleichungen. 

J. Weissinger, Hamburg: Zur Fehlerab- 


- schätzung bei der numerischen Integration (Differen- 


zenverfahren) gewöhnlicher Differentialgleichungen. 

E. Weinel, Jena: Bemerkungen zum Verfahren 
von Runge und Kutta, 

L. Collatz, Hannover: Über die Konvergenz- 
kriterien b>i Iterationsvorfahren für lineare Glei- 
chungssysteme, 

H. Sassenfeld, Darmstadt: Zur Iteration bei 


linearen Gleichungssystemen, 


H. Unger, Darmstadt: Nichtlineare Behandlung 
von Eigenwertaufgaben bei Matrizen. 
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J. N. Lehmann, Dresden: Zur Praxis der 
Eig genwertei ngr eNZUNg. 
R. Zurmühl, Darmstadt: Zum Grastfe-Ver- 


fahren bei komplexen Wurzeln. 


L.Sehwarz, Oberwolfach: Methoden zur Theo- 


rie der Integr altransformationen. 

Bei allen Tagungen sind die Aussprachen 
zwischen den einzelnen Teilnehmern mindestens 
ebanso wichtig wie die Vorträge und die anschlie- 
ßenden Diskussionen. Zu solchen Aussprachen war 
reichlich Gelegenheit gegeben. Eine Einladung der 
Stadt Darmstadt vereinte eine große Zahl der Teil- 
nehmer bei einem Empfang, der in wiederhergestellten 
auf der Mathildenhöhe am Mon- 
tag abend stattfand, bis spät in die Nacht, Am 
Die nstag abend fand eine Festvorstellung im Landes- 

{ theater: „Richard Wagner, Die Meistersinger von 
Nürnberg‘ statt, und am Mittwoch nachmittag wurde 
ein Austlug an die Bergstraße mit gemeinsamer 
Kaffeetafel in Seeheim und Abendessen in Jugen- 
heim gemacht. Die Obstbäume standen in voller 
Blüte und das Wetter erlaubte, daß der Kaffee im 
Freien eing:nommen werden konnte, Die noch am 
Donnerstag anwesenden Teilnehmer fuhren zur Be- 
sichtigung “des Werkes der Firma Adam Opel, Rüssel- 
heim oder der MAN, Gustavsburg und trafen sich 
am Nachmittag zur Besichtigung der Sektkellerei 
Burgeff in Hochheim. Diesen Firmen sei auch hier 
für die Einladung zu diesen außerordentlich inter- 
essierenden Besichtigungen herzlich gedankt. Zum 
Schluß wiederhole ich im Namen aller Teilnehmer 
den mehrfach ausgesprochenen Dank für die mühe- 
volle Vorbereitung und die gute Organisation, die 
von Herrn Walt her und seinen Mitarbeitern durch- 
geführt wurde, Die Tagung wurde infolge dieser aus- 

| : gezeichneten Vorbereitung” eine der schönsten Ta- 
= gungen der Gesellschaft für angewandte Mathematik 

' und Mechanik, 


Dresden. 3 
Felix Klein-Feier in Göttingen 
Die Reden, die von Kneser, Lietzmann 
Bi; und Prandtl bei der Felix- Klein-Feier gehalten 
% wurden, sindin einem Heft des Göttinger Universitäts- 
v bundes 'abgedr uckt. Der Preis beträgt einschl. Porto 
ER und Briefumschlag 0,65 DM. Das Mathematische In- 
stitut der Universität Göttingen schickt nach Ein- 


Willers. 


3 zahlung dieses Betrages das Heft zu. 

j Tagung der deutschen Mathematiker-Vereinigung 

er in Köln 

PR, -Bei dem Bsricht über die Tagung der deutschen 
- "Mathematiker-Vereinigung in Köln ist durch ein be- 
E dauerliches Versehen der, Vortrag von F, Löbell- 
Y München: „Einige: Ba griffsbildungen der Theorie der 
Flächenabbildungen‘“ unerwähnt geblieben. Kamke. 


Berlin: Der o. Prof. an der Techn. Universität 
Berlin, Dr. G.Hamel, wurde am]. Juli 1949 emeri- 
tiert. 


Dresden: Der persönliche o. Prof. an der Uni- 
versität Jena, Dr.K. Maruhn, 
fessor mit Lehrstuhlfür das Fach der reinen Mathema- 

tikan der Techn, Hochschule Dresden ernannt. 


Greifswald: Der Dozent, Dr. phil. R.Kochen- 
dörfer, wurde zum Professor mit vollem Lehrauf- 


Verantwortlich für den Inhalt: 


trag für das Gebiet der Mathematik a an Sr 
_ sität Greifswald ernannt. Se E 


“München für reine und angewandte Mathematik. 


wurde zum Pro- 
_ Ebenda, Zeile 1 


Prof. Dr. Fr. A. Willers, Dresden 
Verlag: Akademie-Verlag GmbH., Berlin, NW 7, Schiffbauerdamm 19, Fernspre 
konto: Berlin 35 021, Bestell- und Verlagsnummer dieses Heftes: 1009/30 | 
und Mechanik erscheint monatlich. Bezugspreis: vierteljährlich RM1 
Abbestellungen können nur bis 4 Wochen vor Quartalsende anerkannt v j 
noch geliefert. Verantwortlich für den Anzeigenteil: Curt F. ae 
Julius Beltz, Langensalza. 57273/4022 — Nr.417. — Veröffentli 
Militär-Verwaltung in Deutsch 


_ München: Studienrat Dr rer,nat. K.. Se Made ee 
bach habilitierte sich an der Techn. Hochschule er 


- Am 2, Februar 1950 verstarb im Alter von 76 Jahren 
der 0. Prof. an der Universität München, Dr. Con- 
stantin Carathöodory. 


Würzburg: Prof.Dr. 6. Aumann wurde zum u 
o. Prof. der Mathematik an der le Greifswald -. 
Arnann: E EB 


. 


Hamburg: Mitte März verstarb der asstars 
professor der Univ. Hamburg Dr. Paul Riebe- 
sell, Präsident der Hamburger Feuerkasse. a 


- Druckfehlerberichti gung x, 


P. Funk, Uber id die Berechnung a en 
kritischen Drehzahlen bei homoge- 
a und fasthomogenen Maschinen“ 

Z.angew. Math. Mech. 15 (1935) 113-120. 2 

In Fußnote 3 der genannten Arbeit liegt eine Ver | 
wechslung vor und außerdem kann der Text der Fuß- 7 
note zu einem Mißverständnis Anlaß geben. An Stelle 
des Satzes: 


„Bei dem Beweis beruft sich Grammel gelegentlich ä 
auf die Mitwirkung von O. Perron.‘“ «:; 


Soll es heißen: 


„Bei dem Beweis eines Hilfspatzes zu Bi von Herrn 
Grammel dort aufgestellten und bewiesenen Sätzen 
beruft sich Herr Grammel anläßlich des Beweises eines 
Hilfssatzes auf die Mitwirkung von Herrn Perron“. Ich FE 
sehe mich zu dieser Berichtigung um so mehr veranlaßt, FR 
da die Anlage der algebraischen Beweise ausschließlich 
von Herrn Grammel herrühren und weil dies bei der } 
ursprünglichen Abfassung des Textes nicht genügend N 
deutlich zum Ausdruck. kommt. TE 


‚Wien. iD; Funk 3 


Ne, Er I Ye 

Paul ER Herleitung einer ‘Näherunge . Ba 
formel von Laplace für die Binomial- 
verteilung ohne RN 2, Kerne 
Math. Mech. 29 (1949), S. 368— 374. ES 


$.373 unten. Der Satz „Ferner ist ... 5 Probleme 2 
ist zu ändern in ‚‚Ferner beseitige ich die Glieder 2 


+ Sim Exponenten durch eine "Parallelverschiebung 


ändere aber aus schreibtechnischen Grün nden die 
zeichnung acht." 


Ebenda letzte Zeile, „oz ist mi 
S. 314, Zeile 13, ebenso. - 


Ebenda. Zu der Ya be 
Ebenda, 4 Zeilen ; 
ist einzufügen. 350 


Brez: R 


+ Es entfällt er Re re 
2 im 7% % 


